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ABSTRACT 


ii. 


The  present  study  is  to  theoretically  investigate  the  effect  of  elastically 
restrained  ends  on  the  axial  and  transverse  vibrations  of  expansion  bellows. 
The  work  carried  out  is  relevant,  as  the  end  conditions  of  bellows  in  general  at 
industrial  sites  are  quite  complex  giving  rise  to  unequal  rotations  at  the  ends. 
Physically,  the  system  may  be  represented  as  a  short  pipe  nipple  welded  at 
both  ends  of  bellows  and  then  the  flanges  welded  to  the  other  end  of  pipe 
nipple 


The  bellows  are  modeled  as  an  equivalent  beam  including  rotatory 
inertia  as  has  been  considered  by  earlier  researchers  and  EJMA.  The  general 
method  followed  in  deriving  exact  frequency  expressions  in  the  present  work 
is  based  on  “Separation  of  Variables”  approach  in  order  to  find  the  axial  and 
transverse  natural  frequencies  of  single  bellows.  This  method  is  nearly  well 
adopted  in  analyzing  the  vibration  of  Universal  (double)  bellows  also. 

Surprisingly,  the  previous  works  contributed  by  researchers  Li  Ting-Xin 
et  al  and  Jakubauskas  V.F  et  al.  [4,5,15],  show  that  they  had  paid  more 
attention  to  study  the  effect  of  classical  fixed-fixed  type  of  boundary 
conditions  only  and  obtained  axial  and  transverse  natural  frequencies.  It  is 
also  seen  that  no  provision  has  been  made  to  these  effects  in  EJMA  code 
also  while  computing  the  axial  and  lateral  vibration  frequencies  for  single  and 
Universal  bellows  and  so  may  not  be  treated  as  accurate  estimates. 


The  influence  and  variation  of  the  rotational  elastic  restraint  parameter 
is  studied.  It  is  seen  that  the  effect  of  this  parameter  on  natural  frequencies  is 
quite  significant.  The  effect  of  variation  of  the  internal  pressure  and  flow 
velocity  of  fluid  in  the  bellows  is  also  studied.  The  results  obtained  thereon 
by  solving  the  closed  form  frequency  equations  are  then  confirmed 
numerically  by  using  the  finite  element  method  and  are  found  to  be  in  below 
1 0%  of  error  for  all  the  cases. 

The  particular  case  of  bellows  with  fixed-fixed  boundary  condition  is  a 
reduced  case  wherein  the  elastic  restraints  against  rotation  are  compared 
well  with  the  experimental,  numerical  and  theoretical  results  for  both  axial  and 
transverse  vibrations  respectively.  All  the  results  are  thereon  plotted  for 
variation  of  rotational  stiffness  internal  pressure  and  flow  velocity  versus  axial 
and  transverse  frequencies. 
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y/D  at  T=oc,  Ms=0,  N=1 

123 

7.6 

y/D  at  T=oc,  Ms=0,  N=2 

123 

7.7 

y/D  at  T=oc,  Ms=0,  T=0.01  to  100 

125 
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TABLE  NO. 

DESCRIPTION 

PAGE 

8.1 

Rocking  Mode  Natural  Frequency  of  DBEJ  for  T=oc, 

N=1  and  Ms  =0 

137 

8.2 

Rocking  Mode  Natural  Frequency  of  DBEJ  for  T=oc, 

N=1  and  ks  =0 

138 

8.3 

Comparison  of  Rocking  Natural  Frequencies  of  DBEJ 
at  T=oc 

138 

9.1 

Comparison  of  Frequencies  at  T=cc 

155 

9.2 

Fundamental  Frequency  for  various  values  of  T=0.01  to  cc 

156 

9.3 

Response  of  bellows  to  Fateral  Seismic  Excitation 

157 
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X. 

Nomenclature 

Ac  :  cross  sectional  metal  area  of  one  bellow  convolution,  m2 

A0  :  Bellows  effective  area  corresponding  to  mean  diameter  of  bellow 

convolution,  m2 

Cd,  Cf,  Cp  :  factor  used  in  design  calculations  to  relate  U  shaped  bellows 
convolution  segment 

Cm  :  Material  strength  factor 

Ct  :  temperature  correction  factor  for  bellows  fatigue  life  below  creep 

range 

Cwb  :  longitudinal  weld  joint  efficiency  factor 

Cz  :  transition  point  factor 

Db  :  Inside  diameter  of  cylindrical  tangent  and  bellows  convolution,  m 

Dm  :  Mean  diameter  of  bellows  convolution,  m 

Eb  :  Modulus  of  elasticity  of  bellows  material  at  design  temperature, 

N/m2 

Ebr  :  Modulus  of  elasticity  of  bellows  material  at  room  temperature,  N/m2 

e  :  Total  equivalent  axial  movement  per  convolution,  mm 

fiu  :  Bellows  theoretical  axial  elastic  spring  rate  per  convolution,  N/m 

f  :  frequency,  Hz 

h  :  convolution  depth  less  bellows  material  thickness,  in 

J  :  mass  moment  of  inertia  per  unit  length,  kg/m 

k  :  factor  which  considers  the  stiffening  effect  of  attachment  weld  and 

the  end  convolution  on  pressure  capacity  of  bellows  tangent 

Lb  :  bellows  convoluted  length,  m 


:  bellows  tangent  collar  length,  m 


:  bellows  tangent  length,  m 

total  mass  of  bellows  per  unit  length  includes  mass  of  bellows  and 
fluid  mass 

:  Support  mass  of  supports,  kg 

:  mode  number 

:  number  of  plies  of  thickness  t 
:  number  of  convolutions 
:  fatigue  life,  number  of  cycles  to  failure,  cycles 
:  internal  pressure  force,  MPa 
:  convolution  pitch,  mm 
:  Rotational  Stiffness 
:  mean  radius  of  bellows,  m 

:  allowable  bellows  material  stress  at  design  temperature,  N/m2 
:  allowable  collar  material  stress  at  design  temperature,  N/m2 
:  total  stress,  N/m2 

:  yield  stress  at  design  temperature,  N/m2 
:  bellows  thickness  of  one  ply,  mm 

:  bellows  material  thickness  for  one  ply,  corrected  for  thinning,  mm 
:  bending  stiffness,  Nm2 
:  time 

:  Rotational  non-dimensional  parameter  (Ti,  T2) 

:  velocity  of  flow,  m/s 


w,  y  :  deflection,  mm 

x,  z  :  axial  coordinate 

pb  :  mass  density  of  bellows  material,  kg/cu.m, 

pf  :  mass  density  of  fluid  flowing  through  bellows,  kg/cu. 

co  :  Frequency  in  radians/s 


Note:  Variables  have  also  been  defined  wherever  required  in  chapters 
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CHAPTER  1 
INTRODUCTION 

1.1  Preamble 

Generally,  U-shaped  bellows  expansion  joints  as  shown  in  Fig  1.1,  are  used  as 
flexible  elements  in  piping  systems  to  absorb  the  relative  movements  resulting  from 
thermal  expansion  in  pipe  system.  In  order  to  prevent  excessive  stress  and  avoid  any 
damage  to  either  the  equipment  or  pipeline  expansion  bellows  are  employed. 
However,  it  is  seen  that  bellows  fail  to  give  many  years  of  satisfactory  performance 
when  they  are  not  properly  designed  for  specified  piping  system  conditions. 


They  are  used  in  chemical  plants,  nuclear  power  systems,  heating  and  cooling 
systems  and  cryogenic  plants.  Typical  service  conditions  could  be  as  severe  as 
pressures  ranging  from  full  vacuum  to  lOOOpsig  and  temperatures  from  minus 
420deg  F  to  1800deg  F.  Therefore  the  design  and  fabrication  of  an  expansion  bellow 
requires  serious  attention  as  it  falls  into  the  category  of  a  highly  engineered  product. 


Failures  of  bellows  can  occur  because  of  many  reasons,  but  experience  has 
shown  that  bellows  failure  due  to  mechanical  or  flow  induced  vibration  falls  into 
fairly  distinct  category.  Hence,  the  design  of  bellows  for  vibrations  has  to  be  viewed 
seriously.  For  example,  the  reciprocating  machines  such  as  compressors  are  usually 
connected  to  piping  system  through  flexible  elements  like  expansion  bellows  in  order 
to  prevent  noise  and  stress  propagation  over  the  whole  piping  system.  The  bellows  in 
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such  piping  systems  are  subjected  to  severe  vibrations.  In  piping  systems  designed  for 
the  thermal  and  nuclear  power  plants  located  in  earthquake  prone  regions,  the 
expansion  bellows  are  to  be  designed  to  take  care  of  the  seismic  excitations  also. 
Non-availability  of  sufficient  space  to  route  the  piping  is  commonly  envisaged  in  a 
chemical  plant.  Expansion  bellows  employed  in  such  piping  systems  are  subjected  to 
axial  movement,  angular  rotation,  lateral  deflection  or  any  combination  of  these, 
giving  rise  to  forces  and  moments.  The  influence  of  such  forces  and  moments  shall  be 
considered  in  determining  the  axial  and  transverse  natural  frequencies  of  bellows 
while  carrying  out  dynamic  analysis. 

The  most  important  part  of  the  expansion  bellow  is  the  convoluted  portion. 
The  bellows  convolutions  are  available  in  variety  of  shapes,  the  most  popular  of 
which  consists  of  flat  circular  rings  connected  to  two  toroidal  half  rings,  forming  the 
root  and  tip  of  the  convolution.  These  are  called  as  U-convolutions.  This  shape  is 
most  widely  used  because  it  permits  formation  of  complete  bellows  either 
mechanically  or  hydraulically  from  a  single  piece  of  cylinder  and  this  thesis  is 
restricted  to  the  study  of  such  U-shaped  bellows  only. 

1.2  The  Bellows 

The  bellows  are  the  convolutions  and  flexible  portion  of  the  expansion  joint. 
They  must  be  strong  enough  circumferentially  to  withstand  the  pressure  and  flexible 
enough  longitudinally  to  accept  the  deflections  for  which  they  are  designed.  Hence, 


25 


the  strength  with  flexibility  is  a  unique  design  problem  that  is  not  often  found  in 
other  components  in  industrial  equipment. 

Most  of  the  structures  are  designed  to  inhibit  deflection  when  acted  upon  by 
outside  forces.  Since  the  bellows  must  accept  deflections  repetitively,  the  deflections 
will  give  rise  to  stresses.  These  stresses  must  be  kept  as  low  as  possible  so  that  the 
repeated  deflections  will  not  result  in  premature  fatigue  failures.  Reduction  of 
bending  stress  for  a  given  deflection  is  easily  achieved  by  simply  reducing  the 
thickness  of  convolution.  However,  in  order  to  withstand  the  pressure,  the 
convolution  must  have  a  thickness  great  enough  that  the  pressure  induced  membrane 
stress  is  equal  to  or  less  than  the  allowable  stress  level  of  the  material  chosen  at  the 
design  temperature.  This  conflicting  need,  thickness  for  pressure  and  thinness  for 
flexibility  is  the  unique  design  problem  faced  by  the  designer. 

Bellows  are  normally  categorized  as  -  formed  membrane  (single-ply  un- 
reinforced  /  reinforced  &  multi-ply)  and  formed  head  (only  flanged  or  flanged- 
flued).  This  kind  of  categorization  is  based  on  the  type  of  numerical  method  used 
and  generally  applicable  for  an  elastic  and  small  deflection  analysis.  Single-ply  un- 
reinforced  type  bellows  are  generally  evaluated  using  linear  thin  shell  theory.  On  the 
other  hand,  exact  analysis  of  a  multi-ply  bellow  or  a  ring-reinforced  type  of  bellow 
requires  consideration  of  non-linear  interaction  between  the  plies,  or  between  the 
ring  and  plies  respectively.  The  flanged  and  flued  type  of  expansion  joint  is 
characterized  by  large  wall  thickness  and  is  governed  by  thick  wall  equation  [10]. 
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A  summary  of  the  major  differences  between  various  types  of  bellows  is 
illustrated  in  Table  1.1  [5]. 

Table  1.1  Differences  between  bellow  configurations 


Type  of  Bellow 

Linear  thin  shell 
theory 

Thick  wall 
behavior 

Nonlinear 

interaction 

Sharp  Corners 

Typical  use 

U-shaped 

Yes 

No 

No 

No 

Most  common 

Semi -toroidal 

Yes 

No 

No 

No 

High  pressure 

S -shaped 

Yes 

No 

No 

No 

As  U-shaped 

Omega 

Yes 

No 

No 

No 

High  pressure 

Multi-ply 

No 

No 

Yes 

No 

Low  stiffness 

Reinforced 

No 

No 

Yes 

No 

High  pressure 

Flanged  &  flued 

No 

Yes 

No 

No 

Exchanger  shell 

1.3  Formed  Membrane 

The  formed  membrane  type  is  popularly  known  as  the  bellows  type.  It  is  made 
by  forming  convolutions  from  a  thin  (a  few  thousands  of  an  inch)  shell  of  a  corrosion 
resistant  material,  such  as  austenitic  steel  or  a  nickel  base  alloy.  Figures  1.1  and  1.2 
show  typical  views  of  solid  model  and  cross-section  of  the  un-reinforced  that  are  used 
in  low-pressure  applications. 


Fig  1.1  Solid  Model  View  of  Bellows 
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Reinforcing  rings  are  added  when  instability  or  Squirm  of  the  bellows  is  of 
concern.  All  bellows  have  a  critical  pressure  at  which  they  become  unstable. 
Instability  occurs  in  two  modes  -  column  instability  (or  squirm),  and  in-plane 
deformation  of  the  convolution  side  wall.  Squirm  is  the  phenomena  whereby  the 
centerline  of  a  straight  bellows  develops  a  sideways  or  lateral  bow. 


Db 


Oja 

Fig  1.2  Cross-sectional  views  of  Bellows 
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1.4  Basis  of  Design  and  Theory 

Differential  longitudinal  expansion  between  the  shell  and  the  tube  in  a  heat 
exchanger  or  in  a  pipeline  is  a  problem  that  occurs  either  because  of  temperature  or 
pressure.  The  effect  of  pressure  induced  differential  growth  is  often  over-looked  in 
the  design  of  fixed  tube  sheet  heat  exchangers.  Physically,  it  means  that  there  is  a 
mismatch  in  the  axial  deformation  of  shell  and  tube  bundle  that  is  caused  by  the 
difference  in  the  state  of  their  pressure  loading. 

1.5  Type  of  Deflections 

In  order  to  properly  apply  expansion  joints  to  piping  systems  it  is  necessary  to 
visualize  the  deflections  that  result  from  thermal  expansion  or  movements  and 
vibrations  of  equipment  and  structures.  All  expansion  joints  do  not  accept  the  same 
types  of  deflection.  Many  can  accept  certain  loads  and  moments,  while  others  are 
incapable  of  resisting  externally  applied  forces.  Hence,  understanding  the  type, 
magnitude  and  direction  of  forces  and  deflections  is  critical  from  safety  point  of 
view. 

1.5.1  Axial  Deflection 

Axial  refers  to  being  parallel  to  the  centerline  of  the  expansion  joint. 
COMPRESSION  is  the  axial  deflection  that  will  shorten  the  bellows  length.  Often, 
thermal  expansion  in  the  piping  will  cause  the  expansion  joint  to  be  compressed. 
EXTENSION  is  the  axial  deflection  that  stretches  the  expansion  joint.  Piping  that 
operates  at  temperatures  lower  than  ambient  such  as  in  cryogenic  systems  -will 
contract,  causing  the  expansion  joint  to  stretch  as  shown  in  Fig  1.3. 
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I 

COMPRESSED  FREE  EXTENDED 

Fig  1.3  Axial  deflection 

1.5.2  Lateral  Deflection 

Lateral  or  Transverse  refers  to  the  direction  perpendicular  to  the  centerline  of 
the  expansion  joint  and  shown  in  Fig  1.4.  This  movement  occurs  with  both  the  ends 
of  the  expansion  joint  remaining  parallel  to  each  other,  with  their  centerlines  being 
displaced.  It  is  not  uncommon  to  find  different  lateral  deflections  can  occur  in  more 
than  one  plane.  As  the  expansion  joint  is  circular  in  cross-section,  all  the  various 
deflections  must  be  resolved  into  a  single  resultant  lateral  deflection  -which  is 
square  root  of  the  sum  of  the  squares  of  the  individual  deflections. 


Fig  1.4  Lateral  deflection 

1.5.3  Angular  Deflection 

When  an  expansion  joint  experiences  bending  about  its  center,  which  is  on 
the  centerline  and  halfway  between  the  ends  of  the  bellows,  the  deflection  is  referred 
to  as  angular  deflection  and  is  shown  in  Fig  1.5.  It  can  occur  in  any  plane  that  passes 
through  the  centerline  and  is  the  maximum  of  the  various  deflections,  and  not  the 


vector  sum  as  in  the  lateral  case. 


Fig  1.5  Angular  deflection 


1.5.4  Torsional  Deflection 

Figure  1.6  depicts  torsion  deflection  of  bellows.  Torsion  refers  to  twisting  one 
end  of  the  bellows  with  respect  to  the  other  end,  about  the  bellows  centerline.  The 
expansion  joints  are  normally  not  expected  to  accept  torsion  deflection. 


Fig  1.6  Torsional  deflection 


1.5.5  Cyclic  Deflections  and  Cyclic  Life 

Most  deflections  are  repeated  a  number  of  times  during  the  life  of  the  piping 
system,  since  the  deflections  usually  are  produced  by  change  in  temperature  that 
occur  each  time  the  system  is  started  and  stopped.  Each  time  a  deflection  occurs  it  is 
a  cycle.  The  number  of  cycles  is  important  to  assure  the  proper  design  of  an 
expansion  joint,  since  each  design  has  a  finite  predictable  life. 


Vibrations  cause  repetitive  deflections  and  can  cause  a  premature  failure  of 
expansion  joint.  Even  though  these  deflections  are  small  in  magnitude,  they  usually 
accumulate  huge  number  of  cycles  in  a  short  period  of  time.  Since  bellows  are 


31 


metallic  structures,  they  have  a  specified  resonant  natural  frequency  and  when 
driven  by  outside  vibrations  of  same  frequencies,  they  can  magnify  the  deflections 
until  they  exceed  the  yield  strength  of  the  bellows  material  and  induce  early  fatigue 
failure.  For  example,  when  a  piping  system  is  known  to  have  equipment  that  can 
produce  vibrations  such  as  pumps,  compressors  and  other  motor  or  turbine  driven 
devices  whose  rotational  speeds  or  frequencies  are  of  importance. 

1.6  Types  of  Expansion  J oints 

There  are  different  types  of  expansion  joints  -  each  has  a  definite  advantage 
and  limitation  and  a  unique  design  style  of  its  own.  The  thesis  cover  only  two 
common  types  of  expansion  joints  that  are  as  follows  - 

1.6.1  Single  Bellows  Expansion  Joint 

The  single  bellows  expansion  joint  as  shown  in  Fig  1.7,  has  a  simple  bellows 
element  with  end  connections.  It  will  deflect  in  any  direction  or  plane,  but  requires 
the  piping  to  be  controlled  to  the  direction  of  the  movement.  This  expansion  joint 
will  not  resist  any  deflections  with  any  force  other  than  the  resistance  of  bellows, 
which  is  a  function  of  the  spring  rate  and  deflection. 


_jm 
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Fig  1.7  Single  Bellows 
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1.6.2  Universal  OR  Double  Bellows  Expansion  Joint 

The  Universal  expansion  joint  as  shown  in  Fig  1.8,  consists  of  two  bellows  separated 
by  a  pipe  section  or  spool.  The  primary  purpose  of  this  kind  of  an  arrangement  is  to 
have  a  unit  that  will  accept  large  amounts  of  lateral  deflection.  The  amount  of  lateral 
deflection  that  it  can  accept  is  a  function  of  the  amount  of  angulation  of  each  bellows 
can  absorb  and  the  distance  between  the  bellows.  So,  for  a  given  bellows  element, 
the  amount  of  lateral  deflection  capability  can  be  increased  or  decreased  by  simply 
changing  the  length  of  the  center  spool. 

The  design  analysis  of  an  expansion  bellow  is  quite  complex  because  it  is 
stressed  due  to  operating  pressure,  axial  deflection  and  at  the  same  time  has  to 


remain  stable  and  work  up  to  the  designed  fatigue  life.  Any  assumed  bellow 
geometry  which  is  defined  by  diameter,  depth  of  convolution,  pitch,  number  of 
convolutions,  thickness  and  number  of  plies  for  a  given  material  of  construction  and 
manufacturing  technique  employed  etc.,  has  to  satisfy  all  the  above  mentioned 
criteria.  All  these  factors  are  interdependent  and  change  in  any  one  of  the  variables 


affect  the  overall  bellow  design. 
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1.7  Literature  Survey  and  Review  of  Vibrations  of  Fluid  Filled  Bellows 

Literature  survey  was  carried  out  to  review  the  present  state-of-the-art  in  the 
area  of  design  of  metallic  bellows  and  fluid-filled  bellow  vibrations.  The  aim  is  to 
identify  the  direction  required  for  further  theoretical  investigations  to  be  carried  out, 
in  order  to  present  a  fairly  good  and  acceptable  theoretical  model  for  predicting  the 
dynamic  behavior  of  fluid-  filled  bellows. 

A  thorough  survey  of  literature  reveals  that  most  of  the  early  investigators  on 
bellows  have  studied  at  great  depth  the  static  response  of  bellows  -  stress  patterns, 
design  aspects  of  bellows  and  the  bellow  failure  modes  like  root -bulge  or  collapse, 
evaluation  of  creep  rupture  and  creep-fatigue  and  stability  analysis.  Quite  often 
bellows  are  also  used  to  arrest  the  propagation  of  vibrations  in  the  piping  systems. 
Being  very  flexible,  bellows  are  susceptible  to  vibrations  and  more  so  can  be  easily 
excited  by  internal  flow. 

It  was  only  in  the  early  1990,  when  some  sporadic  failures  of  bellows  were 
reported  due  to  high  flow  velocities  of  internal  fluid  -  the  failure  of  bellows  due  to 
Flow  Induced  Vibrations  (FIV)  came  to  the  fore.  The  significance  of  dynamic 
response  on  bellows  reliability  did  not  attract  immediate  attention,  until  bellow 
failures  became  a  common  place  particularly  in  power  industry.  The  last  decade  has 
therefore  witnessed  a  vigorous  research  activity  in  the  field  of  flow  induced  vibration 
of  bellows  -  to  understand  the  phenomena  [2].  The  goal  herein  is  to  present  the  early 
developments  on  the  subject  and  to  attempt  a  physical  insight  into  the  vibration 
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problem  and  summarize  the  various  methods  proposed  in  the  literature-  theoretical 
and  numerical  and  present  further  direction  on  flow  induced  vibrations  in  bellows. 

1.7.1  General  Information  about  Flow  Induced  Vibrations 

Flow  induced  vibrations  (FIV)  of  elastic  bodies  involve  the  interaction 
between  the  elastic  and  inertia  forces  of  the  elastic  body  and  the  fluid.  It  is  seen  that 
the  vibration  excitation  mechanisms  are  not  well  understood,  since  the  FIV 
phenomena  is  complex  and  diverse.  The  exact  determination  of  the  nature  of  the 
interaction  between  the  structure  and  fluid  and  the  magnitude  of  force  of  the 
interaction  is  extremely  difficult.  The  FIV  appears  in  tubes  of  heat  exchangers  and 
steam  generators  -investigated  by  Paidoussis  (1979),  Shin  and  Wambsgnass  (1975); 
in  nuclear  fuel  assemblies,  Oldaker  et  al.  (1973);  in  hydro  technical  components  and 
structures,  Weaver  (1976);  in  naval  and  space  industry,  Sainsbury  and  King  (1971). 
The  wind  excitation  of  buildings  and  aboveground  structures  is  another  area  where 
FIV  is  applied  [13]. 

The  list  of  published  works  in  FIV  is  increasing  rapidly  because  of  the  diversity  of 
the  subject.  Naudascher  (1967),  Toebes  (1965)  and  others  attempted  to  classify  the 
flow  induce  vibrations.  However,  according  to  Weaver  (1989)  FIV  is  characterized  as 
follows- 

❖  Forced  vibrations 

Self-controlled  vibrations 
♦>  Self-excited  vibrations 
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Figures  1.9  {(a),  (b),  (c)  &  (d)}  represent  the  various  modes  of  vibration  of  bellows. 


b)  Nonsymmetric  longitudinal  mode 


c)  Local  convolute^bending  mode 
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Forced  vibrations  arise  because  of  various  sources.  Some  of  them  include  the 
vibrations  due  to  turbulent  flow,  vibrations  as  a  response  of  tall  buildings  or  aircraft 
structures  to  wind  gusts,  vibrations  of  ship  propeller  blades  excited  due  to  the 
periodical  flow  near  the  ship  hull,  and  the  response  of  a  pipe  carrying  fluid.  It  is 
believed  that  whether  the  excitation  phenomena  are  random  or  periodic,  the  motion  of 
the  structure  has  no  feed  back  effect  on  the  fluid  forces.  Therefore  the  excitation  force 
can  be  studied  separately  from  the  vibration  problem  using  a  rigid  model.  So  this 
separation  of  the  whole  response  problem  into  two  independent  ones  greatly 
simplifies  the  solution. 

In  self-excited  vibration  problems  some  periodicity  in  the  flow  exists  in  case 
of  stationary  structure.  When  this  periodicity  coincides  with  the  natural  frequency  of 
the  structure  resonance  takes  place.  The  vibration  amplitude  increases  until  the 
structural  motion  starts  to  control  the  fluid  excitation  force.  Under,  these  conditions  in 
some  fluid  velocity  range,  the  vibration  response  of  the  system  is  controlled  not  by 
fluid  velocity,  but  by  vibrating  structure.  The  fluid  velocity  range  is  the  lock-in 
region.  A  common  example  of  periodicity  in  the  flow  is  the  vortex  shedding, 
vibrations  of  stacks  /towers  or  turbine  blade.  Such  vibrations  can  be  controlled  or 
prevented  by  changing  the  stiffness  of  the  structure  in  order  to  alter  the  natural 
frequency  or  more  effectively  to  change  the  geometry  of  structure  to  alter  the  fluid 


excitation. 
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Self-excited  vibrations  appear  in  such  systems  when  the  motion  of  the 
structure  itself  creates  the  periodic  fluid  force  that  in  turn  amplifies  the  vibration  of 
structure.  The  periodic  force  doesn’t  exist  in  the  absence  of  structural  motion,  as  in 
case  of  self-controlled  vibrations  and  this  is  the  main  distinction  between  the  two 
types.  Examples  of  self-excited  vibrations  are  bending-torsion  vibrations  of  aircraft 
wings,  the  oscillation  of  gate  seals  and  the  vibrations  of  vertical  lift  gates.  Because  of 
the  interaction  of  fluid  and  elastic  forces,  both  self-controlled  and  self-excited 
vibrations  are  called  fluid-elastic  vibrations. 

In  FIV  of  liquid  flows,  it  is  important  to  take  into  account  the  inertia  of  the 
fluid  because  this  is  usually  great  enough  to  change  the  vibration  frequencies 
considerably  and  the  modes  of  vibration  of  the  system.  The  motion  of  any  rigid  or 
elastic  body  in  a  fluid  is  accompanied  by  flow  of  fluid  around  the  body.  Considering 
the  fluid  to  be  perfect  and  incompressible  and  its  flow  as  steady  and  irrotational,  both 
the  drag  and  lift  forces  for  the  body  possess  symmetry  in  center  according  to  D’ 
Alembert  and  equal  to  zero.  However,  if  the  motion  is  not  uniform,  the  flow  is  not 
steady.  In  this  case  the  flow  generates  a  drag  force  on  the  body  and  this  force  appears 
as  if  the  inertia  of  the  body  has  been  increased.  Since  the  vibration  motion  is  non- 
uniform,  the  mass  of  the  vibrating  body  appears  to  have  increased  by  an  “Additional 
mass”.  This  phenomenon  is  significant  because  it  can  lead  to  a  drop  in  the  natural 
frequency  of  the  vibrating  body,  if  the  fluid  density  is  sufficiently  high.  The  first 
attempt  to  evaluate  theoretically  the  additional  mass  concept  in  bellows  was  made  by 


Gerlach  (1969). 
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1.7.2  Early  Developments 

♦♦♦  Theoretical  Methods 

Love  [3]  was  the  first  to  propose  the  theory  of  elastic  shells  -  which  resulted 
in  two  formulations  for  axi-symmetric  shells-  Love-Meissner  and  E. Meissner 
formulations.  The  earliest  stress/deflection  analysis  of  bellows  relied  mainly  on 
simple  beam  [11]  or  beam  and  cylinder  approximations.  It  is  observed  that  these  gave 
only  qualitative  but  not  quantitative  assessments  of  bellows  stress/deflection.  The 
analytical  solutions  were  based  on  energy  methods  and  asymptotic  integration  - 
considered  only  axial  loading.  However,  a  number  of  investigators  later  used  energy 
methods  to  evaluate  the  response  of  bellows  to  axial  force  /  deflection.  Salzman 
evaluated  the  pressure  stresses  of  a  U-shaped  bellow  using  a  three-term  series  to 
represent  the  meridoinal  bending  deformation.  Dahl  [9]  used  a  four-term  series  for  an 
Omega  shaped  bellows.  Turner  and  Ford  [44]  used  a  five-term  series  to  represent 
circumferential  stress  and  evaluate  semi-toroidal  bellows.  Turner  extended  this  work 
to  U-shaped  bellows  -  evaluated  the  pressure  loading  by  approximating  the  bellows 
to  beam  theory.  Laupa  and  Weil  [39]  evaluated  U-shaped  and  semi-toroidal  bellows 
using  a  five  term  series  and  expressed  the  bending  deformation  of  a  toroid  and  a 
circular  plate  to  represent  the  convolution  sidewall.  They  repeated  the  same  work 
using  shell  theory.  Clark  [8]  developed  solutions  for  axial  deflection  of  Omega  type 
bellows  using  asymptotic  integration. 

Anderson  [1,2]  used  the  asymptotic  solution  given  by  Clark,  to  develop 
solutions  for  deflection  stresses  in  U-shaped  bellows.  He  also  introduced  corrections 
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to  compensate  for  the  approximations  suggested  by  Clark.  Pressure  stress  was  also 
evaluated  using  the  beam  theory  and  correction  factors  were  introduced  as  charts  in 
order  to  relate  the  simplified  equations  to  the  shell  behavior.  This  work  by  Anderson 
later  formed  the  basis  for  the  stress  equations  for  un-reinforced  /  reinforced  bellows  in 
the  Expansion  Joint  Manufacturers  Association  Handbook  [EJMA  10].  The  state-of- 
the-art  in  bellow  analysis  was  clearly  presented  by  Becht  IV  [5]  and  given  in  Table 
1.2. 

Tablel.2  State-of-the-art  in  Bellow  Analysis 


Type  of 
Bellow 

Elastic 

Elastic/Plast 

ic 

Creep 

External 

pressure 

Internal 

Pressure 

U-shaped 

a 

a 

a 

a 

m 

Multi-ply 

□ 

□ 

□ 

E 

E 

Reinforced 

0 

□ 

□ 

0 

0 

Flanged  & 
flued 

0 

□ 

□ 

0 

0 

S  Simplified  methods  of  evaluation  are  available 
a  Rigorous  numerical  analysis  performed 

□  Has  not  been  evaluated 


❖  Numerical  Methods 

Numerical  methods  gained  importance  only  with  the  advent  of  digital 
computation.  These  methods  are  the  finite  difference,  numerical  integration  and  finite 
element.  Krauss  [35]  solved  the  shell  equations  using  the  first  two  numerical 
techniques.  The  finite  element  method  involves  partitioning  of  the  shell  into  finite 


elements  to  get  more  exact  approximations. 
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♦>  Finite  Difference 

In  this  method  the  second  order  differential  equations  be  replaced  by  a  set  of 
linear  algebraic  difference  equations.  Their  derivatives  are  written  over  a  finite  length 
instead  of  infinitesimal  length  and  the  check  for  convergence  on  exact  solution  is 
attempted  by  solving  the  problem  with  two  different  mesh  densities.  The  finite 
difference  energy  method  is  an  improvement  on  the  finite  difference  method  -  uses 
the  principle  of  minimum  potential  energy. 

❖  Numerical  Integration 

In  this  method  the  second  order  differential  equations  be  converted  to  first 
order  differential  equations  by  using  the  finite  difference  equations.  Initial  values  are 
assumed  for  unknown  boundary  conditions,  and  Runge-Kutta  method  integrates 
equations  across  the  shell. 

❖  Finite  Element 

It  uses  the  element  stiffness  to  get  an  overall  stiffness  matrix.  The  shell  is 
approximated  by  a  number  of  finite  elements  -  having  known  displacement 
polynomials  or  Shape  Functions.  The  stiffness  of  each  element  is  found  out  by  energy 
methods  -  either  by  closed  form  solutions  or  by  numerical  integration. 


(F)  =  (K).  (8) 


(1.1) 
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(F)  Column  matrix  of  forces  at  node  points,  (K)  is  stiffness  matrix  and  (8)  column 
matrix  of  nodal  displacements. 

1.7.3  Flow  Induced  Vibration  Analysis 

The  first  detailed  study  of  flow  induced  vibrations  in  bellows  was  conducted  by 
Gerlach  [13].  He  observed  that  the  source  of  fluid  excitation  was  due  to  vortex 
shedding  from  the  convolution  tips  and  presented  a  simplified  method  of  calculating 
the  natural  frequencies  of  bellows  under  flow  induced  vibrations.  He  argued  that  a 
bellows  with  ‘N’  convolutions  could  be  represented  as  a  system  of  (2N-1)  identical 
masses  connected  by  ‘2N’  identical  springs.  Relatively,  simple  expressions  were 
given  for  elemental  masses  and  springs  in  terms  of  bellows  geometry.  He  assumed 
that  the  elemental  fluid  mass  could  be  considered  as  one-half  the  mass  of  the  fluid 
contained  between  bellows  convolutions  for  the  lowest  modes  of  vibration.  Gerlach 
derived  an  expression  for  calculating  fluid  added  mass  for  higher  modes  of  vibration 
in  which  the  fluid  is  pressed  out  and  sucked  in  between  convolutions.  Then  the 
natural  frequencies  of  axial  vibrations  were  determined  by  calculating  the  elemental 
mass  including  the  added  mass  and  elemental  stiffness  for  a  given  bellow. 

In  a  subsequent  paper,  Gerlach  [14]  developed  a  concept  of  ‘Stress  Indicator’-  based 
on  the  assumptions  of  linear  forced  vibration  theory  to  provide  an  index  of  the 
severity  of  vibrations.  Gerlach  also  suggested  that  there  was  no  ‘vortex  shedding’ 
occurring  in  the  absence  of  bellows  vibration  and  compute  Strouhal  number  by 


equation  (2.2) 
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00. 

Vp 

(1.2) 

Where  ‘L’  is  the  length  of  the  free  shear  layer  and  is  assumed  equal  to  convolution 
pitch,  co  is  the  response  frequency  and  V  p  is  the  mean  flow  velocity.  Gerlach  argued 
that  the  characteristic  length  used  for  computation  shall  be  convolution  pitch  instead 
of  convolution  width.  Using  1  =  3mm  and  velocities  corresponding  to  amplitude  peaks 
in  four  modes,  he  computed  the  Strouhal  number  as  0.45.  These  results  obtained  by 
Gerlach  assume  ideal  up-stream  and  uniform  flow  through  the  bellows.  However,  for 
higher  flow  velocities  and  non-uniform  flow  conditions  which  are  observed  in  actual 
service,  it  has  been  noticed  that  the  resonant  peaks  occurred  at  Strouhal  number  of 
0.58  -  about  29%  higher  than  that  found  for  uniform  flow.  Also  it  is  seen  from  the 
above  equation  (2),  that  as  the  mean  flow  velocity  (Vp)  of  the  fluid  through  the 
bellows  increases  the  value  of  Strouhal  number  (S)  decreases  -  and  hence  increasing 
the  frequency  (f).  Gerlach  also  concluded  that  Strouhal  numbers  based  on  the  mean 
flow  velocity  are  not  valid. 

In  1972,  Bass  and  Holster  [6]  extended  the  work  of  Gerlach  to  study  the 
vortex  excitation  of  metallic  bellows  with  internal  cryogenic  flows.  It  was  found  that 
internal  cavitation  or  boiling  due  to  heat  transfer  in  the  formation  of  frost  or 
condensation  on  outside  of  bellows  convolutions  had  the  effect  of  damping  the 


vibrations. 
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Later,  Rockwell  and  Naudascher  [41]  in  their  review  paper  suggested  that  the 
actual  excitation  mechanism  in  the  bellows  is  probably  the  free  shear  layer  instability 
over  the  periodic  cavities  created  by  the  bellows  convolutions.  They  also  noted  that 
the  Strouhal  number  reported  for  bellows  was  less  than  one-half  for  a  rectangular 
cavity  and  speculated  that  the  effect  of  rounded  corners  in  case  of  the  bellows  was  to 
reduce  the  predominant  oscillation  frequency. 

However,  Franke  and  Carr  [12],  has  shown  that  damping  the  upstream  and 
downstream  corners  of  rectangular  cavities  was  very  effective  in  reducing  the  free 
shear  layer  of  oscillations.  Indirectly,  this  would  mean  that  coherent  shear  layer 
structures  should  not  be  made  use  of  over  bellows  convolutions.  This  observation  was 
found  to  be  true  by  Gerlach,  if  the  convolutions  were  rigid. 

The  EJMA  also  provided  a  simplified  method  for  calculating  the  natural 
frequencies  of  bellows.  The  approach  followed  by  EJMA  is  on  the  assumption  that 
bellows  can  be  considered  as  a  continuous  elastic  rod.  Added  mass  is  accounted  for 
the  mass  of  fluid  trapped  between  the  bellows  convolutions  as  was  assumed  by 
Gerlach  [14]  for  the  lowest  modes  of  vibration.  No  provision  is  made  in  this  approach 
for  modes  of  vibration  in  which  simple  translation  of  fluid  between  un-deformed 
convolutions,  and  hence  the  assumptions  made  are  not  quite  reasonable.  The  EJMA 
method  therefore,  may  not  be  very  good  for  estimating  the  natural  frequencies  of 
relatively  short  bellows  or  the  natural  frequencies  of  higher  modes  of  vibration. 
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In  1986,  Becht  [5]  has  brought  out  an  excellent  survey  on  the  developments  of 
numerical  and  theoretical  methods  for  predicting  the  response  of  bellows.  The  state- 
of-the-art  was  assessed  covering  the  earlier  work  of  1940  and  that  the  directions 
requiring  further  development  was  discussed. 

In  1989,  Weaver.D.S  and  Ainsworth.  P  conducted  experiments  to  demonstrate 
the  fatigue  failure  of  a  number  of  bellows.  It  was  observed  that  with  high  internal 
flow  velocities,  significant  flow  induced  vibrations  developed  that  led  to  fatigue 
failure  of  bellows. 

V.F.Jakubauskas  and  D.S. Weaver  [16]  studied  the  natural  vibrations  of  fluid 
filled  bellows  modeling  the  bellows  using  axi-symmetric  shell  finite  elements  and 
discretizing  the  region  using  axi-symmetric  triangular  elements.  The  in- vacuo  bellows 
modes  were  used  as  boundary  conditions  on  the  potential  flow  model  for  the  fluid  and 
the  added  mass  determined  for  each  mode  of  vibration.  The  added  mass  thus 
computed  was  used  to  determine  the  natural  frequencies  of  fluid  filled  bellows. 
Results  were  also  obtained  by  conducting  suitable  experiments  and  the  theoretical 
models  were  verified  against  experimental  results. 

In  1997,  Jakubauskas  and  D.S. Weaver  [17]  studied  the  transverse  vibrations 
of  fluid  filled  bellows  through  a  theoretical  model  based  on  Timoshenko  Beam  theory 
and  included  the  effect  of  added  mass  on  an  internal  fluid.  They  developed  an 
analytical  expression  for  comprising  the  natural  frequencies  of  vibration  of  bellows  in 
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the  form  of  Raleigh  Quotient  in  a  way  suitable  for  hand  calculations.  The  results  for 
the  first  four  modes  of  vibration  are  compared  with  experimental  results  as  well  as  the 
results  from  expressions  given  by  EJMA.  It  was  found  that  the  approach  by  EJMA 
was  in  substantial  error,  as  it  neglected  the  effect  of  rotatory  inertia  and  convolution 
distortion  component  of  fluid  added  mass. 

Considering  only  axial  vibrations,  Jakubauskas  and  Weaver  (1996)  [32] 
modeled  the  bellow  as  a  shell  of  revolution  containing  ideal  fluid  and  then  conducted 
finite  element  analysis.  The  results  showed  that  the  EJMA  predictions  were 
reasonably  good  for  the  lowest  modes  and  more  truly  for  long  bellows.  However, 
significant  errors  were  observed,  as  the  bellows  became  shorter.  The  effect  of 
convolution  shape  distortion  on  fluid  added  mass  was  neglected  in  the  EJMA  model 
and  so  the  prediction  of  natural  frequency  by  this  method  gave  erroneous  results. 

Having  demonstrated  the  validity  of  the  theoretical  model,  the  investigators 
examined  the  degree  to  which  the  theoretical  model  represents  an  improvement  over 
the  existing  EJMA  model.  To  this  effect,  calculations  were  made  for  all  the  above 
cases  using  the  Bernoulli-Euler  theory  with  and  without  the  contribution  of  added 
mass  and  also  using  the  simplified  approach  as  given  in  EJMA  code. 

It  was  seen  that  the  theoretical  calculations  using  the  Bernoulli-Euler  theory 
agreed  well  with  EJMA.  However,  the  transverse  frequencies  of  bellows  obtained 
using  EJMA  method  was  an  overestimate,  because  the  rotary  inertia  is  neglected  and 
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it  was  observed  that  this  trend  increased  with  mode  number.  On  the  other  hand 
transverse  frequencies  of  bellows  in  water  were  seen  to  be  better  matching  with 
theoretical  results  than  predicted  in  air  -  especially  for  the  first  mode. 

Hence,  Jakubauskas  and  Weaver  concluded  the  following- 

♦♦♦  Excellent  agreement  between  the  theoretical  methods  and  experiment. 

♦♦♦  Effect  of  shear  component  of  the  bellows  neglected  in  the  analysis. 

♦♦♦  EJMA  model  overestimates  the  transverse  natural  frequencies  of  the  bellows  in 
both  air  and  water. 

♦♦♦  Effect  of  internal  pressurization  on  bellows  natural  frequencies  is  minimal  on  1st 
transverse  mode  and  decreases  with  increasing  mode  numbers. 

♦♦♦  Fluid  flowing  through  the  bellows  up  to  mean  fluid  velocity  of  lOm/s  has 
negligible  effect  on  bellows  natural  frequency. 

♦♦♦  Effect  of  rotary  inertia  and  convolution  distortion  component  of  fluid  added  mass 
depends  on  geometry  of  the  bellows. 

Jakubauskas  [32]  studied  the  axial  vibrations  of  bellows  and  presented  an 
improved  method  of  computing  the  added  fluid  mass  for  bellows.  The  added  mass, 
thus  computed  is  shown  to  consist  of  three  parts  -  one  due  to  convolution  translation 
in  axial  direction,  the  second  associated  with  convolution  distortion  and  the  third 
associated  with  return  flow  in  central  area  of  cross-section  of  a  bellows.  The 
distortion  component  for  a  half  convolution  was  determined  using  finite  element 
analysis  and  results  were  presented  for  a  typical  range  of  bellows  geometry. 
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Azar.R.C,  Chadrashekar.S,  and  Szaban.P  [3]  developed  a  computer  program 
in  determining  the  axial  vibrations  of  metallic  bellows  to  improve  the  performance  of 
bellows  by  altering  the  spread  of  its  natural  frequencies  thereby  avoiding  resonance. 
They  argued  that  the  earlier  techniques  for  determining  the  bellows  stiffness  use 
either  the  beam  theory  or  the  theory  of  plates  and  shells.  The  first  approach  according 
to  them  was  over  simplification  while  the  second  method  derives  equations  that  are 
difficult  to  use.  Hence,  they  derived  equations  based  on  the  Curved  Beam  Theory. 
The  computer  program  developed  by  them  uses  the  Holtzer  method  in  determining 
the  natural  frequencies  and  mode  shapes  for  the  axial  vibration.  They  have  also 
conducted  vibration  tests  on  several  bellows  by  fixing  one  end  of  the  bellows  to  an 
Electro-dynamic  shaker  and  the  other  end  free.  Results  obtained  from  the  computer 
program  and  the  single  degree  of  freedom  model  were  compared  and  validated  with 
the  experimental  results  and  were  found  to  be  satisfactory. 

Jakubauskas,  V.F  [19]  derived  expressions  for  column  instability.  Internal 
pressure  may  cause  a  bellows  to  become  unstable  in  two  different  ways  in-plane 
instability  and  column  instability.  As  the  column  instability  is  most  frequent  case  in 
bellows,  the  same  has  been  dealt  in  detail.  Expressions  are  derived  for  both  single  and 
double  expansion  bellows  in  lateral  and  rocking  modes  respectively. 

Jakubauskas,  V.F  [20]  has  dealt  a  new  approach  of  calculating  the  fluid  added 
mass  for  regular  shaped  plates  and  beams.  He  strongly  advocated  that  the  added  mass 
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term  is  getting  obsolete  since  any  coupled  problem  of  fluid  elasticity  is  now  solved 
directly  by  considering  two  mathematical  separate  problems  -  of  elasticity  and 
hydrodynamics. 

Jakubauskas,  V.F  [29]  explained  the  theory  behind  the  derivation  of  exact 
estimation  of  the  natural  frequencies  for  double  expansion  joints  based  on  classical 
fixed-fixed  end  conditions.  The  solution  is  based  on  the  modified  Timoshenko 
differential  equation. 

Jakubauskas,  V.F  [34]  briefly  explained  the  work  of  the  first  researchers  in 
the  area  of  flow-induced  vibrations  of  the  bellows.  He  also  dealt  with  the  various 
components  in  the  differential  equation-boundary  conditions  and  the  solution  of 
differential  equation.  He  made  a  comparison  with  other  methods  and  found  that  the 
expressions  for  calculation  of  natural  frequencies  of  bellows  based  on  the  Bernoulli- 
Euler  differential  equation  were  not  precise  for  the  dynamic  calculations  of  bellows 
expansion  joints.  The  added  mass  caused  by  convolution  distortion  was  found  to  have 
a  great  influence  for  higher  modes.  The  inside  pressure  in  bellows  must  also  be 
accounted  for.  He  provided  theoretical  calculations  using  the  derived  formula  and 
confirmed  that  the  predictions  of  Morashita  et.al  about  the  influence  of  rotary  inertia 
on  natural  frequencies  of  bellows  were  correct. 

Jakubauskas,  V.F  [23]  showed  that  convolution  distortion  during  bending 
vibrations  of  bellows,  in  addition  to  the  very  large  fluid  added  mass  contained  in  the 
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cross  section  of  bellows,  can  produce  even  larger  added  mass  effects  -  especially  for 
shorter  bellows  OR  for  higher  vibration  modes.  The  solution  for  the  problem  is  found 
using  the  finite  element  analysis. 

Jakubauskas,  V.F  and  Weaver  D.S  [22]  presents  the  results  of  an  analysis  of 
the  fluid  added  mass  in  bellows  expansion  joints  during  bending  vibrations.  The 
added  mass  consists  of  two  parts  -  one  due  to  the  transverse  rigid  body  motion  and 
the  other  due  to  distortion  of  the  convolutions  during  bending.  The  distortion 
component  that  is  neglected  in  previous  works  has  been  considered  here  and  showed 
to  be  important  for  short  bellows  and  higher  vibration  modes.  The  distortion 
component  is  determined  using  finite  element  analysis  and  also  in  the  graphical  form 
for  a  range  of  bellow  geometry.  The  total  added  mass  is  given  in  a  form  suitable  for 
hand  calculations. 

Jakubauskas,  V.F  [25]  presents  a  method  to  calculate  the  spring  stiffness  of  a 
bellow  using  finite  element  analysis.  The  calculation  results  were  then  compared  with 
the  results  obtained  from  the  analytical  methods. 

1.7.4  Research  Work 

This  section  essentially  deals  with  the  scope  and  summary  of  the  research 
work  carried  out  in  the  area  of  vibrations  and  stability  of  elastically  restrained 


expansion  bellows  in  pipelines. 
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1.7.5  Scope  of  Work 

The  scope  of  the  present  work  is  to  theoretically  investigate  the  effect  of 
elastically  restrained  ends  on  the  axial  and  transverse  vibrations  of  expansion  bellows. 
The  work  carried  out  is  relevant,  as  the  end  conditions  of  bellows  in  general  at 
industrial  sites  are  quite  complex  giving  rise  to  unequal  rotations  at  the  ends. 
Physically,  the  system  may  be  represented  as  a  short  pipe  nipple  welded  at  both  ends 
of  bellows  and  then  the  flanges  welded  to  the  other  end  of  pipe  nipple. 

The  previous  works  contributed  by  researchers  Li  Ting-Xin  et  al  and 
Jakubauskas  V.F  et  al.  [38,28]  show  that  they  had  paid  more  attention  to  study  the 
effect  of  classical  fixed-fixed  type  of  boundary  conditions  only  and  obtained  axial  and 
transverse  natural  frequencies  there  on.  It  is  also  seen  that  no  provision  has  been 
made  to  these  effects  in  EJMA  code,  while  computing  the  axial  and  lateral  vibration 
frequencies  for  single  and  Universal  bellows  and  so  may  not  be  treated  as  accurate 
estimates. 

The  bellows  are  modeled  as  an  equivalent  beam  including  rotatory  inertia  as 
has  been  considered  by  earlier  researchers  and  EJMA.  The  general  method  followed 
in  deriving  exact  frequency  expressions  in  the  present  work  is  based  on  “Separation 
of  Variables”  approach  in  order  to  find  the  axial  and  transverse  natural  frequencies  of 
single  bellows.  This  method  is  nearly  well  adopted  in  analyzing  the  vibration  of 


Universal  (double)  bellows  also. 
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The  influence  and  variation  of  the  rotational  elastic  restraint  parameter  T, 
from  10"2  to  1010  is  studied.  It  is  seen  that  the  effect  of  this  parameter  on  natural 
frequencies  is  quite  significant.  The  effect  of  variation  of  the  internal  pressure  and 
flow  velocity  of  fluid  in  the  bellows  is  also  studied.  The  results  obtained  thereon  by 
solving  the  closed  form  frequency  equations  are  then  confirmed  numerically  by  using 
the  finite  element  method  and  are  found  to  be  in  below  10%  of  error  for  all  the  cases. 

The  particular  case  of  bellows  with  fixed-fixed  boundary  condition  is  a 
reduced  case  wherein  the  elastic  restraints  against  rotation  are  compared  well  with  the 
experimental,  numerical  and  theoretical  results  for  both  axial  and  transverse 
vibrations  respectively.  All  the  results  are  thereon  plotted  for  variation  of  rotational 
stiffness,  T,  internal  pressure,  P  and  flow  velocity  versus  axial  and  transverse 
frequencies.  At  various  stages  of  this  thesis  the  results  obtained  in  the  chapters  have 
been  reported  in  references  [42  to  48]  respectively. 
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CHAPTER  2 

ASSUMPTIONS  AND  DIFFERENTIAL  EQUATION  OF 
TRANSVERSE  VIBRATION  OF  BELLOWS 

2.1  The  Modes  of  Natural  Transverse  Vibrations  of  Bellows 

As  mentioned  earlier  in  Chapter  1,  the  U-shaped  convoluted  type  of  expansion 
bellows  can  vibrate  according  to  longitudinal,  shell  and  beam  vibration  modes. 
However,  the  thesis  is  devoted  to  the  beam  type  of  vibration  mode  only. 

It  is  seen  that  the  previous  researcher  Jakubauskas.V.F,  considered  the  bellows  to  be 
welded  to  pipe  ends  or  flanges  that  are  stiff  in  comparison  to  bellows  itself  and  so  had 
assumed  that  the  expansion  bellows  are  fixed  at  both  ends.  However,  this  assumption 
of  considering  the  bellows  fixed-fixed  at  both  the  ends  result  in  overestimation  of 
natural  frequencies. 

Therefore,  the  aim  of  the  present  work  is  to  compute  both  axial  and  transverse 
natural  frequencies  of  single  and  Universal  type  of  bellows  by  considering  the  ends  as 
elastically  restrained  against  rotation. 

Even  in  the  case  of  a  Universal  bellows  expansion  joint,  that  has  in  between 
two  bellow  elements  with  a  smooth  connecting  pipe,  that  is  generally  many  times 
stiffer  than  the  bellows.  It  is  assumed  that  the  pipe  is  elastic  and  is  not  perfectly  rigid. 
Therefore,  the  system  under  consideration  is  represented  as  a  fixed-fixed  system  with 
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elastic  supports  at  the  ends  and  a  rigid  pipe  section  in  the  middle.  There  are  two  types 
of  transverse  vibration  modes  that  can  be  expected  “Lateral”  and  “Rocking”  and  are 
shown  in  Figures  2.1  &  2.2  respectively. 
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Fig  2.2:  Lateral  Modes  -a)  1st  b)  2nd  (above) 

Rocking  Mode  -a)  1st  b)  2nd  (below) 

In  general,  the  Double  or  Universal  type  of  expansion  joint  vibration  problem 
can  be  modeled  by  considering  two  bending  differential  equations  with  eight 
boundary  conditions.  But,  if  it  is  split  in  to  two  separate  problems  according  to  the 
two  types  of  bending  modes,  the  solution  can  be  simplified  substantially  by 
considering  just  one-half,  either  left  or  right  of  the  whole  system.  These  two  problems 
are  described  in  subsequent  chapters  in  detail. 

2.2  The  Differential  Equation  and  General  Solution  in  case  of  Transverse 
Vibrations  of  bellows 

The  differential  equation  of  vibration  of  bellows  installed  in  a  pipeline  -either 
single  or  double  is  given  by  equation  (2.1).  For  the  circular  cross-section  of  a  thin 
beam,  k’  is  the  cross-section  correction  coefficient  and  is  approximately  given  by 


Paidoussis  et  al.  (1986)  given  in  equation  (2.2) 


oj  Id) 
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4 

W 


4 

X 


+  mtot 


a2w 

dt2 


Jpi+pgi|  g4w  +p!i^=0 

r  Gk'J  dx2.dt2  Gk'  dt4 

(2.1) 


6(l  +  U)(l  +  y2)2 

(7  +  6u)^l  +  Y2J  +  (20  +  12u)y2 

(2.2) 

v  =  Poisson’s  ratio 

y  =  Ratio  of  internal  to  external  radius  of  tube 

A  bellows  as  a  pipe  is  assumed  with  following  dimensions  is  considered  for  analysis 
-Mean  radius  Rm  =  0.03465m,  wall  thickness  t=  0.00028m  and  v=0.25,  then  Rint 


Rint 

Rext 


0.992 


=0. 03446m  and  Rext  =  0.03474m 

From  the  above  data  and  according  to  formula  given  in  equation  (2.2),  k’  =  0.53.  We 
know- 


2(1 +  v) 

The  shear  term  coefficient  in  equation  (2.1)  becomes  - 


pEI 

Gk7 


2(1  +  o) 


I 

k’ 


(2.3) 


(2.4) 
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=  2(1  +  0.25)  -=£-  =  4.72pl 

0.53 

Since  the  rotatory  inertia  coefficient  in  equation  (2.1)  is  equal  to  pi,  the  ratio 
of  shear  to  rotatory  inertia  is  given  by  equation  (2.5) 


( pEI^  1  „ 

vGk'J  pi 


(2.5) 


This  shows  that  for  a  beam  (pipe/tube)  the  influence  of  rotatory  inertia  is  4.72  times 
less  than  the  influence  of  shear.  Therefore,  in  short  beams  i.e.  in  case  of  short  bellows 
the  problems  of  rotatory  inertia  can  be  neglected. 


Now  let  us  consider  a  bellows  as  a  pipe  with  convoluted  surface  as  shown  in 
Fig. 1.2.  Let  the  bellows  have  the  same  dimensions  as  before  with  n  =r2  =0. 00125m, 
convolution  depth,  h  =  0.00571m  and  length  of  bellows,  L=0. 00321m  respectively. 
No  pressure  of  O.OMPa  is  assumed  as  atmospheric  pressure  through  out  the  thesis. 

The  axial  stiffness  of  a  bellows  is  calculated  using  the  expression  given  below- 

k  =  4  R  ra  E  (t  /  h) 3  (2.6) 

Substituting  the  values  and  considering  the  material  of  construction  of  bellow  as 
SS304  Grade  stainless  steel,  E=  2.07  x  10  11  N/m2  at  ambient  conditions  - 
k  =  4  x  0.03465  x  2.07x  10  11  (0.00028/0.00571) 3 
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Bending  stiffness,  El  =  lA  k  q  Rm2 

=  %  x  3.383  x  106  x  0.005  x  0.034652 
=  5.078  Nm2 

G  =  0.828  x  10  11  N/m2  is  calculated  using  the  formula  given  in  equation  (2.3) 
,\pEI  /Gk  =  px  5.078  / 0.828  x  10  11  x  0.53  =  p  x  11.57  x  10  11  kg-m 


The  equivalent  second  moment  of  area  or  moment  of  inertia  of  bellow  cross-section  lb 
is  given  by  - 

_  7lRm3  t.(27in  +  2L) 

lb  —  - 

4n 


Where,  ‘t’  is  the  thickness  of  bellows 


(2.7) 


n  x  0.034653  x  0.00028  (2n  x  0.00125  +  2x  0.00321) 
4x0.00125 

=  1.045  x  10  7  m  4 


.'.  The  ratio  of  shear  rotatory  inertia  is  as  follows- 

'pEf|  J_ 

vGk'J  plb 

px  11.57  x  10  41 

=  -  =  11.072  x  10 -4 

px  1.045  x  10 -7 
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This  ratio  shows  that  for  a  corrugated  pipe  (bellows)  with  a  moderate 
convolution  depth,  the  influence  of  shear  coefficient  is  just  a  small  fraction,  0.00111 
of  the  rotatory  inertia  -without  taking  into  account  the  fluid  rotatory  inertia. 
Therefore,  for  the  investigation  of  transverse  vibrations  of  bellows,  the  influence  of 
shear  can  be  ignored,  and  only  the  rotatory  inertia  is  accounted  for. 

Let  us  now  consider  the  coefficient  of  the  last  term  in  the  differential  equation  (2.1), 

p2I  _  5.078p2 
Gk'~  EG.k' 

(2.8) 

It  is  seen  that  as  the  value  of  EG  in  the  denominator  is  large  compared  to  the 
numerator  and  so  this  coefficient  is  very  small  in  comparison  with  coefficients  of 
other  terms  in  equation  (2.1)  and  so  is  ignored. 


The  influence  of  rotatory  inertia  on  the  frequency  of  natural  vibrations  of 
bellows  is  studied.  For  the  sake  of  comparison  of  natural  frequencies  with  the  results 
obtained  by  Jakubauskas.V.F,  bellows  having  the  same  geometrical  dimensions  are 
considered  here  too. 

The  differential  equation  (2.1)  without  the  shear  term  is  - 


El 


+  ITltot 


8  w 

dt2 


g4w 

dx2.dt2 


=  0 
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Where  J  is  the  mass  moment  of  inertia  of  bellows  cross-section  or  rotatory  inertia  per 
unit  length  and  m  tot  is  the  total  mass  of  the  bellows  that  includes  the  bellow  material 
mass  and  mass  of  the  fluid  flowing  through  the  bellows. 

Let  us  assume  the  first  mode  for  simply  supported  ends  to  be  defined  by  equation 

(2.10)- 

w  =  Asm  re  —  cos cd it 

L 

(2.10) 

Substituting  equation  (2.9)  into  the  differential  equation  (2.8)  we  get  - 


(2.11) 

Leads  to  the  frequency  expression  that  is  as  follows  - 

(2.12) 


K  = 
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Where  K  is  the  coefficient  that  takes  care  of  the  rotatory  inertia  of  bellows.  If  J=0. 
K=1  and  the  frequency  equation  becomes  well  known  frequency  equation  of 
Bernoulli-Euler  equation.  In  order  to  compute  the  value  of  K,  it  is  required  to 
calculate  J  and  total  mass 

m  tot. 

n  Rm3 1  (tu'i  +  7ii'2  +  21)  71  Rm3  (l+ri+r2)  (2r2-t) 

J  =  . . . -pb  +  - pf  (2.13) 

2(ri+r2)  2(n+r2) 

If  ri=r2=0. 00125,  bellows  length  L  =  0.0693m,  pb=  7560kg/m3  and  pf  =  1000kg/m3 
and  substituting  all  the  values  in  the  above  equation  (2.13)  we  get  J  =  0.001 153kg-m. 

The  total  mass  of  the  bellows  is  calculated  according  to  equation  (2.14)  - 

ttRm  (7t  ri  +  n  r2  +2L) 

m  tot  =  - - - pb  +  TrRm2  pf  (2. 14) 

ri  +r2 

Substituting  the  values  above  we  get  mtot  =  5.139  kg/m. 

Now  substituting  the  numerical  values  for  J  and  mtot  in  the  expression  -  we  get  the 
value  of  as  K  =  0.827.  Since  the  value  of  K  is  less  than  unity,  it  is  inferred  that  the 
inertia  of  rotation  of  the  cross-section  lowers  the  natural  frequency  in  comparison 
with  the  Bernoulli-Euler  solution  by  about  17.3%  for  simply  supported  bellows.  A 
similar  work  is  carried  out  for  fixed-fixed  elastically  restrained  ends  against  rotation. 
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2.3  The  Influence  of  Internal  Pressure  and  the  Centrifugal  force  of  flow  on 
Transverse  Vibrations  of  Bellows 

The  differential  equation  is  derived  that  takes  into  account  the  static  inside 
pressure  and  the  centrifugal  force  of  flow  of  fluid  in  addition  to  the  Bernoulli-Euler 
conditions.  The  differential  element  of  such  a  pipe  is  shown  in  Fig  2.1  The  moment 
equation  for  bellows  with  respect  to  point  O  is  given  as  - 


3M 

5x 


.dx  —  Q.dx  +  P7iRm2dw  =  0 


(2.15) 


From  where  the  shear  force  equation  Q  of  bellows  is  obtained  as  follows  - 


Q 


<9M 

<9x 


+  PTTRin2 


5w 

<9x 


(2.16) 


The  translational  differential  equation  with  respect  to  the  y  direction  is- 


f)2  w 

mb  — —  .dx  =  -pf .A  min .  V"  9  .dx  — —  dx 

at2  dx2  5x 


(2.17) 


a2w  ^ 

nti - —  .dx 

at2 


a2w 


aQ 


pf.Amin.V~  o  .dxH - dx  =  O 

ax2  ax 

(2.18) 


A  min  is  the  cross-section  of  the  bellows  and  is  given  by- 
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A  min  =  n  {R  m-  L/2  -n  — t/  2 } 2 


(2.19) 


From  beam  theory,  we  know  - 


M  =  El 


d2w 


d 


x 


(2.20) 


Substitution  of  equation  (2.20)  into  equation  (2.15)  and  subsequently  substituting  the 
equation  (2.15)  into  equation  (2.17)  yields  the  differential  equation  for  transverse 
vibration- 


EI 


8  w 


+  PtiRitl2  +  P  'A  mi  n  V 


x 


2  S2  W 

^X2 


-+-  nn. 


8~  ™ 

St2 


(2.21) 


Therefore,  it  is  seen  from  the  above  equation  that  both  the  pressurization  and 
centrifugal  force  effects  on  the  vibrating  bellows  are  similar  to  that  of  both 
coefficients  of  the  curvature  term  S2w  /dx2 


El 


cf  w 

<9x4 


-I-  TT\ 


3"  w 


-I-  mt. 


c?  w 

3  t2 


(2.22) 

Where  q  =Pti  R2m  the  pressure  coefficient  or  pf  A  min  V2  the  centrifugal  coefficient  or 
both  taken  together. 


Division  of  the  above  equation  (2.22)  by  El  gives- 
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c)4  w  T“i  0~  w  mb  w 

ax4  +  m  +  m  at2 


if 


Then  the  differential  equation  2.23  can  be  written  as  follows  - 

04W  52w  02w 

- +  2c2  - h  a4  —  =  0 

dx4  dx2  dt2 


(2.23) 


(2.24) 


Let  w  =  X  (x).  T  (t) 

Then,  if  T  (t)  is  some  harmonic  function,  the  derivatives  of  (2.24)  needed  in  equation 
(2.23)  are  as  follows- 


And 


<94w 

dX4 

=T  (t) - 

dx4 

dx4 

02W 

dX2 

-  T  O') 

—  I  {l). 

dx2 

dx2 

02W 

=  -  co2X  .T  (t) 

0t2 


Substituting  the  above  derivatives  into  equation  (2.24)  gives 
d4X  d2X 


- h  2c2  . 

d  x4  dx2 


+X4.  X  =  0 


(2.25) 
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Where  A,4  =  a'4  co2  (2.26) 


Let  X  =  C  e  sx 

d4X  d2X 

- =  Cs4  e sx  and  =  Cs2esx  (2.27) 

d  x4  dx2 

Substitution  of  the  above  derivatives  into  equation  (2.24)  we  get  a  quartic  equation 
s4  +  2c2s2  -  A,4  =  0  and  the  roots  of  which  are  as  follows  - 
s  1,2  =  ±  a 

S  3,4  =  +  iP 

Where  a  =  V-c2  +Vc4  +  A,4  (2.28) 

p  =  Vc2  +  Vc4+X4  (2.29) 


Solution  of  the  equation  (2.24)  is  written  as  - 

X  =  Ci  e  a'x  +  C2  e  “a’x  +  C3  e  ^  +  C4  e  4px  (2.30) 

B  +  A 

Ct=  - 

2 

B  -  A 

C2=  - 

2 

D-iC 

C3  =  — . 

2 


D  +  iC 

C4=  - 

2 


.‘.Equation  (2.30)  can  be  written  as  follows- 
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X  =  A  sinh  ax  +  Bcosh  ax  +  CsinPx  +  Dcos  Px  (2.31) 

And  the  derivative  of  which  becomes 
dX 

. =  Aa  cosh  ax  +  Ba  sinh  ax  +  CP  cosPx  -  Dp  sin  Px  (2.32) 

dx 

The  procedure  adopted  to  derive  the  solution  of  differential  equation  is  given  by 
considering  an  example  for  fixed-fixed  case  of  bellows  and  is  assumed  to  be  same  as 
pipe.  The  four  required  boundary  conditions  are  given  by- 

At  X=0 

X  (0)  =0 

dX(0)  _  Q 
dx 

At  X=L 

X(L)  =  0 

dX(L)  Q 
dx 

Substitution  of  the  boundary  conditions  results  in  a  system  of  linear  simultaneous 
equations  with  respect  to  A,  B,  C  &D  respectively. 


B  +  D  =  0  (2.33) 

Aa  +  Cp  =0  (2.34) 

A  sinh  aL  +  B  cosh  aL  +  C  sinpL  +  D  cospL  =  0  (2.35) 

Aa  coshaL  +  Ba  sinh  aL  +  CP  cospL  -  Dp  sinpL  =  0  (2.36) 
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For  a  non-trivial  solution,  the  determinant  formed  by  the  coefficients  of  this  equation 
must  be  equal  to  zero- 


0  1  0 

a  0  p 

sinh  aL  coshaL  sin  PL 

acoshaL  asinhaL  PcosPL 


1 

0 

cosPL 
-Psin  PL 


=  0 


(2.37) 


The  expansion  of  the  determinant  will  result  in  the  frequency  equation  of  the 
differential  equation  - 

(P/a  -a/p)  sinh  aL  .  sin  pL  +  2  cosh  aL  .  cosPL  -2  =  0  (2.38) 

The  influence  of  inside  pressure  is  also  investigated  - 

r|  =  Ptt  Rm2  (2.39) 

According  to  EJMA,  the  critical  pressure  of  bellows  is 

Per  = 

L2 

Where  k  is  the  axial  stiffness  of  the  bellows  and  q  is  the  pitch 

.'.  P  max  =  P  cr  /  6.666  =  n  k  q  /  6.666  L2  (2.40) 

We  know  that  El  =  lA  k  q  Rra2  (2.41) 

Substituting  equation  (2.40)  into  equation  (2.41)  and  subsequently  into  expression  for 
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c  =  0.5477  n  /  L 

Bellows  length,  L  =  2  Rm.  So  in  this  case  L  =0. 0693m 
,\c  =  24.831  L/nr 

Now,  using  the  numerical  value  of  c  and  expressions  a  and  [!  the  frequency  equation 
(2.38)  can  be  solved  numerically. 

The  first  frequency  gave  X  =  65.61205 

From  2.26,  co  =  X2  /  a'2 

Substituting  the  expression  for,  a'  and  the  calculated  value  of  X  into  the  above 
equation  -  the  first  mode  frequency  of  the  bellows  for  fixed-fixed  case  is  calculated  - 
co  =  4305  a/  {El /  mb} 

The  natural  frequency  of  bellows  when  inside  and  outside  pressure  difference  is  equal 
to  zero,  then  the  second  term  in  the  differential  equation  (2.24)  is  cancelled  and  the 
simplified  equation  for  the  first  mode  of  fixed-fixed  beam  becomes  - 

co  =  4.732  /  L2  a/  {El  /  m  b} 

For  Bellows  length,  L  =  0.0693m:  co  =  4658.68  a /  {El  /  mb} 

So  comparison  of  both  the  frequencies  shows  that  taking  into  account  the  maximum 
pressure  allowed  by  EJMA  standard,  the  natural  frequency  becomes  lower  by  about 


7.6% 
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CHAPTER  3 

AXIAL  VIBRATIONS  OF  SINGLE  BELLOWS  EXPANSION 
JOINT  RESTRAINED  AGAINST  ROTATION 

3.1  Investigations  of  Axial  Vibrations  in  Bellows 

The  previous  work  contributed  by  Li  Ting-Xin  et  al,  in  the  area  of  axial 
vibrations  of  bellows  dealt  with  classical  fixed-fixed  end  conditions.  The  present 
work  considers  the  ends  as  elastically  restrained  against  rotation  and  aims  at  finding 
out  the  effect  of  the  rotational  restraint  parameter,  T  on  the  axial  natural  frequencies. 
The  analysis  considers  finite  stiffness  axial  restraints  on  the  bellows,  i.e.,  solving  the 
set  of  equations  with  non-homogeneous  boundary  conditions.  The  bellow  specimen 
considered  here  for  comparison  has  the  same  dimensions  as  taken  by  Li  Ting-Xin  in 
his  analysis.  The  transcendental  frequency  equation  is  derived  and  exact  as  the  first, 
second  and  third  mode  frequencies  computed  are  in  close  agreement  to  the  ones 
obtained  by  Li  Ting-Xin. 

The  paper  contributed  by  Li  Ting-Xin, et  al  [38],  presents  equations  to 
calculate  the  axial  and  lateral  natural  frequencies  of  single  bellows  with  three  types  of 
end  conditions-  one  end  fixed  and  other  end  free,  one  end  fixed  and  other  end 
attached  to  weight  and  both  ends  of  bellows  fixed.  The  theoretical  results  were  then 
compared  with  experiments  performed  on  bellow  specimens  having  different 
geometrical  parameters.  Though  the  error  corresponding  to  the  experimental  value 
was  found  to  be  reasonably  close  to  theoretical  value,  the  end  conditions  represented 
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by  Li  Ting-Xin  do  not  represent  most  of  the  practical  situations.  It  is  seen  on  several 
occasions  that  bellows  ends  are  welded  to  a  small  pipe  spool  that  has  a  lumped  mass 
such  as  a  valve  or  an  instrument.  The  physical  system  considered  for  the  present  work 
in  axial  vibrations  is  shown  in  Fig  3.1.  A  weight  W0  is  attached  as  shown  and  is 
connected  to  a  spring  of  stiffness  Ri.  It  is  assumed  that  a  straight  pipe  is  partially 
fixed  at  the  bottom  and  connected  to  a  spring  of  stiffness  R2.  The  various  geometrical 
parameters  of  the  U  shaped  bellows  used  in  the  analysis  are  shown  in  Fig  1.2. 


wo  hwwH^ 


V 


Fig  3.1:  Physical  System  showing  End  Supports  Elastically  Restrained 


The  geometrical  dimensions  of  U  shaped  bellows  are  given  below  in  Table  3.1.  In 
order  to  compare  the  frequencies,  the  same  geometrical  dimensions  taken  by  Li  Ting- 
Xin  are  considered  here  [38]. 

Table  3.1  Geometrical  Dimensions  of  bellows  (mm) 


D  b 

h 

q 

n 

t 

N 

322.5 

24.5 

22.4 

T 

0.49 

9 

Where  Db  is  outside  diameter  of  bellows,  h  -depth  of  convolution,  q  -pitch,  n- number 
of  plies,  t-thickness  of  bellows,  n-number  of  plies  and  N-number  of  convolutions 


respectively. 
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3.2  Differential  Equation  for  Axial  Vibration  of  Bellows 


The  bellows  are  considered  as  straight  pipe.  The  differential  equation  to  express  the 


axial  vibration  for  the  straight  pipe  is  given  by 


=  2  3  u 

Dt  a 


Where  ‘u’  is  the  axial  displacement  of  the  pipe  (mm)  and 


‘t’  is  the  time  (s),  E  is  the  elastic  modulus  of  bellows  material  (MPa),  g  is  the 
gravitational  acceleration  (9806.65mm  s"2),  v  is  the  weight  per  unit  volume  of  the 
bellows  material  (N/mm3)  respectively. 

3.3  Derivation  of  Natural  Frequency  Equation  for  Axial  Vibrations  of  Elastically 
Restrained  Single  Bellows 

Considering  the  axial  vibrations  of  a  single  bellows  expansion  joint  the  two  boundary 
conditions  for  the  system  are  as  follows- 
We  know, 


At  x  =  0, 


u(x)  -0 
u(t)=0 


At  x  =  0  AE  —  =  kp  u 

dx 
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3u  W0  d2u  , 

At  x  =  L  AE - 1 - -  — —  +  ksu  =  0 

5x  g  at2 

The  exact  solution  of  the  differential  equation  is  written  as  follows  - 
u  (x,  t)  =  {C  sin  p+  D  cos  Px}.e  ltot 
Where  p  =  coi  /  a 


(3.5) 


(3.6) 

(3.7) 


And  coi  is  circular  frequency  in  radian/s,  C  &  D  are  integration  constants.  Now 
applying  the  boundary  conditions  and  substituting  t=0  we  get  - 


u  (x)  =  C  sin  Px  +  D  cos  Px 


(3.8) 


=  P(Ceos  J3x  1) sin  (3x) 

<9x 


(3.9) 

Applying  the  first  boundary  condition  and  substituting  for  du/dx  and  u  (x)  in  equation 
(3.4),  we  get  - 


AE 


<9u 

<9x 


kp.u 


(3.10) 


<9u(x)  _  cos  px  _  | )  sin  px:) 

<3x 


(3.11) 


AE  P(C  cos  Px-D  sin  px)  =  k  P(C  sinPx  +  D  cosPx) 


(3.12) 


At  x  =0,  i.e. 
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AE  011  (x  =  0) 

Ox 


kp.u 


.‘.Substituting  x=0  in  equation  (3.13)  we  get, 


(3.13) 


(p.  AE)  C  =  kp.  D 


( 


.-.c  = 


kp 


A 


p.AE 


D 


(3.14) 


(3.15) 


Now  applying  the  second  boundary  condition  and  substituting  for  u  (x),  5u /8x  and 
52u/5t2  in  equation  (3.5)  we  get- 


*.*  u(x,  t)  =  u(x)  e1CDt 

(3.16) 


AE.p(C  cos  px  -  D  sin  px)  + 


( 


ks 


Wo 


2A 


CD 


(C  sin  (3x  +  D  cos  Px)  =  0 


(3.17) 


At  x=  L  equation  (3.17)  becomes- 
Where  L  is  the  bellows  length  we  get- 


AE.p(C  cos  PL  —  D  sin  PL) 


f 


ks 


Wo 


2A 


CD 


g  ; 


(C  sin  PL  +  D  cos  PL)  =  0 
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AE.p(C  cos  Px  —  D  sin  Px)  — 
ks(C  sin  Px  +  D  cos  Px)  =  0 


Wo  co2 
g 


(C  sin  Px  +  D  cos  Px)  + 


(3.18) 


C 


AE.pcos  pL 


ks 


Wo 


2  A 


CO 


sin  PL 


Di  AEpsin  pL 


ks 


Wo 


2  A 


CO 


cos  pL 


(3.19) 


C 


AE.pcospL  + 


,  Wo©  .  rt 
ks - sin  pL 


V 


sin  PL 


J 


D  AEpsin  pL- 


ks  - 


Wo 


2) 


CO 


cosPL 


(3.20) 


Substituting  (3.15)  in  (3.20)  we  get 


f  1  ) 

kp 

4 

UaeJ 

1L 

AE.pcosPL  + 


ks  -  — °b:>  sin  PL 


y 

1] 

rD 

)_ 

JJ 

D  AEpsin  pL- 


Aks_Wo©A 


v 


cos  PL 


J 


(3.21) 


Cancellation  of  D  on  both  sides,  we  get- 


'  kp  ^ 
VP-AEy 


AE.p  cosPL + 


ks  ■ 


Wo©2 


sin  pL 


■  =  \  AEP  sin  PL 


Aks_Wo©A 


v 


cos  PL 


J 


(3.22) 


Now  multiplying  throughout  by  AEP  we  get  the  final  frequency  equation  as- 
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f  Woco2' 


AEp  kp  +  ks  — 

tan  (}L  = 


f 


kp 


ks  — 


Wo  or 


+ 


(pae)2 


(3.23) 


P 


kp 


ks 


tan  PL  = 


AE  AE 


Wo  or 
g.AE  _ 


P‘ 


kp  (  ks 
AE  AE 


Wo  or 
g.AE 


AEP  kp  +  ks  - 


Wo 


2  A 


co 


tan  PL  = 


g  ) 


(pAE)2  +  kp 


Aks_  w°or" 

V  g  y 


Dividing  throughout  by  AE  we  get- 


(3.24) 


(3.25) 


tan  PL  = 


PL 


^kp.L  ks.L  Wo.Lco2^ 


AE  AE 


g.AE 


(PL)2 


kp.L  f  ks.L  Wo.Loo 


2  A 


AE 


AE 


g.AE 


Let 


(3.26) 


kp.L  ks.L 

AE  b_  AE 

(3.27) 
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Where  Tp  and  Tb  are  the  rotational  restrained  parameters  of  the  pipe  and  bellows 
respectively  and  are  non-dimensionless.  Also,  we  refer  Tp,  Tb  as  Ti  and  T2  and  kp,  ks 
as  Ri  and  R2  in  the  subsequent  chapters  of  the  thesis,  all  meaning  the  same. 


(PL)  =  coi 


a  ) 


LAu 


AE 


G 

g 


Q 

<9u  /  <9x 


.L 


And  we  know, 


a 


& 


Where 
G=LAv  and 


k  = 


du/  dx 


G_ 

Wo 


(3.28) 


Where  a  is  the  ratio  of  weight  of  bellows  to  the  weight  attached  in  pipeline,  G  is  the 
total  weight  of  bellow  material  with  fluid  in  KG  and  k  is  the  axial  spring  rate  of  pipe 
in  (N/m). 


---  CPD 


OOi 


O 

g-kp 


(3.29) 
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3.4  Results  for  Different  End  Conditions 

Case  1  -  The  frequency  equation  for  one  end  of  bellows  fixed  and  elastically 
restrained  and  the  other  end  free  (R2=0,  Tb=0,  Wo=0  &  Ri  =  oc)  is  given  by- 


Pi{(a(Tb  +  Tp)-Pi2)} 

tan  Pi  =  - 

{a(l-TbTp)  +Tppi2} 


Pt  =  couL^/p/E 


TP 


Ri.L 

E.A 


lp  = 


EV 

K3r 


Pi  tan  pi  =  a 


We  know, 


f ,  = 


Wi 

2tt 


(3.30) 


(3.31) 


.-.  f  i  =  49.5  (i-0.5)  Vkn/G 


(3.32) 


Where  Tb  and  Tp  are  dimensionless  rotational  restraint  parameters,  W0  is  the 
weight  attached  to  the  pipe  and  Ri  and  R2  are  rotational  stiffness  of  pipe  and  spring 
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respectively.  The  values  of  [1,  for  the  modes  of  vibration  N=l,  2  &  3  and  increasing 
order  of  Tp  from  0.01  to  109  are  found  out  and  the  frequency  equation  solved. 

A  computer  program  in  FORTRAN  is  developed  to  find  the  solution  of  the 
frequency  equations  and  compute  the  value  of  Pi.  The  values  of  P;  for  the  modes  of 
vibration  N=l,  2  &  3  and  increasing  order  of  Tpis  given  in  Table  3.2. 


Table  3.2:  Case  1  (One  end  fixed  and  other  Free 


C/5 

=0,  Ts=0 

,  Wo=0  &  Tp=  oc 

TP 

0i,  N=1 

Pi,  N=2 

Pi,  N=3 

0.01 

1.502 

3.17331 

6.3465 

0.1 

1.5023 

3.4761 

6.8862 

1.0 

1.5022 

4.4934 

7.7725 

10.0 

1.5044 

4.6910 

7.8412 

102 

1.5644 

4.71026 

7.8527 

103 

1.5701 

4.71217 

7.853 

104 

1.5707 

4.71236 

7.8539 

105 

1.5707 

4.71238 

7.8539 

106 

1.5707 

4.71238 

7.8539 

107 

1.5707 

4.71238 

7.8539 

10s 

1.5707 

4.71238 

7.85398 

109 

1.5707 

4.7123 

7.8539 

Fig  3.2.  Axial  Vibrations  of  Single  Bellows  Expansion  Joint 


The  value  of  pi  obtained  at  Tp=  oc  is  1.5707  for  N=l;  4.7123  for  N=2  and 


7.8539  for  N=3  respectively.  For  Case  1,  Wo=0,  a=0  and  R2  and  Ts=0.  It  is  seen  that 
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as  the  value  of  elastic  restraint  Tp  of  the  pipe  increases  from  0.01  to  109,  the 
frequencies  tend  to  increase  for  all  the  modes  of  vibration.  As  the  restraint  parameter 
T  approaches  oo,  the  frequency  increases  by  about  54%  for  N=1  (P=0.0MPa)  and  by 
69%  for  N=1  (P=166.0MPa)  respectively.  However,  it  is  observed  that  there  is  no 
change  in  frequency  and  it  becomes  constant  from  T=104  onwards.  Table  3.3  gives 
the  values  of  axial  frequencies  in  air  and  water.  P  is  the  internal  pressure  in  the  bellow 
and  O.OMPa  represents  atmospheric  pressure  and  P=166.0MPa  is  the  maximum 
critical  pressure  up  to  which  the  bellows  can  be  subjected. 


Tab 


le  3.3  Axial  frequencies  in  Air  and  Water  for  Case 


SP 

Air 

Water 

1 

Hz 

N=1  N=2  N=3 

Hz 

N=1  N=2  N=3 

40.5  121.5  202.5 

25.5  79.95  124.5 

1 


Case  2  -  The  frequency  equation  for  both  ends  of  the  bellows  elastically  restrained 
(Wo=0  and  R2=oc)  is  given  by- 

pi  {( (1/Tb  +l/Tp)-|3i2)  /a  (Tb.Tp) 

tan  Pi  =  - 

{ [ l/(Tb  .T  p)  -1]+  (Pr/Tb.a)} 

fi  =  49. 5i  *  A  n/G  (3.33) 

Where  ‘i’  is  the  order  number  of  frequency,  i  =  1,  2,  3..  The  values  of  Pi  for  the  mode 
numbers  N=l,  2  &  3  and  different  values  of  Tp  are  given  in  Table  3.5. 
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The  values  of  Pi  for  the  mode  numbers  N=l,  2  &  3  and  different  values  of  T  p 


(0.01  to  109)  are  obtained.  The  value  of  pi  at  Tp=cc  is  3.1415  for  N=l,  6.2831  for  N=2 


and  9.4247  for  N=3  respectively.  For  Case2,  Wo=0,  a  =  oc  and  R2  &  Tb=  oc.  The  axial 


natural  frequencies  for  this  case  are  presented  in  Table  3.4. 


Table  3.4  Axial  frequencies  in  Air  &  Water  for  Case  2 


SP 

Air 

Water 

Hz 

Hz 

N=1 

N=2 

N=3 

N=1 

N=2 

N=3 

1 

81.0 

162.5 

243.0 

51.1 

101.2 

149.4 

Tab 


le  3.5  Case  2  (Both  end 

Is  fixed  A 

f 

■nm 

■HESS 

IE125B 

0.01 

1.57713 

4.71451 

7.85525 

0.1 

1.63199 

4.733518 

7.86669 

1.0 

2.02875 

4.91318 

7.97866 

10.0 

2.86277 

5.76055 

8.70831 

102 

3.11049 

6.22105 

9.33172 

103 

3.13845 

6.27690 

9.41536 

104 

3.14127 

6.28257 

9.42383 

105 

3.14156 

6.28312 

9.42468 

106 

3.14158 

6.28317 

9.42476 

107 

3.14159 

6.28318 

9.42477 

108 

3.14159 

6.28318 

9.42477 

109 

3.14159 

6.28318 

9.42477 

The  value  of  Pi  for  elastically  fixed-fixed  condition  and  for  T  p=>oo  is  3.1415  for  N=l, 


6.2831  for  N=2  and  9.4247  for  N=3  respectively.  The  weight  Wo=0,  a=  oc  and  k  s  & 
Ts=cc. 


Case  3  -  One  end  fixed  and  other  end  attached  to  weight  represents  a  situation  where 
one  end  of  the  bellows  is  welded  to  a  pipe  and  the  other  end  is  attached  to  a  lumped 
weight  Wo  (Wo^O,  and  R2=0)  either  in  form  of  a  valve  or  an  instrument.  Therefore, 
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the  values  of  pi  for  this  type  of  end  condition  are  different  from  the  other  two  cases 

and  are  obtained  by  varying  Tp  and  keeping  1/a  (0.1,  1,  100  &  1000)  a  constant 

value.  It  depends  on  the  value  of  a.  The  frequency  equation  for  this  case  is  given  by- 

pi  {Tb  +Tp-(l/a)pi2} 

tan  pi  =  - 

{ 1-Tb.Tp  +  (1/a)  Tp.pr} 

The  values  of  Pi  are  found  for  different  values  of  1/a  and  are  as  follows. 

For  1/a  =0.1,  pi  =3. 14159  (N=l),  6.28318  (N=2)  &  9.4247  (N=3) 

For  1/a  =1.0,  pi  =3. 14159  (N=l),  6.28318  (N=2)  &  9.4247  (N=3) 

For  1/a  =100,  p;=3. 14159  (N=l),  3.14159  (N=2)  &  6.28318  (N=3) 

For  1/a  =1000,  pi=1.5671  (N=l),  3.14159  (N=2)  &  6.28318  (N=3) 

By  substituting  the  values  of  Pi  in  the  frequency  equation,  the  axial  natural 
frequencies  are  obtained. 


Table  3.6  Axial  frequencies  in  Air  &  Water  for  Case  3 


SP 

a 

Air 

Water 

Hz 

N=1  N=2  N=3 

Hz 

N=1  N=2  N=3 

1 

0.1&1.0 

81.0  162.5  243.0 

51.1  101.2  149.4 

1 

100&1000 

40  81  161 

25  50  99.5 

Table  3.4  presents  axial  frequencies  for  different  values  of  1/a  =  0.1,  1.0,  100 
and  1000  respectively.  It  is  seen  that  the  values  of  Pi  in  both  the  cases  tend  to  be 
constant  from  Tp  =  105  onwards.  However,  for  lower  values  of  1/a  of  0.01  and  0.001, 
which  means  the  bellows  weight  is  greater  than  the  weight  attached  (G>Wo),  the 
values  of  Pi  almost  become  constant  for  N=2  &  3  and  vary  for  N=1  only.  Axial 


81 


natural  frequencies  in  air  and  water  are  obtained  for  specimen  number  1  by  exact  are 
given  in  Table  3.6.  It  is  seen  that  the  frequencies  obtained  by  exact  method  has  a 
percentage  of  error  of  7.2%  and  7.4%  in  air  and  water  respectively  in  comparison  to 
the  experimental  and  analytical  results  obtained  by  Li  Ting’s  on  specimen  (SP-1) 

[15]. 


Table  3.7  Comparison  of  Axial  Frequencies  in  Hz 


A 

r 

SP 

EXP 

Li  Ting 

Exact 

%  error 

[38] 

[38] 

[42] 

1 

37.5 

40.4 

40.5 

7.4 

Water 


SP 

EXP 

Li  Ting 

Exact 

%  error 

[38] 

[38] 

[42] 

1 

27.5 

27.5 

25.5 

7.2 

Fig  3.3.  Axial  Vibrations  of  Single  Bellows  Expansion  Joint 
(Both  ends  elastically  fixed) 

Bta(i)  for  N=l,  2  &3,  Wo=0,  Tb  &  alpha=infinity) 


Tp 


Tables  3.8  &  3.11  present  Pi  values  for  1/a  of  0.1  and  1.0  respectively.  It  is  seen  from 
Figures  3.5  &  3.6  that  the  values  of  Pi  in  both  the  cases  tend  to  be  constant  from  Tp  = 
105  onwards.  However,  for  lower  values  of  1/a  of  0.01  and  0.001,  which  means  the 
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bellows  weight  is  greater  than  the  weight  attached  (G>Wo),  the  values  of  pi  are 
almost  constant  for  N=1  &  2  and  vary  for  N=1  only.  The  same  is  shown  in  Figures 
3.7  &  3.8  respectively. 

Table  3.8  Case  3 

(One  end  fixed  and  other  attached  to  weight 
_ l/q=0.1  and  T i,—  oc) _ 


TP 

Pi.N=l 

Pi.  N=2 

Pi.  N=3 

0.01 

1.435 

4.3076 

7.2289 

0.1 

1.4899 

4.3239 

7.2366 

1.0 

1.8964 

4.4899 

7.3172 

10.0 

2.8417 

5.5990 

8.228 

102 

3.1101 

6.2185 

9.323 

103 

3.13845 

6.2768 

9.4152 

104 

3.14127 

6.2825 

9.4238 

105 

3.14156 

6.28312 

9.4246 

106 

3.14158 

6.28317 

9.4247 

107 

3.14159 

6.28318 

9.4247 

10s 

3.14159 

6.28318 

9.4247 

Fig  3.4.  Axial  Vibrations  of  Single  Bellows  Expansion  Joint 
(One  end  of  bellows  fixed  and  other  end  attached  to  weight) 

(Btai  for  N=1, 2&3)  and-alpha=10> 

1.E+01 


9.E+00 

8.E+00 

7.E+00 


6.E+00 

5.E+00  2 

CD 

4.E+00 

3.E+00 

2.E+00 

- Btai  (N=1) 

1  .E+00 - Btai<N=2> 

-  -  -  BtaiJN=3) 

0.E+00 

1  .E+00  1.E+01  1.E+02  1.E+03  1.E+04  1.E+05  1.E+06  1.E+07  1.E+08 

Tb 
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Fig  3.5.  Axial  Vibrations  of  Single  Bellows  Expansion  Joint 
(One  end  fixed  and  other  end  attached  to  weight) 

(Btai  for  N=l,2&3  and  alpha  =1.0) 


Table  3.9  Case  3 

(One  end  fixed  and  other  attached  to  weight 
_ l/q=1.0  and  Tb=oc) _ 


TP 

3i,  N=1 

3i,  N=2 

Bi,  N=3 

0.01 

0.8645 

3.42583 

6.43733 

0.1 

0.9014 

3.42775 

6.43765 

1.0 

1.2077 

3.44823 

6.44095 

10.0 

2.5293 

3.82916 

6.48289 

102 

3.1072 

6.18340 

8.98748 

103 

3.1384 

6.27665 

9.42382 

104 

3.14127 

6.28312 

9.42468 

105 

3.14156 

6.28317 

9.42468 

106 

3.14158 

6.28317 

9.42476 

107 

3.14159 

6.28318 

9.42477 

10s 

3.14159 

6.28318 

9.42477 
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Table  3.10  Case  3 

(One  end  fixed  and  other  attached  to  weight 


l/a=10( 

)  and  Tb=oc) 

TP 

pi. 

N=1 

Pi, 

N=2 

Pi, 

N=3 

0.01 

0.1003 

3.1447 

6.28477 

0.1 

0.1047 

3.1447 

6.28477 

1.0 

0.1411 

3.1447 

6.28477 

10.0 

0.3311 

3.1448 

6.28478 

102 

1.0031 

3.1455 

6.28481 

103 

3.0809 

3.2236 

6.28531 

104 

3.1412 

3.1412 

6.28311 

105 

3.1415 

3.1415 

6.28311 

106 

3.1415 

3.1415 

6.28317 

107 

3.1415 

3.1415 

6.28318 

108 

3.1415 

3.1415 

6.28318 

Table  3.11  Case  3 

(One  end  fixed  and  other  attached  to  weight 
l/q=1000  and  Ti,=qc  ) 


TP 

pi.  N=1 

Pi,  N=2 

Pi,  N=3 

0.01 

0.0317 

3.1419 

6.28334 

0.1 

0.0331 

3.14190 

6.28334 

1.0 

0.0447 

3.14190 

6.28334 

10.0 

0.1048 

3.14190 

6.28334 

102 

0.3177 

3.14191 

6.28334 

103 

1.0003 

3.14191 

6.28339 

104 

1.1123 

3.14151 

6.28308 

105 

1.3112 

3.14155 

6.28317 

106 

1.4234 

3.14158 

6.28318 

108 

1.5671 

3.14159 

6.28318 
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Fig  3.6.  Axial  Vibration  of  Single  Bellows  Expansion  Joint 
(One  end  fixed  and  other  attached  to  a  weight) 

(Btai  for  N=l,  2  &  3  and  alpha=0.01) 


Tb 


Fig  3.7.  Axial  Vibrations  of  Single  Bellows  Expansion  Joint 
(One  end  fixed  and  other  end  attached  to  weight) 

(Btai  for  N=1 , 2  &  3  and  alpha=0.001 ) 


Tb 
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CHAPTER  4 

FINITE  ELEMENT  ANALYSIS  OF  AXIAL  VIBRATIONS  OF 
SINGLE  BELLOWS  EXPANSION  JOINT  RESTRAINED 
AGAINST  ROTATION 

4.1  Theoretical  Background 

The  work  presented  here  are  the  results  of  investigation  of  axial  vibrations  of 
single  bellows  expansion  joint  restrained  against  rotation  on  either  end  using  the  finite 
element  method.  The  aim  is  to  model  flexible  U  shaped  bellows  using  bar/rod 
elements  and  compare  the  results  obtained  by  finite  element  analysis  with  the 
theoretical  frequencies  for  the  first  fundamental  mode  of  vibration. 

Rao.C.K  &  Radhakrishna  M  [42]  analyzed  U-shaped  bellows  subjected  to 
axial  force  and  obtained  an  analytical  solution  for  the  axial  vibrations  of  single 
bellows  that  are  elastically  restrained  against  rotation  on  either  end  using  the  Euler  - 
Bernoulli  beam  theory. 

As  mentioned  before  the  expressions  presented  either  in  EJMA  code  or  those 
derived  by  investigators,  [38]  cover  classical  boundary  conditions  only.  It  is  observed 
that  the  flange-bellow-flange  junction  -  as  fixed  end  conditions  are  supposed  to  be 
infinitely  stiff  compared  to  the  bellows  stiffness  and  so  is  not  a  practical  situation 
compared  to  the  bellows  welded  on  either  end  in  many  of  the  piping  systems.  Hence, 
this  assumption  of  fixed-fixed  will  lead  to  an  over  estimation  of  natural  frequencies. 


87 


4.2  Characteristics  of  Bellows 

The  various  geometrical  parameters  of  U-shaped  bellows  are  given  in  Fig  1.2  n 
represents  the  meridional  radius  of  the  convolution  root,  r2  is  the  meridional  radius  of 
the  convolution  crown  and  h  is  the  convolution  height.  R  m  is  the  mean  radius  of  the 
bellows  -  the  distance  from  the  bellows  centerline  to  mid  convolution  height,  and  t  is 
the  bellow  material  thickness. 

It  is  assumed  that  t<SCn,  r2  and  h  <SC  Rm  as  shown  in  Fig  1.2.  The  total  length  of  the 

bellows  is  L  =  2(ri+r2)  N.  Based  on  the  assumptions  given  above  the  bellows  are 
considered  to  be  as  an  equivalent  rod  or  bar  of  radius  R  m  and  wall  thickness,  t. 

The  other  dimensions  of  the  bellows  are  Db=322.5mm,  depth  of  convolution,  h 
=24. 5mm,  pitch,  q=22.4mm,  number  of  plies,  n=l,  overall  length  of  two  ends  of 
cylindrical  tangent,  L=30mm,  wall  thickness,  t  =  0.49mm  and  number  of 
convolutions,  N  =  9  [2]. 

4.3  Theory  of  Free  Vibrations  of  Bellows 

The  matrix  equation  for  the  free  vibration  of  bellows  can  be  written  as  - 

[M]  { q }  +  [K]  { q }  =0  (4.1) 

Where 

{q}  -Generalized  coordinates 
[M]  -  Mass  matrix 


[K]  -  Elastic  stiffness  matrix 
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4.4  Formulation  of  Elastic  Stiffness  Matrix 

The  strain  energy  U  of  a  bellow  element  of  length  ‘L’  is  given  by- 

1  i 

U=  (Vi  AE/L2)j(0u/5ri)  2.dri  +  Yi  k.  L  J(u)  2.dri=0 
o  o 


(4.2) 


Where- 

u  is  axial  deflection,  k  -  axial  stiffness,  A  -area  of  cross-section  and  E  is  the  elastic 
modulus  of  material  of  bellows  respectively. 

Now,  assuming  a  polynomial  expression  for  ‘u’  to  be  of  the  form  u  =  ai  +  bi  x-  the 
generalized  stiffness  matrix  for  axial  vibrations  is  given  by- 


AE 


1 

1 


1 

1 


1 1 1 
1 1 


(4.3) 

4.5  Formulation  of  Mass  Matrix 

The  mass  matrix  is  found  out  by  using  the  expression  for  kinetic  energy  ‘V’- 


V=  -*/2  (pAl+Wo)  J(5u/0t)  2.dr|  (4.4) 

If  G=pAL,  then  the  above  expression  becomes- 

V=  !/2(G+W0)  l(ou/ot)  2.dr|  (4.5) 

Considering  u  (rj,  t)  =  u  (r|).  e  ipt  we  get,  (4.6) 

1 

V=-  ^(G+Wo)  roV.dri  (4.7) 

0 


The  generalized  mass  matrix  [M]  for  axial  vibrations  is  given  by- 


u2.dr,  =[M] 


pAL 

6 


(4-8) 
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The  matrix  equation  for  free  vibration  of  single  bellows  is  given  by- 

[[K]  -XA  [M]]  {^}  =  0  (4.9) 


Where  £,  represent  the  amplitude  of  displacement  called  the  mode  shape  or  eigen 
vector  and  X  denotes  the  natural  frequency  of  vibration.  It  will  have  a  non-zero 
solution  for  c,  provided  that  the  determinant  of  the  coefficient  matrix 
I  [K]  -co2  [M]  |  =0 


The  bellows  are  modeled  using  the  rod/bar  elements  -we  get  the  resultant  matrix  - 


\  4  \  4 

l  +  Ri-4^-  -1  +  4*-  0 

Ui 

3  6 

-iX  2-^  -i+^ 

U2 

< 

6  3  6 

\  4  \  4 

0  -1  +  4*-  l  +  R2-4^ 

U3 

6  3 

-  - 

(4.10) 

Where  Ri  and  R2  are  the  stiffness  of  supports  at  both  ends  of  the  bellows  and 
X4=(p.L2C0n2)/E 

Now  the  matrix  equation  (4.10)  can  be  written  as- 


"(1  +  Ri) 

-1 

0 

UI 

.  4 

'  2 

-1 

0  " 

Ul 

-1 

2 

-1 

U2 

6 

-1 

4 

1 

U2 

0 

0 

(-I  +  R2) 

U3 

0 

0 

-1 

U3 

(4.11) 


Since  ui  ^0  and  U2  ^  0,U3  ^0.  the  determinant  of  the  matrix  is  equal  to  zero- 
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{(1  +Ri-X4/3)[(2-1  /3X4)(1  +R2-A4/3)-(-1  +r76)2]}-(-1  +X4/6)[(-1  +X4/6)(1  +R2- 
^4/3)]}=0 

(4.12) 


Expansion  of  the  determinant  equation  (4.12)  leads  to  the  final  frequency  equation- 


Rl.R2 

f  1 


1  2 
- 1 - 1 - 

vRi  R2  R1.R2 


36  a. 


4 


r 


1 


H - 1 - hi 

\Ri  R2  R1.R2 


10 


+  —  +  2 


Ri  R2 


=  0 


(4.13) 


Where  Ri  and  R2  are  the  rotational  restraint  parameters  and  have  been  defined  earlier 
It  considers  two  types  of  elastically  restrained  end  conditions  that  are  as  follows  - 


♦♦♦  Case  1  -One  end  is  fixed  and  the  other  end  is  free 
♦♦♦  Case  2  -  Both  the  ends  are  fixed 


The  frequencies  are  obtained  for  the  two  cases  by  considering  geometrical 
dimensions  of  the  bellows  mentioned  above  and  results  are  then  compared  with  those 
obtained  from  closed  form  solutions. 

Applying  the  boundary  conditions  for  both  the  cases  we  get  A  . 


Case  1-  represents  an  end  condition  of  bellows  where  one  end  is  fixed  and  other  free 
and  both  ends  elastically  restrained  i.e.  where  Ri  =oc  or  1/Ri=0,  W=0  and  R2=0. 
Substituting  the  values  for  Ri  and  R2  in  equation  (4.13)  we  get  Ai=1.5507,  A2=4.6123, 
^3=7.6539  respectively.  Also  the  values  for  A  are  found  out  for  three  modes  of 


vibration  by  increasing  the  value  of  Ri  from  0.01  to  oc. 
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Table  4.1  Case  1  (One  end  fixed  and  other  free) 


k  p 

Xi.N=l 

7.2  N=2 

7,„N=3 

0.01 

1.4902 

3.17331 

6.3465 

0.1 

1.4923 

3.3761 

6.4862 

1.0 

1.5022 

4.3934 

7.5725 

10.0 

1.5044 

4.5910 

7.6412 

102 

1.5144 

4.61026 

7.6527 

103 

1.5146 

4.61217 

7.6539 

104 

1.5507 

4.61236 

7.6539 

105 

1.5507 

4.61238 

7.6539 

106 

1.5507 

4.61238 

7.6539 

107 

1.5507 

4.61238 

7.6539 

10s 

1.5507 

4.61238 

7.6539 

109 

1.5507 

4.61238 

7.6539 

Case  2  -represents  a  case  where  both  the  ends  of  the  bellows  are  fixed  and  elastically 
restrained  i.e.  where  Ri  =oc  and  R2=oc  or  1/Ri=1/R2=0  &  W=0  we  get  Xi=3.1316, 
^2=6.283 18  and  7,3=9.42477  respectively.  By  varying  Ri  from  0.01  to  cc  or  109  and 
keeping  R2  constant  the  values  of  X  are  found  out  for  N=l,2&3  respectively. 
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Table  4.2  Case  2 

(Both  ends  fixed) 

TP 

Xi,  N=1 

7.2  N=2 

7.3  N=3 

0.01 

1.57713 

4.71451 

7.85525 

0.1 

1.63199 

4.72518 

7.86669 

1.0 

2.02875 

4.91318 

7.97866 

10.0 

2.86277 

5.76055 

8.70831 

102 

3.11049 

6.22105 

9.33172 

103 

3.12845 

6.27690 

9.41536 

104 

3.13127 

6.28257 

9.42383 

105 

3.13156 

6.28312 

9.42468 

106 

3.13158 

6.28317 

9.42476 

107 

3.13159 

6.28318 

9.42477 

108 

3.13159 

6.28318 

9.42477 

109 

3.13159 

6.28318 

9.42477 

Fig  4.2.  Axial  Vibrations  of  Single  Bellows  Expansion  Joint 
(Both  ends  fixed) 

LMDA  for  N=l,  2  &3  &  Rj  and  R2=infinity) 


< 

Q 


h-3 


- LMDA( 

N=1 ) 

-  -  LMDA 

(N=2) 

-  -  -  LMDA 
(N=3) 


Ri 


4.6  Results  and  Discussion 

It  is  seen  that  the  values  of  X  obtained  using  finite  element  method  are  found 
to  be  in  close  agreement  with  the  values  of  [1,  obtained  from  closed  form  solution  [3]. 
The  value  of  pi  obtained  by  closed  form  solution  is  1.5707  (N=l)  for  Case  1  and 
3.1415  for  Case  2  respectively.  The  number  of  elements  considered  in  the  present 
work  is  two.  By  increasing  the  number  of  elements  to  3  it  is  seen  that  the  values  of  X 
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and  pi  would  converge  and  the  error  percentage  in  frequency  would  decrease 
correspondingly.  Figures  4.1  &  4.2  represent  the  same.  Tables  4.3  &  4.4  presents  a 
comparison  between  frequencies  obtained  by  experiments,  exact  and  finite  element 
methods. 


Table  4.3  Comparison  of  Frequency  Solutions  in  Air 


(One  end  fixed 

l  and  other  free  (R2=0  and  R  i=oc) 

Mode 

N 

EXP 

[38] 

Hz 

Exact  [42] 
©,  Hz 

FEA  [43],  go 
Hz 

Error,  % 

1 

37.5 

40.5 

39.99 

6.4 

Table  4.4  Comparison  of  Frequency  Solutions  in  Air 


(Both  ends  fixed  and  elasticall 

[y  restrained  IC=a:  and  Ri=oc 

Mode 

N 

EXP[38] 

Hz 

Exact  [42] 
go,  Hz 

FEA  [43],  go 
Hz 

Error,  % 

1 

83.0 

81.0 

80.77 

2.6 

The  theoretical  and  experimental  comparisons  confirm  that  the  finite  element 
method  developed  and  used  in  the  study  is  well  adapted  to  the  dynamic  response  of 
bellows  and  gives  fairly  exact  results. 
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CHAPTER  5 

TRANSVERSE  VIBRATIONS  OF  SINGLE  BELLOWS 
EXPANSION  JOINT  RESTRAINED  AGAINST 

ROTATION 

This  chapter  deals  with  the  theoretical  investigation  of  transverse  vibrations  of 
single  bellows  subjected  to  internal  pressure,  rotatory  inertia  and  considering  the  ends 
as  elastically  restrained  against  rotation.  The  mathematical  model  represents  the 
expansion  joint  with  the  ends  restrained  against  rotation  as  shown  in  Fig  5.1.  A  non- 
dimensional  rotational  restraint  parameter  is  defined  in  order  to  derive  a  closed  form 
general-purpose  transcendental  frequency  equation. 


Fig  5.1  Mathematical  Model  of  Single  Bellows  Expansion  Joint 

(Lateral  Mode) 


5.1  The  Differential  Equation 

The  general  form  of  the  differential  equation  of  vibration  of  bellows  for  single  or 
Universal  double  bellows  expansion  joint  is  given  by  - 
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El 


d4w 

5x4 


+  P7tRm2 


d2W 

dx2 


j  a4w 

dx2  dV 


+  ncitot 


a2w 

dt2 


=  o 


(5.1) 

Where  El  is  bending  stiffness,  P-internal  pressure,  J-  total  mass  moment  of  inertia  per 
unit  length,  w  -  deflection,  x-axial  coordinate,  Rm  -  mean  radius  of  bellows,  m  tot  is 
the  total  mass  of  bellows  per  unit  length  includes  mass  of  bellows  and  fluid  mass  and 
t-time  [28].  Using  the  technique  of  Separation  of  variables,  the  transverse  deflection 
of  the  bellows  axis  ‘w’  can  be  expressed  as  - 


w  (x,  t)  =  X  (x)  .  T  (t)  (5.2) 

Where  X  (t)  is  a  parameter  of  x  only  and  T  (t)  is  the  time  harmonic  function  (T  is 
associated  with  time  here)  such  as  - 

T (t)  =  A  e  icot  (5.3) 

Differentiating  the  above  equation  (5.2)  and  substituting  in  the  differential  equation 
(5.1)  we  get  - 

<34X  (Pti  Rm2  + J.  w2)  <32X  2  mtot 

- 2 — I - - - w  - A  = 

a  x4  ei  a  x2  ei 


|f  c=  /( P7tRin2  +  Jc02) 

V  2EI 


d4X 

d>a 


+  2c2 


d2X 

dx2 


(5.4) 

(5.5) 


A.4X  =  0 


(5.6) 


96 


The  general  solution  of  the  equation  is  given  by  - 

X  (x)  =  A  sinh  ax  +  B  cosh  ax  +  C  sinPx  +  D  cospx  (5.7) 

The  first  three  derivatives  of  equation  (5.7)  are  given  as  follows- 


dX(x) 

dx 


=  A  or  cosh  otx  +  Bex  sinh  cxx  +  Cp  cos  Px  +  Dp  sin  Px 


(5.8) 


dX(x)2 

dx2 


2  2  2  ^  2 

Aq^  sinh  ax  +  B^coshax-Cn  sinPx -Dp  cospx 


(5.9) 


Where  the  roots  of  the  equation  are  a  &  P  and  their  values  are  given  by  - 


(5.10) 

P  =  A / G2  ■+■  V G4  ■+■  AJ7 


(5.11) 


A  = 


mtot.  oo“ 


El 


(5.12) 


Where  A,  B,  C,  D  are  arbitrary  constants 
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5.2  Derivation  of  the  Frequency  Equation  for  Transverse  Vibrations  of 
Elastically  Restrained  Single  Bellows  Expansion  Joint 


At  end  ‘A’  the  bellows  are  connected  to  a  pipe  nipple  and  is  considered  to 
have  a  rotational  stiffness  of  Ri  and  R2  at  either  end  as  shown  in  figure  5.1. 


At  x  =  0  the  boundary  conditions  are  given  by- 


X  (0)  =  0  and 


EI  5~x(Q)  =  _T  co) 


8  —  2 


X. 


At  end  B,  the  boundary  conditions  are  given  by  - 
X  =  L 

X  (L)  =  0  and 

x(L)  ....  ax(L> 

=  ' 

Ti  =  Ri  L/EI 

t2  =  r2  l/ei 


Or  can  also  be  written  as- 


X(0)  =0 


d2x(0) 

dx2 


R 


dx(0) 

dx 


X  (L)  =0 


Ej  dx(L) 
dx 1 


R 


dx(L ) 
dx 


Applying  the  boundary  conditions  we  get- 
B  +  D  =0 

EI  [B  (aL)2  -  D(pL)2]  =Ti  [A(aL)  +  C(pL)] 


(5.13) 


(5.14) 

(5.15) 

(5.16) 

(5.17) 


(5.18) 

(5.19) 
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{B  (aL)2-D(pL)2}  =^{A(aL)  +  C(pL)}  (5.20) 

El 

A  SinhaL  +  B  coshaL  +  C  sinPL  +  D  cosPL  =  0  (5.21) 

A{((aL)2  sinhaL  +  T2(aL)  coshaL)}  +  B{(aL)2  coshaL  +  T2  (aL)  sinhaL)}  + 
C{(T2pL 

cosPL  -  (PL2)  sinPL)}  -  D{((pL)2  coshpL  +  T2  pL  sinPL)}  =  0  (5.22) 


Re-writing  the  above  equation  by  re-arranging  the  terms  we  get- 
{ A  (aL)2  sin  haL  +  B(aL)2coshaL  +  C(PL)2sinPL  +  D(PL)2coshpL} 
{AaL  coshaL  +  BaL  sinhaL  +  CpL  coshpL  -SinpL} 


(5.23) 


A  {(aL)  (aLsinhaL  +  T2  coshaL}  +  B  (aL)  {(aL  coshaL  +  T2  sinhaL] 

+  C  (pL)  {T2coshpL  -  (pi)  sinpL}  -  D  (pL)  {(T2sinpL  +  (pL)  cospL}  =  0  (5.24) 
Writing  the  above  set  of  equations  in  matrix  form,  we  get- 


0 

1 

1 

1 

'A' 

Ti(aL) 

-aL2 

Ti(pL) 

PL2 

B 

sinh  aL 

cosh  aL 

sin  PL 

cosPL 

C 

Ci 

C2 

C3 

C4 

D 

(5.25) 

Where,  ci  =  (aL)  [(aL)  sinhaL  +  T2CoshaL] 

C2  =  (aL)  [(aL  coshaL  +  T2sinhaL] 

c3  =  (pL)  [T2  cospL  -  (pL)  sinpL] 

c4  =  (pL)  [T2  sinpL  +  (pL)  cospL] 


(5.26) 
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1 

-S 

+  (P1)2  _ 

Tj  (al)  -(al)2 

T,(P1) 

Sinh  (al)  Sin  (pi) 

Cos(pl) 

+ 

Sinh  (al)  Cosh(al) 

Sin  (pi) 

C1  C3 

-c4 

Cl  c2 

c3 

(5.27) 


Expanding  the  above  equation  -we  get 

Ti(aL)  [-C4  sin(PL)  -  C3  cos(PL)]  -  Ti(pL)[-C4  sinh  (aL)  -  ci  cos(PL)]  +  (PL)2  [C3 
sinh(aL)  -  ci  sin(PL)]  +  Ti(aL)  [C3  cosh(aL)  -  C2  sin(PL)]  +  (aL)2[c3  sinh  (aL)  -  ci 
sin(PL)]  +  Ti  (pL)  [C2  sinh(aL)  -  ci  cosh(aL)]  =0  (5.28) 


Substituting  ci,  C2,  C3  &  C4  in  above  equation  and  simplifying  we  get- 

-C4  sinPL  -  C3  cosPL  =  -  (PL)  sinPL  [T2  sinpL  +  (PL)  cosPL]  -  (pL)  cosPL  [T2  cosPL 

-  (pL)sinpL]  (5.29) 

=  -  (PL)  T2  sin2pL  -  (PL)2  cospL  sinpL  -  (PL)  T2  cos2pL  +  (pL)2  sinpL  cospL 

.'.  [C4  sinpL  -  C3  cosPL]  =  -  T2(PL)  (shrpL  +  cos2pL)  =  -T2(PL)  (5.30) 

Now, 

-C4  sinh(PL)  -  ci  cospL)  =  -(PL)  sinhaL  [T2  sinpL  +  PLcosPL]  -  (aL)  cospL  [(aL) 
sinhaL  +  T2  coshaL]  (5.31) 

=  -T2  (pL)  sinhaL  sinpL  -  (PL)2  sinhaL  cospL  -  (aL)2  sinhaLcosPL- 
T2(aL)coshaLcosPL  (5.32) 

=  -T2  [(aL  coshaL  cospL  +  (pL)  sinhaL  sinpL]  -  [(aL)2  +  (PL)2]  sinhaL  cospL 
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C4  sinhaL  -  cicospL  =  -  T2  [(aL)  coshaL  cospL  +  (PL)  sinhaL  sinpL]  -  [(aL)2  + 
(PL)2]  sinhaL  cosPL  ]  (5.33) 

Now, 

[C3  sinh(aL)  -ci  sin(pL)]  =  (PL)  sinh(aL)  [T2  cosPL  -(PL)  sinpL  -  (aL)  sinpL  [aL 
sinhaL  +  T2  coshaL]  =  T2  [(aL)  coshaL  sinpL  +  (PL)  sinhaL  cospL]  -  [(aL)2  + 
(pL)2  sinhaL  sinpL  (5.34) 

C3  sinh(aL)  -  ci  sin(PL)  =  {T2  [(aL)  coshaLsinPL+pLsinhaLcosPL]  - 
[(aL)2+(PL)2]sinhaLsinPL}  (aL)2  (5.35) 

C3  cosh(aL)  -  C2  sin(PL)  =  pL  coshaL  [T2  cosPL-PL  sinpL]  -aL  sinpL  [aL 
coshaL+T2  sinhaL]  (5.36) 

C3  cosh(aL)  -  C2  sin(PL)  =  {T2  [(pL)  coshal  cospi  -  (aL)  sinhaL  sinpL]  -  [(al)2  + 
(PI)2  ]  coshaL  sinpL}.  Ti(aL)  (5.37) 

C2  sinhaL  -  ci  coshaL  =  (aL)  sinhal  [aLcoshaL  +T2sinhaL] 
-aL[coshaL[aLsinhaL+T2CoshaL]  (5.38) 

Collecting  and  rearranging  all  the  terms  we  determine  frequency  equation  as  follows- 

-  Ti  T2  (aL)  (pL)  +  Ti  T2  (pL)  [(aL)  coshaL  cospL  +  (pL)  sinhaL  sinpL]  +  Ti  (pL) 
[(aL)2  +  (PL)2]  sinhaL  cospL-  T2  [(aL)2  +  (PL)2]  [(aL)  coshaL  sinaL  -  (pL) 
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sinhaL  cospL]  -  [(aL)2  +  (PL)2]2  sinhaL  sinPL  -T2[(aL)2+(PL)2].[aLcoshaLsinPL- 
pLsinhaLcosPL]-[(aL)2+(pL)2].  [(aL)2+(PL)2]  sinhaL  sinPL-Tl.T2  (aL.pL) 

(5.39) 

=  {-  2(aL)  (pL)  +(aL)  (pL)  coshaL  coshpL  +(PL)2  sinhaL  sinPL  +(aL)  (PL) 
coshaL  cosPL(aL)2  sinhaL  sinPL]  T1T2  +  Ti  {(PL)  [(aL)2  +  (PL)2]  SinhaL  cosPL 
-(aL)[(aL)2  +  (PL)2]  coshal  sinPL}  -T2[(aL)2  +  (pL)2]  [(aL)  coshaL  sinpL  -  (PL) 
sinhaL  cosPL]  -  [(aL)2  +  (PL)2]  sinhaL  sinpL  =  0  (5.40) 

-Ti  T2  {[(aL)2(PL)2]  sinhaL  sinpL  +  2  (aL)  (PL)  (1-coshaL  cosPL)}  +  Ti  [(aL)2  + 
(PL)2]  [aL  coshaL  sinpL  -(pL)  sinhaL  cospL]  +  T2  [(aL)2  +  (PL)2]  [aL  coshaL 
sinpL  -  (pL)  sinhaL  cospL]  +  [(aL)2  +  (PL)2]2  sinhaL  SinpL  =  0  (5.41) 

{[(aL)2  +  (pL)2]2  +  [(aL)2  -  (pL)2[  T1T2}  sinhaL  sinpL  +  2(aL)  (pL) 

T 1T2  (1-coshaL  cospL)  =  0  (5.42) 

The  final  frequency  equation  is  given  by- 

{[(aL)2  -(PL)2]  +  — —  [(aL)2  +  (PL)2]}  sinhaL  sinpL  +  2aL.  pL  (1-coshaL  cospL) 
T|T2 

\  1  ^ 

+ - I - [(aL)2  +  (PL)2  [  (aL)  coshaL  sinpL  -  (PL)  sinhaL  cospL]  =  0  (5.43) 

V^2  hj 

If  T 1  =  T2  =  oc  we  get 

[(aL)2  -  (PL)2]  sinhaL  sinpL  +  (pL)  (1-  coshal  cospL)  =  0 

[(PL)2  -  (aL)2]  sinh(aL)  sinh(PL)  +  2(aL)  (PL)  [cosh(aL)  cos  (pL)  -  1]  =0 


(5.44) 
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5.3  Derivation  of  the  Mode  Shape  Expression  for  Transverse  Vibrations  of 
Elastically  Restrained  Single  Bellows  Expansion  Joint 

We  know 

y  =  A  sinh(  aL)^  +  B  cosh(  aL)^  +  C  sin(  pL)^  +  D  cos(pL)^ 

(5.45) 


Where 


The  boundary  conditions  for  a  single  bellows  are  given  by- 


y  =  0;  at^  =  0 


cA=_t  3y 

d^2  1  dt, 


^=-T,^at^0^T,= 


'R,Ln 

v"eTv 


y  =  O  ;  at  £,  =  L 

=  T2  ^  at  £  =  L 

.  2  z  ~ 


d't; 


T 

1  2 


'  R  2 1  -*  " 

l  El  J 


—  =  (aL)  [A  cosh(aL)£,  +  BSinh  (aL)^]+(pL)[Ccos(pL)^-DSin(pL)^]  (5.46) 

dt, 


^4=(aL)2  [ASinh  (aL)^ +BCosh(aL)^]-(pL)2  [CSin(pL^+  DSin (pL)^]  (5.47) 


At  £,  =0;  y=  0;  B  +  D=0;  (aL)2  [B]-(pL)2  [d]=-T,  [(aL)A  +  (pL)c] 

At  £,  =l;(aL)2  [Asinh  (aL)+Bcosh(aL)]-(pL)2  [Csin  pL  +  DcospL]  + 

=T2  {(aL)[Acosh(aL)+Bsinh  (aL)]+(pL)[Ccos(pL)-DSin(pL)] }  (5.49) 
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B  +  D=  0 

D  =  -B 

B(aL)2  -  D(PL)2  =T,  [A(aL)A  +  (f3L)c] 

(5.50) 

A  sinh  (aL) + B  cosh(aL)  +  Csin  (pL)+  Deos(pL)=0 

(5.51) 

At  £,=l;(aL)2  [Asinh  (aL)+Bcosh(aL)]-(pL)2  [Csin  pL  +  DcospL] 

=T2  {(aL)[Acosh(aL)+Bsinh  (aL)]+(pL)[Ccos(pL)-Dsin(pL)]} 

(5.52) 

“=-i 

D 

(5.53) 

[(aL)]A + [(PL)]C= + ^  [(aL)2 + (PL)2  ] 

1 

(5.54) 

[sinh  (ocL)]A + [sin  PL]C = -  D  [cos  PL + cosh  aL] 

(5.55) 

|  (aL)^sinh  (aL)  A-  (pL)^sinPL  C -[T2(aL)coshaL]A-[T2(pL)cos 

,pL]c| 

=  {t2  (aL)[sinh  (aL)]B  -  [T2 (pL)sin  pL]D - B  [(aL  f  cosh  aLj+  D  [(pL)2  cos pLj }  (5.56) 

A{(aL)[(aL)sinh  aL-T2  coshaL]}-C{(pL)[(pL)sin  PL+T2  cosPL]} 

=  D  [-T2(aL)sinh  aL-T2(pL)sin  pL  +  (pL)2  cosPL  +  (aL)“  coshaL] 

(5.57) 

A(aL)[(aL)sinh  aL-T2  coshaL]  -C  (pL)[(pL)sin  pL+T2  cosPL] 

=  D  {(aL)2  coshaL+(pL)“  cosPL-T2  [aLsinh  aL+pLsin  pL]} 

(5.58) 

104 


“=-i 

D 


(aL)  A  +  (pL)C  =  —  [(aL)2  +  (pL)2 
Ti 


(5.59) 


[sinh  aL]A  +  [sin  PL]C  =  -D[cospL  +  cosh  aL] 


(5.60) 


Multiplying  (5.59)  by  sinpL  and  (5.60)  by  pL  and  canceling  the  like  terms  we  get, 


5.59xsin  pL  :  [(aL)  sin  pL]A  +  [(pL)  sin  pL]C  =  °  [(aL)2  +  (pL)2 


5.60x  PL  :  [(PL)  sin  olL]A  +  [(pL)  sin  pL]C  =  ±D(pL)  [cosh  aL  +  cos  PL] 


On  subtracting,  we  get- 


A[c/L  sin  pL  -  pLsinh  aL]  =  D<!  S'n  ^  (al )2  +  (PL)2  +  (pL)[cosh  aL  +cosh  pL]  J 


(5.61) 


A  (aL)2  +  (pL)2  sin  PL  +  Tj  (pL)[cosh  al  +  cos  PL] 
D  Tt[aL  sin  PL  -PL  sinh  aL] 


(5.62) 


To  find  the  ratio  C’=C/D,  we  multiply  (5.59)  by  sin  haL  and  (5.60)  by  aL 


(5.59) x  sinh  aL  :  [(aL)  sinh  aL]A  +[pL  sinh  aL  ]C  =  ^ [(aL)  +  (pL)" 


(5.60)  x  aL  :  [(aL)  sinh  aL]A  +[aLsinh  PL  ]C  =  -D(aL)[cosh  aL  +  cos  PL] 


Subtracting  we  get- 
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-C[(aL)  sin  |3L-((3L)sinh  aL] 

=  D  <  ]  gj^  aL  +  (aL)[coshaL  +  cos 


[f(aL)2+(pL)2' 

sinh  aL  +  T,  (aL)[cosh  aL  +cos  pL]  ] 

l 

aL  sin  pL  -  pL  sinh  aL] 

(5.63) 


(5.64) 


C'= 


B 

D 


(aL)2+(pL)2 

sinh  aL  +  T,  (aL)[coshaL  +cos  pL] 

T, 

aL  sin  pL  -  pL  sinh  aL] 

(5.65) 


A  (aL)9  +  (pL)9  sinh  pL  4-  T:  (pL)[coshaL  4-  COS  J3L] 
D  T[  [aL  sin  PL  -  pi  sinh  aL] 


y(£,)  =  D]  [cos  pL  -  cosh  (aL)^]  +  —  sin  pL^+  —  sinh  aL<; 

I  vDj  vDj  I 


m 

D 


[cos  (5L  \  -  cosh  aL  £,]  +  {  — f - [(aL)  +(P^)  ] - j  ^  pjj:  _ 

[  T,  [  aL  sin  f)L  -  (3L  sinh  aL]J 

[(aL)2  +(pL)‘]sinh  aL  +  T,  (aL)[cosh  aL  +  COS  pL]l  .,  .  . 

T,  [  aL  sin  (3L  -  (3L  sinh  aL] 


(aL)A  +  (PL)C  =  ^-[(aL)2  +  (PL)2 

A(aL)[(aL)sinh  aL  -  T2  cosh  aL]  -  C(pL)[(pL)sin  PL  +  T2  cos  PL] 
=  D|(aL)2  coshaL  +  (pL)2  cos  pLj-T2[aL  sinh  aL  +  pLsin  PL]} 
(2)  x  [aL  sinh  aL  -  T2  cosh  aL] 

A(aL)[aL  sinh  aL  -  T2  cosh  aL]  +  (pL)c[(aL)sin  aL  +  T2  cos  aL] 
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=  —  [(aL)2  +  (PL)2  JaL  sinh  aL  -  T2  cosh  aL] 

Ti 

(pL)c[aL  sinh  aL  -  T2  cosh  aL]  +  (pL)c[PL  sinh  PL  -  T2  cosh  PL] 

(PL) 


—  [(aL)’  +  (pL)’  JaL  sinh  aL  -  T2  cosh  aL] 

T| 


D|(aL)2  coshaL  +  (pL)"  cosh pLj-T2[aL sinh  aL+pLsin  PL] 


c 

D 


(pL){[aL  sinh  aL  +  PL  sin  pL]  +  T2  [cos  PL  -  cosh  aL]} 


=  [(aL)"  +  (pL)2  JaL  sinh  aL  -  T2  cosh  aL] 

-  Tj  {[(aL)2  cosh  aL  +  (pL)2  cosh  pL  J  -  T2  [  aL  sinh  aL  +  pL  sin  PL]} 


D 


[(aL)2  +  (pL)2]aL  sinh  aL  -  T2  cosh  aL]  -  T,  [(aL)2  cosh  aL  +  (pL)2  cosh  Pl] 
+  T[To  [  aL  sinh  aL  +  PL  sin  PL] 


T,(pL)[aLsinh  aL+  pLsin  PL]+  T,T,(cospL  -  coshaL) 


|[(aL)2  +  (PL)2  ](aL)  sinh  aL  -  T2  [(aL)2  +  (PL)2  ]cosh  aL  -  T, 

(aL)2  cosh  aL  +  (PL)2  cos  PL  +  | 
T, T2  [aL  sinh  aL  +  PL  sin  PL]  i 

[  D  )  {T,  (pL)[aL  sinh  aL  +  PL  sin  PL]  +  T,T2  [cos 

;  PL  -  cosh  aL]} 

2  x  |j3L  sin  PL  +  T2  cos  PL] : 

(aL)A  [|3L  sin  PL  +  T2  cos  pL]  +  (pL)c[pL  sin  PL  +  T2  cos  pL] 

=  —  [(aL)2  +  (pL)2  ]j3L  sin  pL  +  T2  cos  pL] 

Ti 


(aL)A[aL  sinh  aL  -  T2  cosh  aL]  -  (j3L)c[pL  sin  PL  +  T2  cos  PL] 
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=  Dj[(aL)2  cosh  aL  +  (pL)2  cosPlJ-  T2  [aLsinh  aL  +  pLsin  PL]} 
Adding- 

(aL)A{[aL  sinh  aL  +  PL  sin  PL]  +  T2  [cos  PL  -  cosh  aL]} 


=  —  [(aL)2  +  (pL)2  ])3L  sin  PL  +  T2  cos  PL  ] 

L 

+  D|(aL)2  cosh aL  +  (pL)2  cos pL]-T2[aLsinli  aL  + PLsin  PL]] 

f  \\ 

—  KT,  (aL)[aL  sinh  aL  +  PL  sin  PL]  +  TjT2  (aL)[cos  PL  -  cosh  aL]} 

vDJ 

_  |[(aL)2  +  (PL)2  JpL  sin  PL  +  T2  cos  PlJ  +  T,  [(aL)2  cosh  aL  +  (PL)2  cosPl]} 
T,T2[aL sinh  aL  +  PLsin  PL]  j 


|[(aLp  +  (pL)2]j3Lsin  PL  +  T2  cospL]  +  T,  [(aL)2  cosh  aL  +  (PL)2  COS  Pl]- 

A  [t,L [aL sinh  aL  +  PLsin  PL] _ ; 

D  T,  (aL){[aL  sinh  aL  +  PL  sin  PL]+T2  [cos  PL  -  cosh  aL]} 

Mode  shape  can  be  therefore  be  written  as- 


D 


A' sinh  (aL)£,  -  cosh(aL)£,  +  C'sin  (PL)£,  +  cos(pL)^ 


Or  can  be  written  as- 


"yfe)~ 

D 

=  [C  sin  (pL)c,  +  A'  sin  (aL)c,]  +  [cos(pL)c,  -  cos(haL)c,] 
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5.4  Results  and  Discussion 

Equation  (5.44)  is  a  closed-form  general-purpose  transcendental  frequency 
equation  for  finding  out  the  natural  frequencies  of  transverse  vibrations  of  single 
bellows  expansion  joints  restrained  against  rotation  on  either  ends  and  defined  by 
non-dimensional  rotational  restraint  parameters  such  as  T iand  T2. 

By  providing  different  values  of  non-dimensional  parameters  ‘c’  which 
depends  on  internal  pressure  ‘P’,  bellows  moment  of  inertia  ‘J’  and  Ti  &  T2  that 
depend  on  rotational  spring  stiffness  Ri  and  R2,  the  transcendental  equation  can  be 
solved  for  values  of  ‘  X  ’  for  various  modes  of  vibration  ‘N’. 

The  bisection  method  of  trial  and  error  is  utilized  in  this  study  in  obtaining 
various  values  of  frequency  parameter  ‘  X  ’  for  a  given  value  of  c,  Ti  and  T2. 

To  show  the  effect  of  rotational  restraint  parameters  T 1  and  T2  on  frequencies 
of  vibration  of  bellows,  a  single  bellow  with  length  L  =  0.0693m,  mass  moment  of 
inertia  per  unit  length,  J  =  0.001 153kgm,  flexural  stiffness  El  =  5.078Nm  and  total 
mass  m  tot  =  5.138kg/m  is  considered.  The  above  dimensions  of  bellows  are  same  as 
those  considered  by  Jakubauskas.V.F. 

For  T=Ti=T2=0,  the  general  frequency  equation  (5.44)  reduces  to  a  case  of 
single  bellows  with  simply  supported  ends.  Similarly,  for  Ti=T2=cc,  the  equation 
reduces  to  a  single  bellows  with  fully  fixed  and  elastically  restrained  end  condition. 
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The  frequency  equation  derived  here  is  exact  and  can  be  applied  to  all  types  of  single 
bellows  with  unequal  rotational  restraints  on  either  end. 

According  to  [1],  the  maximum  internal  pressure  in  bellows  is  given  by  - 

(5.45) 

El  =  0.25  k.P.Rm2  (5.46) 

Substituting  (5.45)  and  (5.46)  and  the  dimensional  and  material  property  data  of  the 
bellows  expansion  joint  given  above  into  expressions  (5.6)  and  (5.7),  we  get  - 

C  =  V(61 6.49  +  0.0001 135co2)  (5 .47) 

A. =1 .0029V©  (5.48) 

In  order  to  study  the  effect  of  variation  of  internal  pressure  on  the  lateral 
natural  frequency  is  investigated.  The  first  four  modes  of  transverse  natural 
frequencies  are  obtained  for  internal  pressures  of  P=0.0MPa  and  P=166.0MPa  and  are 
given  in  Tables  5.1  &  5.2.  The  same  are  graphically  represented  in  figures  5.2  &  5.3 
respectively.  The  effect  of  rotational  restraint  on  the  natural  frequency  is  studied  for 
an  initial  value  of  0.01  to  a  maximum  of  1010  and  is  assumed  to  be  equal  on  either 
end. 

Firstly,  it  is  seen  that  as  internal  pressure  of  the  bellows  increases  the 
transverse  natural  frequency  decrease.  For  example  at  P=0.0MPa  and  mode  of 
vibration  N=l,  natural  frequency  obtained  is  ©=1690.477  radians/s  and  at 
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P=166.0MPa  and  N=l,  and  (0=1069.4297  radians/s  a  drastic  drop  by  about  37%.  The 
same  is  true  for  even  higher  order  of  modes  of  vibration  N=2,  3  &  4. 

It  is  also  seen  that  as  the  value  of  rotational  restraint  T  -  increases  from  0.01 
to  1010,  the  frequencies  tend  to  increase  for  all  modes  of  vibration.  As  the  rotational 
stiffness  value  T— »cc,  the  frequency  value  increase  by  about  54%  for  N=1  at 
P=0.0MPa  and  by  69%  for  N=1  at  P=166.0MPa  respectively.  However,  it  is  observed 
that  there  is  no  change  in  frequency  and  it  becomes  constant  from  T=104  onwards. 

Table  5.1  gives  a  comparison  of  the  exact  frequency  solutions  obtained  using 
the  bisection  method  for  a  single  bellows  fixed-fixed  at  both  ends  with  rotational 
restraint  vis-a-vis  to  the  results  presented  in  the  reference  [5]. 


Table  5.1  Comparison  of  Frequency  Solutions  for  T=a  and  P  max 


Mode# 

Exact  [44]  (o, 

rad/s 

[28]  to,  rad/s 

Error,  % 

1 

3400.159 

3400.334 

0.005 

2 

6890.181 

6890.356 

0.002 

The  first  two  natural  frequencies  are  calculated  for  the  bellows  data  as  given 
above  using  the  frequency  equation  that  is  derived.  It  is  seen  that  the  error  %  of  the 
frequency  obtained  from  the  exact  method  is  about  0.005%.  Therefore,  as  the 
percentage  of  error  is  well  within  the  engineering  accuracy  it  is  found  to  be  precise 
enough  to  estimate  the  natural  frequency  of  single  bellows  expansion  joint  using  this 


method. 
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Table  5.2:  Natural  Frequencies  at  P=O.OMPa  for  T=Ti=T2  &  N=1  to  4 


T 

N=1 

N=2 

N=3 

N=4 

0.01 

1690.477 

4845.901 

8103.03 

14431.45 

0.1 

1692.61 

4847.431 

8104.168 

14432.18 

1 

1713.655 

4862.628 

8115.502 

14439.46 

10 

1899.581 

5005.034 

8224.826 

14511.13 

102 

2751.07 

5851.544 

12067.52 

15109.97 

103 

3523.861 

6909.163 

13444.12 

16590.34 

104 

3666.924 

7123.951 

13808.88 

17018.88 

105 

3682.946 

7154.077 

13848.72 

17066.6 

106 

3683.946 

7156.511 

13852.74 

17071.41 

107 

3684.102 

7156.754 

13853.16 

17071.89 

108 

3684.118 

7156.778 

13853.19 

17071.94 

109 

3684.119 

7156.781 

13853.19 

17071.95 

1010 

3684.119 

7156.781 

13853.19 

17071.95 

Fig  5.2:  Transverse  Vibrations  Single  Bellows 
(Lateral  Frequency  for  P=0.0MPa  &  N=1  to  4  &  T=T1=T2 
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Table  5.3:  Natural  Frequencies  at  P=166.0MPa  &  T=Ti=T2  &  N=1  to  4 


T 

N=1 

N=2 

N=3 

N=4 

0.01 

1069.425 

4467.753 

7828.299 

11080.62 

0.1 

1072.798 

4469.413 

7829.476 

11081.53 

1 

1105.697 

4485.891 

7841.207 

11090.6 

10 

1375.915 

4639.817 

7954.354 

11179.13 

102 

2408.615 

5539.906 

8754.995 

11868.26 

103 

3235.598 

6636.899 

10063.38 

13270.83 

104 

3382.595 

6862.794 

10364.6 

16872.66 

105 

3398.391 

6887.421 

10397.61 

16920.74 

106 

3399.982 

6889.904 

10400.94 

16925.59 

107 

3400.142 

6890.153 

10401.27 

16926.08 

108 

3400.142 

6890.179 

10401.3 

16926.13 

109 

3400.16 

6990.181 

10401.31 

16926.13 

1010 

3400.16 

6890.181 

10401.31 

16926.13 

Fig  5.3  Transverse  Vibrations  Single  Bellows 
(Lateral  Frequency  for  P=166.0MPa,  T=T!=T2  and  N=1  to  4) 
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The  mode  shapes  are  obtained  for  single  bellows  by  varying  the  rotational 
restraint  parameter  T  from  0.01  to  10  for  both  first  and  second  modes  of  vibration 
(N=l&2)  respectively.  The  ratio  y/D  is  found  out  for  c,  values  between  0  to  0.5.  The 


response  for  the  first  half  of  the  mode  is  obtained  and  is  found  to  be  in  symmetry  for 
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the  next  half  mode.  The  effect  of  rotational  restraint  T,  on  the  response  of  single 
bellows  with  elastically  restrained  ends.  Firstly,  we  observe  the  effect  of  variation  of 
rotational  restraint  T  on  response  of  bellows.  It  is  clearly  seen  that  the  response  shows 
upward  trend  as  L,  increases  from  0  to  0.5  for  the  first  half  cycle  of  vibration.  It  is 
observed  that  at  T=0.01,  y/D  is  -175  and  at  T=10.0,  y/D  is  -83  a  rise  by  about  50%. 


Table  5.4  Ratio  y/D  for  different  values  of  T  and  N=1 


y/D,  T=0.01 

y/D  T=0.1 

y/D,  T=1 

y/D,  T=10 

0 

0 

0 

0 

0 

0.1 

-175 

-44 

-62 

-83 

0.2 

-351 

-87 

-124 

-165 

0.3 

-526 

-133 

-186 

-248 

0.4 

-701 

-177 

-248 

-330 

0.5 

-877 

-221.5 

-310 

-412 

Table  5.5  Ratio  y/D  for  difl 

ferent  values  of  T  and  N=2 

S 

y/D,  T=0.01 

y/D,  T=0.1 

y/D,  T=1 

y/D,  T=1 0 

0 

0 

0 

0 

0 

0.1 

4538 

4462 

4536 

8330 

0.2 

9076.4 

8925 

9071 

16659 

0.3 

13614.5 

13387 

13107 

24989 

0.4 

18512.7 

17849 

18142 

33318 

0.5 

22690.9 

22312 

22678 

41647 

Tables  5. 4. and  5.5  presents  values  of  y/D  by  varying  the  rotational  restraint  T  and 


^=0  to  0.5.  The  same  are  plotted  and  shown  in  Fig.  5.4  and  5.5  respectively. 


The  response  changes  signs  for  the  second  mode  of  vibration.  The  response  for  T=10 


and  N2  is  almost  double  when  compared  to  the  response  at  T=1.0  at  c=0.5. 
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Fig  5.5  Mode  Shape  of  Single  Bellows  with  Elastically  Restrained  Ends 
for  varying  values  of  rotatioanl  restraint  T=T1=T2,  P=().0MPa 

&  N=2 


zeta=x/L 
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CHAPTER  6 

FINITE  ELEMENT  ANALYSIS  OF  TRANSVERSE  VIBRATIONS 
OF  SINGLE  BELLOWS  EXPANSION  JOINT  RESTRAINED 

AGAINST  ROTATION 

6.1  Theoretical  Background 

This  chapter  presents  the  application  of  finite  element  method  to  the 
determination  of  the  transverse  frequencies  of  single  type  of  bellows  restrained 
against  rotation  on  either  end.  The  aim  of  the  work  here  is  to  model  the  U  shaped 
bellows  using  beam  elements.  The  effect  of  rotatory  inertia  on  the  natural  frequency 
is  included  in  the  beam  element  matrices. 

Results  obtained  from  the  finite  element  method  are  then  compared  with  the 
frequencies  obtained  by  exact  method  for  the  first  two  modes  of  vibration. 
Experimental  results  are  also  used  for  comparison  purpose.  The  effects  of  variation  of 
the  rotational  restraint  parameter,  internal  pressure  and  velocity  on  the  natural 
frequencies  are  studied.  The  geometrical  dimensions  assumed  in  the  present  work  are 
the  same  as  considered  by  Jakubauskas.V.F. 

6.2  Characteristics  of  Bellows 

The  various  geometrical  dimensional  parameters  of  U-shaped  bellows  are 
given  in  Fig  1.2.  Where  n  is  the  meridional  radius  of  the  convolution  root,  T2  is  the 
meridional  radius  of  the  convolution  crown  and  h  is  the  convolution  height.  Rm  is  the 
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mean  radius  of  the  bellows,  that  is,  the  distance  from  the  bellows  centerline  to  mid 
convolution  height,  and  t  is  the  bellow  material  thickness.  It  is  assumed  that  t «r  1 ,  r2 
and  h  «  Rm.  With  ‘N’  as  number  of  convolutions  the  total  length  of  the  bellows  is  L  = 
2(ri+r2)  N.  Therefore,  with  these  assumptions  bellows  are  considered  as  equivalent 
pipe  /  bar  of  radius  Rm  and  wall  thickness  t. 

6.3  Theory  of  Free  Vibrations  of  Bellows 

The  matrix  equation  for  the  free  vibration  of  bellows  can  be  written  as  - 

[M]  { q }  +  [K]  { q }  =0  (6.1) 

Where 

{q}  -Generalized  coordinates 

[M]  -  Mass  matrix 

[K]  -  Elastic  stiffness  matrix 

6.4  Finite  Element  Formulations  of  Elastic  Stiffness  Matrix 

The  strain  energy  U  of  a  bellow  element  of  length  ‘L’  including  the  effect  of  rotatory 
inertia  is  given  by  the  relation- 

i  i 

U=  Vi  El  I(d2w/dx2)2.dx  -  Vi  P7iRm2!(dw/d  x)2.dx  (6.2) 

o  o 

Where 

w  -  deflection  of  the  bellow 

El  -bending  stiffness  of  bellow 

P  -internal  pressure  and 

R  m  -mean  radius  of  bellows 

Now  on  non-dimensionalizing  and  substituting  the  following  relations,  we  get 
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r)  =  x  /  L  and  cj)  =  y  /  L  the  above  expression  of  strain  energy  becomes  as  follows  - 


i  i 

U=  y2  EI/L  /(dtydrf)  2.dq  -  Vi  PnRm2/  Lj(d(p/dri)  2.dq  (6.3) 

0  0 


Assuming  a  cubic  polynomial  expression  for  cp  to  be  of  the  form  - 

(p  =  X  ar  rj  r  (6.4) 

Now  substituting  in  to  equation  6.1)  and  replacing  the  coefficient  of  ar  (r=0, 1,2,3),  the 
strain  energy  expression  becomes  - 

U=  Vi  EI/1  {£,}  r  [K]  (6.5) 

Where  t,  is  the  degree  of  freedom  (DOF) 

.-.  [K]  =  [k  e]  -  A  2  [k  v]  (6.6) 

Where 

A2  =  Pn  Rm2  L2  /  El 
[K]  -Elastic  stiffness  matrix 

[ke]  -Velocity  stiffness  matrix  and  are  defined  as  follows  - 


12  —  6L  -12  6L 
2 El  6L  4l2  -6L2l2 

1 ~Y  -12  -6L  12  -6L 

6L  -2l2-6L  4l2 


(6.7) 


-36  —  3L  -36-3L 

[k.]=v!  3L  4L2  —  3L  — l2 

L  J  30  -36-36L  36  +  3L 
-3L-3l2  +  3L  4l2 


(6.8) 
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Where  £, T  =  { cpi,  cpf,  cpi+i,  cp’  i+i}  (6.9) 

6.5  Finite  Element  Formulation  of  Mass  and  Rotatory  Inertia  Matrices 

The  kinetic  energy  ‘V’  of  a  bellow  element  of  length  including  the  effects  of  rotatory 
inertia  is  given  by  - 


V=  Vi  p  A  L3  J  (9)  2  dr)  +  Vi  p  J  L  J  (cp’)  2  dr|  (6.10) 

0  0 

The  mass  of  bellows  per  unit  length  is  computed  by  using  the  formula  given  in 
equation  (6.11)- 


M  = 


pb.27i.Rm{7r(rt  +  n)  +  2(h 
2(rt  +  n) 


n  —  r2)}t 


(6.11) 


Where 

pb  -  mass  density  of  bellows  material 

pf  -  mass  density  of  fluid  flowing  through  bellows 

A  -  area  of  cross-section  of  the  bellow 

Now  substituting  for  (p  from  equation  (6.2)  and  replacing  the  coefficient  a  r(r=0, 1,2,3) 
by  the  nodal  coordinates  the  expression  becomes- 

T=  Vz  pAL3  {^}T  +  [M]{^}  (6.12) 

Where 

[M]  =  [Mi]  +  R  k4  [M2] 

(6.13) 
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Where 


R  = 


J 


mtot.  L2 


mt0tL4Q)n 

El 


The  mass  moment  of  inertia  of  the  bellow  per  unit  length  is  given  by- 


(6.14) 


J  =  Jxx  =  Jyy  =  tur 


17  h  1 

\  2— +  0.571 

t.pb  +  —  (2r2  -  t)pf 

11  q  J 

q  J 

(6.15) 


Where  Rm.,  h,  q,  t,  ri  and  ri  are  already  been  defined  earlier 
Now  the  mass  matrices  Mi  and  M2  are  defined  as  follows- 


[M,] 


156-22L  54  -13L 

pAL  22L  4L2  13L  -3l2 

420  -54-13L  156  -  22L 
13L  —  3  L2  —  22L  4l2 


[M2] 


1 

420 


36-3L-36-3L 
3L4l2-3L-l2 
36-3L  36  —  3L 
3L  —  L2  —  3L  4  L2 


(6.15) 


(6.16) 


The  matrix  equation  for  free  vibration  of  single  bellows  is  given  by- 
{ [K]  -  A,  4  [M] }  {£,}  =  0 


Where 

}.4  =  {p  A14P„2}  /  El 


(6.17) 
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The  matrix  eigen  value  equation  is  solved  for  the  bellows  restrained  against  rotation 
using  the  following  boundary  conditions 

6.6  Boundary  Conditions 

At  the  end  ’A’  -  bellows  are  connected  to  a  pipe  nipple,  and  are  considered  to  have  a 
rotational  stiffness  of  Ri  and  R2  at  either  end  as  shown  in  figure  5.1. 


At  the  end  A  (x  =  0) 

X  (x=0) 

=  0  and 

(6.18) 

d2  X  (0) 

ax(0) 

-Ti  - 

(6.19) 

dx  2 

dx 

At  the  end  B  (x  =  L)  the  boundary  conditions  are  given  by  - 
X  (L)  =  0  and  (6.20) 

82  X  (L)  ax  (L) 

- =  x2  -  (6.21) 

dx 2  dx 

Where, 

Ti  =Ri  L  /  El  (6.22) 

T2  =  R2  L  /  El  (6.23) 

6.7  Results  and  Discussion 

A  single  bellow  with  the  following  geometrical  and  physical  parameters  is 
considered  -bellows  length  L  =  0.0693m,  mass  moment  of  inertia  per  unit  length,  J  = 
0.001 153kgm,  El  =  5.078Nm  and  total  mass  m  tot  =  5.138kg/m  respectively.  The 
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geometrical  dimensions  of  bellows  used  are  same  as  that  considered  by 
Jakubauskas.V.F  (28). 

The  matrix  eigenvalue  equation  is  solved  using  the  Jacobi’s  method  for  the 
bellows  restrained  against  rotation  and  using  the  relevant  boundary  conditions.  The 
stiffness  matrix  is  notably  modified  to  include  the  unequal  restraints  against  rotation 
considered  at  either  end. 

The  first  two  modes  of  transverse  natural  frequencies  are  obtained  for 
maximum  critical  pressure  of  P=166.  Firstly,  it  is  seen  that  as  internal  pressure  of  the 
bellows  increases,  the  lateral  mode  natural  frequencies  decrease.  At  P=0.0MPa  and 
mode  number  N=l,  natural  frequency  obtained  is  co  =1690.20rad/s  and  at 
P=166.0MPa  and  N=l,  co=1069.06  radian/s  a  drastic  drop  by  about  37%.  The  same  is 
true  for  even  for  other  modes  of  vibration  N=2  &  3  respectively. 

It  is  also  seen  that  as  the  value  of  rotational  restraint  T  (T=Ti=T2)  -increases 
from  0.01  to  1010,  the  frequencies  tend  to  increase  for  all  modes  of  vibration.  As  the 
rotational  stiffness  value  T— »cc,  the  frequency  value  increases  by  about  54.1%  for 
N=l(P=0.0MPa)  and  by  68.56%  for  N=1  (P=166.0MPa)  respectively.  However,  it  is 
observed  that  there  is  no  change  in  frequency  and  it  becomes  constant  from  T=104 


onwards. 
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Table  6.1:  Natural  Frequencies  for  T=Ti=T2  and  N=1  to3  &  P=0.0MPa 


T 

N=1 

N=2 

N=3 

0.01 

1690.2 

4846.58 

8112.22 

0.1 

1717.17 

4866.0 

8126.83 

1.0 

1927.65 

5028.8 

8252.68 

10 

2821.56 

5938.55 

9104.02 

102 

3544.04 

6944.05 

10343.6 

103 

3669.43 

7138.60 

10608.4 

104 

3682.84 

7159.61 

10637.1 

105 

3684.19 

7161.73 

10640.0 

106 

3684.33 

7161.94 

10640.3 

107 

3684.34 

7161.94 

10640.3 

108 

3684.34 

7161.94 

10640.3 

109 

3684.34 

7161.94 

10640.3 

10'° 

3684.34 

7161.94 

10640.3 

Fig  6.1 :  Transverse  Vibrations  of  Rotationally  Restrained  Single  Bellows 
T=0.01  to  1010,  P=0.0  MPa,  RIX=0. 04673  &  V=0  m/s 


frequency,  rad/s 
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Table  6.2:  Natural  Frequencies  for  T=Ti=  l  2  and  N=1  to  3  &  P=166.0MPa 


T 

N=1 

N=2 

N=3 

0.01 

1069.06 

4468.55 

7837.89 

0.1 

1111.2 

4489.69 

7852.98 

1.0 

1414.32 

4665.5 

7983.22 

10 

2486.92 

5631.27 

8862.16 

102 

3256.51 

6672.85 

10131.0 

103 

3385.25 

6871.85 

10400.1 

104 

3398.95 

6893.31 

10429.2 

105 

3400.33 

6895.47 

10432.1 

106 

3400.47 

6895.69 

10432.4 

107 

3400.48 

6895.71 

10432.5 

108 

3400.48 

6895.71 

10432.5 

109 

3400.48 

6895.71 

10432.5 

1010 

3400.48 

6895.71 

10432.5 

Fig  6.2:  TransverseVibrations  of  Rotationally  Restrained  Singles 
Bellow  for  T=0.01  to  1 010  ,  P  =1 66.0  MPa  &  V=0  m/s, 
RIX=0. 04673 
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Rotational  Stiffness,  T 
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Table  6.3  gives  a  comparison  of  frequencies  obtained  by  exact  method  with 
frequencies  obtained  using  the  finite  element  method  for  a  single  bellows  elastically 
restrained  at  both  ends. 


Table  6.3  Comparison  of  Frequency  for  T-cc  and  Pmax=166.0MPa 


Mode 

# 

[28]  q, 
rad/s 

Exact  co, 
rad/s  [44] 

FEA  [45] 
o>,  rad/s 

1 

3400.334 

3400.159 

3400.48 

2 

6890.356 

6890.181 

6895.71 

The  first  two  natural  frequencies  are  calculated  for  the  bellows  using  the  finite 
element  method.  It  is  seen  that  the  results  obtained  by  the  finite  element  method  are 
found  to  be  in  close  agreement  with  the  exact  method.  The  percentage  error  in  the 
frequency  obtained  is  less  is  precise  enough  to  estimate  the  natural  frequency  of 
single  bellows  expansion  joint. 


The  effect  of  variation  of  velocity  of  flow  at  V=  1.0,  5.0  and  10  m/s  on  the 
first  three  frequencies  is  also  studied.  It  is  seen  that  for  a  constant  value  of  rotational 
stiffness  T  =  1010  the  frequency  increases  with  increase  in  the  velocity  of  flow  and 
also  with  increase  in  mode  numbers. 


Table  6.4  Frequency  at  different  velocities  and  constant  1=  101(> 


Mode 

Number 

V=1.0m/s 

V=5.0m/s 

V=10.0m/s 

Frequency,  Hz 

Frequency, 

Hz 

Frequency, 

Hz 

N=1 

0.08 

2.0 

7.4 

N=2 

0.16 

4.0 

15.5 

N=3 

0.25 

6.0 

23.7 
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It  is  seen  that  the  effect  of  variation  of  pressure  and  velocity  of  flow  on  the 
natural  frequency  is  marginal  because  the  dimensions  of  bellows  given  in  reference 
(4)  are  small.  The  effect  of  variation  of  pressure  and  flow  velocity  on  the  transverse 
natural  frequencies  for  modes  N=l,  2  &  3  is  significant. 


In  order  to  demonstrate  this  clearly  a  large  dimensioned  bellow  is  considered 
that  has  the  following  dimensions  -Db=  1.2m,  Rm  =  0.6395m,  L=  0.254m, 
EI=5449Nm2,  J=1485.36kgm  and  mtot=  912.13kg/m  respectively. 

Firstly,  the  effect  of  velocity  of  flow  for  V=1.0  m/s  and  5.0  m/s  is  studied  on 
the  natural  frequency.  It  is  found  that  for  a  constant  value  of  T— >oc  and  velocity  of 
flow  increasing  the  natural  frequency  decreases.  It  is  also  observed  that  the  frequency 
increases  with  increase  in  the  modes  of  vibration  N=l,  2  &  3.  Table  6.5  shows  the 
frequencies  in  radians/s  obtained  for  different  velocities  of  flow  of  liquid  inside  the 
bellow. 


Table  6.5  Frequency  in  radian/s  at  T=oc  &  Varying  Velocities 


Velocity 

V=1.0 

V=2.0 

V=3.0 

V=4.0 

V=5.0 

N=1 

46.71 

11.53 

5.054 

2.8 

1.76 

N=2 

67.38 

16.75 

7.39 

4.13 

2.62 

N=3 

94.77 

23.61 

10.46 

5.867 

3.742 
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Fig  6.3:  TransverseVibrations  of  Single  Bellows  Expansion 
Joint  for  V  =  l.Om/s  and  N=l,2  &3  &  T=infinity 


1.E-02  1.E+01  1.E+04  1.E+07  1.E+10 


Rotational  Stiffness,  T 


- N  =  1 

- N=2 

-  -  N=3 


Fig  6.4:  Transverse  Vibrations  of  Single  Bellows  Expansion  joint 
for  T=infinity  and  V=5.0m/s,  N=l,  2  &3 


Rotational  Stiffness,  T 


u 

u 


- N=1 

-  -  -  N=2 


Similarly,  the  effect  of  variation  of  internal  pressure  on  natural  frequency  is 
investigated  for  the  same  geometrical  dimensions  of  bellows.  It  is  found  that  for 
T— He  and  pressure  increasing  from  lO.OMPa  to  150.0MPa,  the  transverse  natural 
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frequencies  for  the  first  three  modes  of  vibration  decreases.  Table  6.6  shows  the 
frequencies  obtained  at  different  internal  pressures  of  bellows. 


It  is  seen  that  for  P=10.0MPa  &  N=l,  the  percentage  increase  in  frequency  for 
T=0.01  to  T=106  is  50%  and  then  the  frequency  remains  constant  as  T  approaches 
infinity.  At  same  pressure  of  lO.OMPa  and  N=2,  the  percentage  increase  in  frequency 
for  SR=0.01  to  106  is  30%  and  for  N=3  the  percentage  increase  in  frequency  is  25%. 
Similarly  for  a  maximum  pressure  of  150.0MPa  and  N=l,  the  percentage  increase  in 
frequency  for  T=0.01  to  106  is  69%.  At  same  pressure  of  160.0MPa  and  N=2,  the 
percentage  increase  in  frequency  for  T=0.01  to  106  is  33%  and  for  N=3  the 
percentage  increase  in  frequency  is  26%. 


Table  6.6  Frequency  in  radian/s  at  T=oc  &  Varying  Pressures 


Pressure 

P=10MPa 

P=50MPa 

P=100MPa 

P=150MPa 

N=1 

105.35 

20.58 

9.97 

6.43 

N=2 

151.41 

29.94 

14.75 

9.69 

N=3 

212.60 

42.28 

20.98 

13.81 

Fig  6.5:  TransverseVibrations  of  Single  Bellows  Expansion  Joint 
for  T=infinity  and  P=1 0.0  MPa,  N=1 ,2&3 
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150  ' 


100  o5 


50 


1.E+10 


Rotational  Stiffness,  T 


- N=1 

-  -  N=2 
-  N=3 


Fig  6.6:  Transverse  Vibrations  of  Single  Bellows  Expansion  Joint 
for  T=infinity,  P=50MPa  and  N=1 ,2  &3 


Fig  6.7:  Transverse  Vibrations  of  Single  Bellows  Expansion  Joint 
for  T=infinity,  P=1 00.0  MPa  and  N=1 ,2&3 


-  -  -  N=2 

-  N=3 


Rotational  Stiffness,  T 


Frequency,  radian/s 
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Fig  6.8:  Transverse  Vibrations  of  Single  Bellows  Expansion  Joint 
for  T=infinity  and  P=1 50.0  MPa  and  N=1,2  &3 


Therefore,  the  finite  element  method  developed  gives  fairly  accurate  results 
and  is  in  close  agreement  to  the  exact  and  experimental  results.  The  work  leads  to  a 
logical  conclusion  that  is  as  follows- 


i.  As  the  internal  pressure  of  the  bellow  is  doubled  the  transverse  frequency  of 
vibrations  for  the  first  three  mode  numbers  is  reduced  by  nearly  50%. 


ii.  It  is  seen  that  as  the  flow  velocity  is  doubled  the  frequency  decreases  by  about 


4  times  for  all  the  first  three  modes  of  vibration. 
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CHAPTER  7 

TRANSVERSE  VIBRATIONS  OF  UNIVERSAL  EXPANSION 
JOINT  RESTRAINED  AGAINST  ROTATION 
(LATERAL  MODE) 

7.1  Universal  Expansion  Joint 

A  Universal  double  bellows  type  of  expansion  joint  has  two  single  bellows  in 
series  and  separated  between  with  a  short  pipe  spool.  Most  studies  of  transverse 
vibration  of  these  expansion  joints  consider  the  classical  fixed-fixed  end  conditions 
for  investigation  of  transverse  mode  shapes  and  natural  frequency  [46].  However,  in 
real  practice  the  bellow  is  subjected  unequal  rotations  at  the  ends.  Hence,  the  work 
aims  at  investigation  of  transverse  vibrations  of  double  bellows  expansion  joint  by 
considering  elastically  restrained  ends  against  rotation. 


Fig  7.1  Mathematical  Model  of  Double  Bellows  Expansion  Joint 


(Lateral  Mode) 
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A  theoretical  model  as  shown  in  Figure  7.1  is  developed  based  on  the 
Timoshenko  beam  equation  for  lateral  mode  of  vibration  where  L  is  the  length  of 
bellow  and  ‘a’  is  the  length  of  pipe.  It  neglects  the  effect  of  shear  and  includes  the 
effect  of  rotatory  inertia.  An  exact  frequency  equation  is  derived  for  the  transverse 
vibrations  of  double  bellows  expansion  joint  that  includes  the  variation  of  rotational 
restraint  parameter,  T. 

The  bisection  method  is  used  in  computing  the  natural  frequencies.  They  are 
then  compared  with  bellow  showing  the  effects  of  variation  of  elastic  restraints  and 
internal  pressure  on  the  first  two  modes  of  vibration. 

7.2  The  Differential  Equation 

The  general  form  of  the  differential  equation  of  vibration  of  bellows  for  single  or 
double  bellows  expansion  joint  s  given  by 


El 


d4  w 

ax4 


+  PTC  Rm2 


d2  w 


j  g4w 

dx2  dt2 


02w 

+  ITTot - — 

dt2 


(7.1) 


Where  El  is  the  bending  stiffness,  P-internal  pressure,  J-mass  moment  of  inertia  per 
unit  length,  mtot  -  total  mass  of  bellows  per  unit  length  includes  bellows  material 
mass  and  fluid  mass,  x-axial  coordinate,  Rm-  is  the  mean  radius  of  bellow,  w- 


deflection  and  t-time. 
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Using  the  technique  of  separation  of  variables,  the  lateral  deflection  of  the  bellows 
axis  ‘w’  can  be  expressed  as- 


w  (x,  t)  =  X  (x).  T  (t)  (7.2) 

T  (t)  =  A  e  itot  (7.3) 


Differentiating  the  above  equations  and  substituting  into  differential  equation  (7.1) 
we  get- 

a4X  (P7iRm2  + J.w2)  52X  2  mtot 

- — l — - - - - - w  - A  =  U 

ax4  ei  ax2  ei 


The  general  solution  of  the  equation  is  given  by- 

X(x)  =  A  sinh  ax  +  B  cosh  ax  +  C  sinPx+  D  cosPx 

The  first  three  derivatives  of  equation  (7.8)  are  given  as  follows- 


(7.6) 


(7.7) 


(7.8) 


^  =  Aa  cosh  ax  +  pa  sinh  ax  +  CP  cos  P  +  Dp  sinPx 

dx 


(7.9) 
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d2(x) 

dx2 


=  Aa2  sinh  ax  +  Ba2  cosh  ax  -  Cp2  sinPx  -  Dp2  cosPx 


(7.10) 


d3(x) 

dx3 


=  Aa3  cosh  ax  +  Ba3  sinh  ax  -  Cp3  cosPx  +  Dp3  sinPx 


(7.11) 


cx  =  V-  O2  -+-  a/o4  -+-  >0 

Where  the  roots  of  the  equation  are  a  &  P  and  their  values  are  given  by  - 

(7.12) 

P  =  Vc2  +  Vc4  +  7C 

(7.13) 

Where  A,  B,  C,  D  are  arbitrary  constants 

7.3  Derivation  of  Frequency  Equation  for  Transverse  Vibrations  of 
Elastically  Restrained  Double  Bellows  Expansion  Joint  in  Lateral  Mode 

A  case  of  vibration  of  bellows  in  lateral  mode  is  shown  in  figure  7.1.  It  is  seen 

that  the  pipe  has  pure  transnational  motion  due  to  the  geometry  and  physical 

symmetry  of  the  system  provided.  The  coriolis  component  of  force  acting  on  the 

bellows  from  the  fluid  flowing  inside  is  neglected. 

As  mathematical  approximation,  one  half  of  the  system  is  considered  with  its 
left  end  having  a  rotational  stiffness  ‘R’  and  right  end  fixed  to  the  vertical  rollers. 
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The  boundary  conditions  for  such  a  system  is  given  by- 


i  :i 


I< 


^cCo) 


Deflection  is  zero  at  w  (0,  t)  =  0  at  x=0 


a2x(o>  _  r  ax(O) 

Q  x 2  El  3x 

If  T=RL/EI,  we  get 


(7.14) 

(7.15) 

(7.16) 


a2x(o>  _T  ax(O) 

d  x2  <9x 

(7.17) 


As  the  end  B  does  not  rotate  the  third  boundary  condition,  slope  is  zero  and  given  by 


<3w(L,t) 

<9x 


(7.18) 


The  fourth  boundary  condition  at  end  ‘B’  is  the  shear  force  Q  (L,  t)  of  the  bellows  and 
is  given  by  - 


d3W(L,t ) 

d? 


R+K+»V3Md2W(E0  |  ksW(L,t) 
El  dt2 


El 


El 


(7.19) 
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d3W  (L)  —W2{m  +(m  +  mf3)a]  k 

- W  = - (  s  \  p - OEl*X(L)  +  ^(L) 

dx'  El  El 


- W 2  {Ms  +  (mp  +  m /3  )a}  ks  X (L) 
El  El 


(7.20) 

(7.21) 


Ms  is  the  equivalent  lateral  support  mass,  ks  is  the  equivalent  spring  stiffness  of  the 
lateral  support,  m  p  is  the  mass  per  unit  length  of  the  connecting  pipe  of  length  ‘a’  and 
m  f3  is  the  mass  of  the  fluid  per  unit  length  in  the  connecting  pipe. 


Substitution  of  the  general  solution  (7.8)  and  its  derivatives  into  boundary  condition 
expressions  (7.14),  (7.15),  (7.16)  and  (7.17)  will  give  a  set  of  linear  equations  with 
respect  to  the  constants  A,  B,  C  and  D. 


B  +  D  =  0  (7.22) 

{B(aL)2  - D(pL)2  - T( A  aL  +  C  pL)}=  0  (7.23) 

Aa  cosh  aL+  Ba  sinh  aL  +  CP  cosPL  -  Dp  sinPL  =  0  (7.24) 

[a3  cosh  aL  +  b  sinh  aL]  A  +  [a3  sinhaL  +  b  cosh  aL]  B  -  [p3  cosPL  -  b  sin  pL]  C  + 
[p3  sinpL  +  b  cos  pL]  D  =  0  (7.25) 


Substituting  ci,  C2,  C3,  C4  for  the  terms  in  the  brackets  in  equation  (7.23)  we  get  - 
ci  =  a3  cosh  aL  +  b  sinh  aL  (7.26) 

C2  =  a3  sinh  aL  +  b  cosh  aL  (7.27) 


C3  =  p3  cospL  -  b  sinpL 


(7.28) 
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C4  =  p3  sinPL  +  b  cos  pL 


(7.29) 


For  a  non-trival  solution  the  determinant  formed  by  the  coefficients  of  the  system  of 
algebraic  equations  above  must  be  equal  to  zero  - 

The  expansion  of  the  above  determinant  results  in  the  frequency  equation  of  double 
bellows  in  lateral  mode- 


0 

1 

0 

1 

A 

+  T(a) 

(a)2 

+  T(P) 

+  (P)2 

B 

a  cosh  (aL) 

asinh(aL) 

P  cospL 

-psinpL 

C 

Ci 

c2 

C3 

c4 

D 

(7.30) 


The  frequency  equation  is 

T  {b  [(a2  -  p2)  sinPL  sinhaL  +  2a,p  (1-  cosPL  coshaL)}  -ap  (a2  +  p2)  (a  cosPL 
sinhaL  +p  sinpL  coshaL)}  +  b  (a2  +  p2)  (a  sinPL  coshaL  -  p  cospL  sinhaL)}  - 
ap(a2  +  p2)  cosPL  coshpL  =  0  (7.31) 


Table  7.1:  Exact  Natural  Frequency  for  T=oc,  N=l,2  and  Ms=0 


Spring  Support 
Stiffness 
ks 

Frequency 

N=l,  rad/s 

Frequency 

N=2,  rad/s 

102 

331.4023 

3500.686 

103 

333.41187 

3500.698 

104 

352.8724 

3500.812 

105 

507.8993 

3501.972 

106 

1256.9545 

3515.098 

107 

3204.5183 

4195.031 

108 

3392.1072 

6856.604 

109 

3399.4055 

13673.99 

1010 

3400.1597 

13678.91 
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Table  7.2:  Exact  Natural  Frequency  for  N=l,  2  and  T=oc  &  ks=0 


Support  mass 

Ms 

Frequency,  rad/s 

N=1 

Frequency,  rad/s 
N=2 

1 

210.103 

3439.571 

10 

83.1873 

3406.248 

20 

59.7721 

3403.299 

30 

49.0707 

3402.274 

40 

42.6135 

3401.754 

50 

38.8901 

3401.144 

60 

34.8901 

3401.129 

70 

32.3276 

3401.074 

80 

30.2577 

3400.963 

90 

28.5405 

3400.875 

100 

27.0859 

3400.804 

Fig  7.2  Exact  Transverse  Vibrations  of  Double  Bellows 


I" 


Support  Mass,  Ms 


120 
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Fig  7.3  Exact  Transverse  Vibrations  (Lateral)  of  Double  Bellows 
For  N=  1  &  2  and  MS=Q 


7.4  Mode  Shape  Expression 

The  boundary  conditions  of  double  bellows  subjected  to  lateral  mode  of  vibration  are 
as  follows- 


AtX(x)  =  0, 

<32x  rl  ax 

d  El  <3x 

a2x_  sx 

d  x2  5x 

where 

RE  _ 

- =  Ti 

El 


140 


dx(L) 


=  Oatx 


=  L 


dx 

A.cx  cosh  otL  +  B cx  sinh  otL  +  CP  cos  PL 


Opsin  PL  =  O 


(7.32) 

And  Substituting,  in  equation  (7.32)  we  get- 
B+D=0  ,\B=  -D 

(7.33) 

L2  (B  a2  -  D  R2)  =  ^  (aLA  +  PLC) 

1  nl 

B  (aL)2  -  D  (PL)2  -  Ti(AaL  +  CpL)  =  0 
Substituting  in  equation  (7.33),  we  get- 


i.e. 

0^3(coshaL  +  bsinh  aL)A  +  (Q^  sinh  aL  +  bcoshaL)B  -  (J3  cos  pL  -  b  sin  pL)C 

3 

+  (p  sin  pL  +  bcos  pL)D  =  0,  where 

3 

a  =  cosh  aL  +  b  sinh  aL 

3 

a  =  0£  sinh  aL  +  bcosh  aL 
C3  =  |3  cospL  +  bsinpL 

3 

C4  =  p  sin  pL  + bcos  pL 
-  D(a  L)2  -  D(p  L)2  -  T,  (AaL  +  CpL)  =  0 


-  D[(a  L)2  +  (p  L)2  j-  T,  (AaL  +  CpL)  =  0 

-  T,  (AaL  +  CPL)  =  D[(a  L)2  +  (p  L)2  j 


(7.34) 
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Now, 

Aa  cosh  aL  +  Ba  sinh  aL  +  Cpcos(3L-Dpsin  PL  =  0 
By  substituting  B=  -D  the  above  equation  can  be  written  as- 

Aa  cosh  aL  -  Da  sinh  aL  +  CpcospL-Dpsin  PL  =  0 
Aa  cosh  aL  +  CPcospL-D[a  sinh  aL  +  P  sin  PL]  =  0 

AacoshaL  +  CPcospL  =  D[asinh  aL  +  Psin  PL]  (7.36) 

Multiplying  equation  (7.34)  by  cos  haL  and  equation  (7.36)  by  T iL  and  canceling  the 
like  terms  A  and  C,  we  get  the  ratios  C/D  and  A/D  respectively- 

cosh  aL  x  -  T,  (AaL  +  CpL)  =  D[(aL)2  +  (pL)2  j 

TJL  x  Aa  cosh  aL  +  CP  cos  PL  =  D(a  sinh  aL  +  Psin  PL) 

-T,  AaL  cosh  aL  +  CpL  cosh  aL  =  D[(aL)2  +  (j3L)2JcoshaL 
— Tj  AaL  cosh  aL  +  CpL  cosh  aL  =  D[asinh  aL  +  Psinh  PLJTjL 


CpL[cosh  aL  +  T)  cos  PL]  =  D|(aL)~  cosh  aL  +  (PL)'  cosh  aLj+  [aL  sinh  aL  +  pLsinh  PL]t,  | 

(7.37) 


(aL)’  cosh  aL  +  (PL)’  cosh  aL  -  [a  sinh  aL  +  P  sinh  PlJTj  1 


pL(cosh  aL  +  T,  cos  PL) 
Now  finding  out  the  ratio  A/D- 


-  r 

> 


cos  PL  X  T,  (AaL  +  CpL)  =D[(aL)2  +  (pL)2J 
L  x  Aa  cosh  aL  +  CPcospL  =  D[a  sinh  aL  +  P  sin  PL] 


-Tj  aLAcospL  +  CpLcospL  =  D[(aL)2  +(pL)2J  cospL 
AaL  cosh  aL  +  CpL  cos  PL  =  D[asinh  aL  +  Psinh  PL]L 
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[ijaL  cos  PL  +  aL  cosh  aL]A  =  D|(aL)  cosPL  +  (pL)?  cospLj-[aLsinh  aL  +  PLsin  PL]} 

(7.38) 


A'=  —  =  • 

D 


(aL)2  cospL  +  (pL)2cospL  -  asinh  aL  +  Psinh  PL] 


(T,aL  cos  pL  +  aL  cosh  aL) 


«=-l 

D 

Substituting 

y  =  D 

J  max 


^  a  cosh  aL£,  +  ~  a  sinh  aL£,  +  ~  P  cos  PL£,  -  Psin  PL^ 


=  0 


D 


=  A'acoshaL(£,)-asinh  aL£,  +  C'PcospL^-psin  PL2, 


Where 


A'=  —  =  • 

D 


(aL)2  cos  PL  +  (pL)2  cos  pL  -  [aL  sinh  aL  +  PL  sin  PL] 


(T,aLcos  pL  +  aLcosh  aL) 


c-£- 

D 


(aL)’  cosh  aL  +  (pL)1  cosh  aL  4-  T,  [aL  sinh  aL  +  pL  sin  pL] 


pL(cosh  aL  +  T,  cos  pL) 
Therefore,  the  final  mode  shape  equation  can  be  written  as  - 


^  =  A'  a  cosh  aL^  —  asinh  aL^  +  C'PcospU;  —  Psin  PL£, 

(7.39) 


Where  A'  and  C'  have  already  been  defined  above. 


7.5  Results  and  Discussion 

A  double  bellows  expansion  joint  having  the  following  geometrical 
dimensions  is  taken  and  is  similar  to  the  bellows  considered  by  Jakubauskas.V.F. 
Bellow  length  L  =  0.0693m,  mass  moment  of  inertia  per  unit  length  J  = 
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0.001 153kgm,  El  =  5.078  Nm2,  total  bellows  mass  mtot  =  5.13kg/m,  total  connecting 
pipe  mass,  mp+mf3  =  5.0kg/m  and  assuming  a  =  L. 


According  to  [1],  the  maximum  allowable  pressure  in  bellow  is  - 

_  7lkP 

6.666L2 

Where  k  is  the  equivalent  axial  stiffness  of  bellows  - 


(7.32) 


Fig  7.4  Exact  Transverse  Vibrations  of  Double  Bellows 
For  N=l,  ks=0 


Support  Mass,  Ms 


Substitution  of  the  above  numerical  values  into  the  expressions  7.4  &  7.5,  we  get  - 
c=^61 6.49  +  0.0001 135co2  (7.33) 

A, =1.00297©  (7.34) 

b  =0.10235co2  (7.35) 
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The  study  involves  the  effect  of  varying  the  equivalent  support  stiffness  and 
mass  on  the  frequency  of  vibration.  Using  expressions  in  (7.12),  (7.13),  (7.32),  (7.33), 
(7.34)  and  (7.35)  the  frequency  equation  (7.31)  is  solved  by  applying  bisection 
method.  The  results  are  presented  in  Tables  7.1  and  7.2  respectively. 

Table  7.1  presents  the  frequencies  obtained  by  varying  the  support  stiffness  ks 
for  modes  of  vibration  N=1  and  2.  It  is  seen  that  as  stiffness  increases  from  a  value  of 
102  to  1015,  frequency  increases  from  331.40rad/s  to  3400.1597rad/s,  a  steep  rise  by 
almost  90%.  The  frequency  attains  almost  a  constant  value  from  108  onwards. 

The  effect  of  varying  the  support  mass  on  the  frequency  of  vibration  is  also 
studied  and  the  results  presented  in  Table  7.2.  It  is  seen  that  there  is  a  steep  drop  of 
87%  in  frequency  as  the  support  mass  increases  from  1.0  to  100  N=1  and  by  a  small 
drop  of  1%  for  N=2. 


Table  7.3  Exact  Lateral  Mode  Frequency  for  T=oc  and  Ms  =0 


S.  No 

ks,  spring 

stiffness 

Frequency,  co  (rad/s)  N=1 

Frequency,  ©  (rad/s)  N=2 

1 

102 

331.4 

3500.68 

2 

1010 

3400.1597 

13678.906 

Table  7.4  Exact  Lateral  Mode  Frequency  for  T=oc  and  ks  =0 


S.  No 

Ms,  support 

mass 

Frequency,  ©  (rad/s)  N=1 

Frequency,  ©  (rad/s)  N=2 

1 

1 

210.0 

3439.5 

2 

100 

27.0 

3400.8 

Fig  7.5  First  Mode  Shape  of  Double  Bellows  with  Elastically  Restrained  Ends 

2 

for  Rotational  RestraintT=infinity,  Support  stiffness  ks  at  10  &  infinity, 
Support  mass,  Ms=0  and  N=1 


ks=1 00 
■  ks=infinity 


zeta=x/L 


Fig  7.6  Second  Mode  Shape  of  Double  Bellows  with  Elastically 
Restrained  Ends  at  T=infinity,  Support  stiffness  ks  at  102  &  infinity, 
Support  mass  Ms=0  and 

N=2  ^ 


0.2  0.3 

zeta=x/L 


■  ks=1 00 
ks=infinity 


Table  7.5:  y/D  at  T=oc,  Ms=0  and  N=1 


% 

y/D,  ks=100 

y/D,  ks=oc 

0.1 

0.9x1 0"8 

-0.1 4x1 0"6 

0.2 

0.1 8x1 0"7 

-0.27x1 0"6 

0.3 

0.4x1 0"7 

-0.39x1 0"6 

0.4 

0.37x1 0"7 

-0.55x1 06 

0.5 

0.66x1 07 

-0.72x1 0"6 

Table  7.6:  y/D  at  T=oc,  Ms= 

0  and  N=2 

5 

y/D,  ks=100 

y/D,  ks=oc 

0.1 

0.83x1 0"8 

0.1 7x1 0"5 

0.2 

0.1 6x1 0"7 

0.34x1 0"5 

0.3 

0.23x1 0"7 

0.51x1 0"5 

0.4 

0.33x1 07 

0.68x1 0"5 

0.5 

0.43x1 0"7 

0.85x1 05 
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Figures  7.5  and  7.6  indicate  the  mode  shapes  of  double  bellows  with 
elastically  restrained  ends.  The  support  stiffness  ks  is  varied  and  its  effect  on  the 
mode  shape  is  studied  at  ks=100  and  ks=oc,  and  by  keeping  the  support  mass,  Ms=0 
and  T  =oc.  The  values  y/D  are  obtained  for  two  modes  of  vibration,  N=1  &  2 
respectively.  The  ordinate  in  Fig  7.5  &  7.6  is  the  ratio  of  minimum  to  maximum  y/D 
plotted  for  different  values  of  £  ranging  from  0  to  0.5  for  the  first  half  mode.  It  is 
found  that  as  £,  varies  between  0  to  1  and  support  stiffness  ks  increases,  the  mode 
shape  follows  an  up  trend  with  a  fall  and  rise  at  £=0.2  and  0.3  and  remains  a  straight 
line  approaching  infinity.  This  phenomenon  is  observed  for  both  the  modes  of 
vibration  N=1  and  2  respectively. 


Fig  7.7  Mode  Shape  of  Double  Bellows  with  Elastically  Restrained  Ends 
for  T=T,=T2  and  varying  from  0.01  to  100,  Support  stiffness  ks=100  and 

Support  mass  Ms=0 
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Figure  7.7,  shows  the  mode  shape  of  double  bellows  by  varying  the  rotational 
restraint  parameter  T,  from  0.01  to  100.  The  support  stiffness  ks  is  kept  constant  at 
100  and  support  mass,  Ms  =0.  The  plot  between  £  and  ratio  y/D-  indicate  that  as  T 
increases,  the  ratio  y/D  decreases.  It  is  observed  that  there  is  a  sharp  rise  in  y/D  at 
T=100,  and  £=0. 1 .  However,  no  significant  change  is  observed  in  the  ratio  y/D  by 
varying  T  from  0.01  to  10. 


Table  7.7:  y/D  for  T=0.01  to  100 


5 

y/D,T=0.01 

y/D,T=0.1 

y/D,T=1.0 

y/D,T=10 

y/D,T=100 

0 

0 

0 

0 

0 

0 

0.1 

-4.73 

-4.73 

-4.98 

-7.74 

-1.04 

0.2 

-9.42 

-9.47 

-9.97 

-15.5 

-2.08 

0.3 

-14.1 

-14.9 

-14.9 

-21.2 

-3.13 

0.4 

-18.8 

-18.9 

-19.9 

-30.9 

-4.17 

0.5 

-23.5 

-23.6 

-24.9 

-39.74 

-5.21 

0.6 

-18.8 

-18.9 

-19.9 

-30.9 

-4.17 

0.7 

-14.1 

-14.9 

-14.9 

-21.2 

-3.13 

0.8 

-9.42 

-9.47 

-9.97 

-15.5 

-2.08 

0.9 

-4.73 

-4.73 

-4.98 

-7.74 

1 

0 

0 

0 

0 

0 
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CHAPTER  8 

TRANSVERSE  VIBRATIONS  OF  UNIVERSAL  EXPANSION 
JOINT  RESTRAINED  AGAINST  ROTATION 
(ROCKING  MODE) 

8.1  Derivation  of  Frequency  Equation  for  Transverse  Vibrations  of  Elastically 
Restrained  Double  Bellows  Expansion  Joint  in  Rocking  Mode 

In  case  of  the  vibration  of  bellows  in  rocking  mode  the  middle  point  of  the 

connecting  pipe  does  not  translate  because  of  the  geometrical  and  physical  symmetry 

of  the  system  with  respect  to  the  imaginary  vertical  axis.  The  coriolis  forces  acting  on 

the  bellows  from  the  fluid  inside  are  neglected.  Therefore,  as  a  mathematical 

approximation,  one  half  of  the  physical  system  can  be  taken  with  its  left  end  fixed  and 

right  end  simply  supported  in  the  middle  of  the  connecting  pipe  as  shown  in  fig  8.3 

Therefore,  one  half  of  the  system  can  be  taken  with  its  left  end  ‘A’  fixed  and 
having  a  rotational  stiffness  ‘R’  while  the  right  end  ‘B’  is  simply  supported  and  held 
at  the  middle  of  the  connecting  pipe. 

8.2  Boundary  Conditions  for  Transverse  Vibrations  of  DBEJ  Rocking  Mode 

As  shown  from  figure  8. 1  the  boundary  conditions  at  the  end  ‘A  &  B’  are  given  by  - 

dX(L)  _  X(L) 
dx  a 
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Fig  8.1  Mathematical  Model  of  Double  Bellows  Expansion  Joint  (Rocking  Mode) 


d3X(L) 

dx3 


d2X(L)  L 
dx3  a 


PUR,;  dX(L) 
El  dx 


J  CQ2  dX(L) 
El  dx 


2  Jtot  ks 

CO  ~ — 2X(L)  +  —  X(L) 


El 


a 


EI 


w  (0,  t)  =  0  at  x  =  0 

(8.1) 

B  +  D  =  0  and 

(8.2) 

EI  (pa2  -  Dp2)  =  R  (Aa  -  Cp) 

(8.3) 

dw  (L,  t)  1 

=  w  (L,  t) 
ax  a 

(8.4) 

a[Aa  cosh  aL  +  Ba  sinh  aL  +  CPcos  pL  -  DPsinPL]  +  [A  sinh  aL  +  B  cosh  aL  +  C 
sinpL  +  D  cosPL]  =  0  (8.5) 
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(a  a  cosh  aL  +  sinh  aL)  A  +  (aa  sinhaL  +  cosh  aL)  B  +  (Pa  cosPL  +  sinpL)  C  -  (pa 
sinpL  -  cosPL)  D  =  0  (8.6) 

Where 


di  =  (aa  cosh  aL  +  sinh  aL)  (8.7) 

d2  =  (aa  sinhaL  +  cosh  aL)  (8.8) 

d3  =  (Pa  cosPL  +  sinpL)  (8.9) 

d4  =  (Pa  sinpL  -  cospL)  (8.10) 


Another  boundary  condition  at  end  ‘B’  was  derived  from  the  rotation  differential 
equation  for  one-half  of  connecting  pipe  OB  - 

g3W  (L,  t)  _  a3W  (L,  t)  1  Pnn  J  dW  (1,  t)  |  kbW  (1,  t)  | 

<9x3  dx2  a  El  dx  El 

J*d3ft>(L,t)  |  JTotd20)(L,t) 

El *dxdt2  El* a2  St2 

The  total  mass  moment  of  inertia  of  fluid  and  metal  of  one-half  of  the  connecting  pipe 
about  point  of  rotation  ‘O’,  including  the  inertial  of  lateral  supports  is  expressed 
analytically  as  follows  [1]  - 


J  =mp+mQ  *a3  + 

J  tot  ~  Cl  -r 


2mp+ms 


*aR2  +M  a2 


(8.12) 


As  in  the  lateral  mode  of  analysis  taking  for  ‘w’  in  the  expression  (8.2),  the  boundary 
conditions  (8.2),  (8.3),  (8.4)  and  (8.5)  can  be  simplified  by  eliminating  the  time 
harmonic  function  and  rewriting  the  expression  as  - 
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X  (0)  =  0 


(8.13) 


d^C(U=_l£X(U_2cIdX(L)b 
dx3  a  dx2  dx 


(8.14) 


Where  ‘c’  and  ‘A,’  are  defined  in  equations  (7.5)  &  (7.6)  respectively  and  b  is  given 

by  b  =W2  (8.15) 

Ela2  El 

Using  the  equation  (7.5)  expression  (8.14)  is  written  as  - 

{a3  cosh  aL  +  (a2/a)  sinh  aL  +  2c2  coshaL  +  b  sinh  aL}A  +  {a3  sinh  aL  +  (a2/a) 
cosh  aL  +  2c2  a  sinhaL  +  b  cosh  aL}  B  -  {p3  cospL  +  (p2/a)  sinPL  -  2c2p  cospL  -  b 
sinpL}C+  {p3  sinpL  +  (p2/a)  cospL  -  2c2psinpL  +  b  cos  PL}D  =  0  (8.16) 

Grouping  all  the  like  terms  together  we  get  - 

{a(a2  +  c2)  cos  haL  +  (a2/a+b)  sin  haL}  A  +  {a(a2  +  2c2)  sin  haL  +  (a2/a  +b)  cosh 


pL}  B  -  { P(p2  -  2c2)  cos  pL  +(p2/a-b)  sin  pL}  C  +  { p(p2 

pL}  D  =0 

-  2c2)  sinpL  -  (P2/a  -  b) 

(8.17) 

If 

ai  =  a2/a  +  b 

(8.18) 

a2  =  a(a2  +  2c2) 

(8.19) 

bi  =  p2/a  -  b 

(8.20) 

b2  =  p(p2  -  2c2) 

(8.21) 

ci  =  ai  sinh  aL  +  a2  cosh  aL 

(8.22) 

c2  =  ai  cosh  aL  +  a2  sinh  aL 

(8.23) 
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C4  =  b2  sin  pL  -  bi  cos  pL  (8.24) 

Then  equation  (8.17)  can  be  written  as  follows  - 

ci  A  +  C2B  -  C3C  +  C4D  =  0  (8.25) 


The  determinant  formed  by  the  coefficients  of  this  system  of  algebraic  equations  must 
be  set  equal  to  zero  - 


0  10  1 

A 

+  T1(aL)  (aL)2  +T1(pL)  -(pL)2 

B 

d,  d2  d3  d4 

C 

Ci  c2  c3  c4 

D 

(8.26) 


Expanding  the  above  matrix,  we  get  the  final  frequency  equation  of  double  bellows  in 
rocking  mode- 

[-aPabi  -  pa2  +  aPaai  +  ofe}  (1  -  cosh  aL  cos  pL)  -  [-P2aai  -  pb2  +  cx,2abi  +  aa2]  (sin 
PL  sinh  aL)  +  [aPab2  +  ap2a2  +  aai  +  abi}  cosh  aL  sinPL  +  [Pai  +  pbi  -  aa2b2 
+aPa2a]  sinh  aL  cos  pL  =  0  (8.27) 


8.3  Derivation  of  Mode  Shape  Expression 

At  X(0)  =  0 


Hax(0) 

<9x~ 


dx(0) 

R. -  at  the  end  A 

1  dx 


dx(L) 


1 

—  x 


(L) 


dx 


a 
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-S-W(L,t)  +  — 


j  a2w(L,T)  ,  jwt  .  a2w(L,t) 


El 


El  5x2. at2  El.a- 


at 


The  above  equation  can  be  written  as- 

d3(L)  =  1  d2x  (L)  2c2dx(L 

dx3  a  dx2  dx 


-  bx  (L) 


Where 


PTtRm7  +  J  00 
El 


=  c 


At  X  (0)  =0 

X  =  A  sinhaL  +  B  coshaL  +  C  sinpL  +  D  cosPL 
Applying  X  (0)  =0  in  above  equation,  we  get- 
B  +D  =  0 
B  =  -  D 


d  x  (0) 

Now  - —  =  0 

dx 


(8.28) 

(8.29) 

(8.30) 


dX 

- =  Aa  CoshaL  +  Bp  SinhaL  +  Cp  Cos  pL  -  Dp  SinpL 

dx 


(8.31) 


d  X 
dx2" 


Aa2  CoshaL  +  Bp2  SinhaL  +  Cp2  Cos  pL  -  Dp2  SinpL 


(8.32) 


El 


d  2X  (0) 
dx2 


=  Ba2  -  Dp2 


(8.33) 


R[  =  (Aa+Cp)  R, 

dx 


(8.34) 
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Equatmg  the  above  two  equations,  we  get- 


Bot2-DB2  =  —  (Aa  +  CB) 

El 

Multiplying  by  ‘L’  on  both  sides 

B(aL)2  -  D(PL)2  =  ^t(Aa+CP) 

El 

B(aL)2  -  D(pL)2  =  T,  (Aa  +  Cp) 

T,  =  Substituting  B  =  -  D 

-D(aL)2-  D(pL)2  =  Tj  (Aa  +  Cp) 
-  D[(aL)2  -  D(pL)2  ]  -  T,  (Aa  +  Cp) 
We  know. 


(8.35) 

(8.36) 

(8.37) 

(8.38) 

(8.39) 

=  0 


d  X  (L) 
dx 


And, 

a 


=  Aa  coshaL  +  Ba  sinhaL  +  CP  cos  PL  -  Dp  sinpL  (8.40) 


x  (L)  =  —  [A  sinh  aL  +  B  cosh  aL  +  C  sin  pL  +  D  cos  PL] 

a 


(8.41) 

.'.  a[Aa  Cosh  aL  +  Ba  SinhaL  +  Cp  Cosh  aL  +  D  Sin  pL]  +  [  A  Sinh  aL  + 
B  Cosh  aL  +  C  Sin  pL  +  D  CospL]  =0  (8.42) 


A[aacCosh  aL  +  sinh  aL]  +  B[aa  sinh  aL  +  cosh  aL]  +  C[aP  cospL  +  sin  pL]  - 
D[aP  sin  pL  -  cos  pL)]  =  0  (8.43) 
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Writing  B  =  -  D, 

A  [aa  cosh  aL  +  sinh  aL]  +  D  [aa  sinh  aL  +  cosh  aL]  +  C  [ap  cosPL  +  sin  pL]- 
D  [ap  sin  pL  -  cos  pL)]  =  0  (8.44) 

=  A  [aa  cosh  aL  +  sinh  aL]  +  C  [ap  cos  pL  +  sin  pL]  +  D  [aa  sin  aL  +  coshaL  + 
ap  sin  pL  -  cos  pL]  =  0  (8.45) 

We  know,  -Ti  (Aa  +  CP)  =  D  (aL)2  +  (pL)2] 

A  [aa  cosh  aL  +  sinh  aL]  +  C  [ap  cos  PL  +  sin  pL]  =  D[aa  sin  aL  +  cosh  aL  +  ap 
sin  pL  -  cos  pL]  (8.46) 

Multiplying  LHS  by  (aa  cosh  aL  +  sinh  aL)  and  RHS  by  TiAa  we  get- 

-  TiA  (aa  cosh  aL  +  sinh  aL)  -  TiB  (aa  cosh  aL  +  sinh  aL)  =  D  (aa  cosh  aL  + 
sinh  aL)  [(aL)2  +  (pL)2] 

TiAa  (aa  cosh  aL  +  sinh  aL)  +  C  Tia  (ap  cos  PL  +  sin  pL)  =  DTia  (aa  sinh  aL  + 
cosh  aL  +  ap  sinpL  -  cos  PL] 


Canceling  ‘A’  terms  and  adding,  we  get- 

-  TiCP(aa  cosh  aL  +  sinh  aL)  +  C  Tia  (ap  cos  pL  +  sin  pL)  =  D  [aa  cosh  aL  +  sinh 
aL)  [(aL)2  +  (PL)2]  -  T ia(aa  sinh  aL  +  cosh  aL  +  ap  sinpL  -cos  PL]  (8.47) 
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C{-TiP(aa  cosh  aL  +  sinh  aL)  +  Tia  (ap  cos  pL  +  sin  pL)  =  D[aa  cosh  aL  +  sinh 
aL)  [(aL)2  +  (PL)2]  -  T  ia(aa  sinh  aL  +  cosh  aL  +  aP  sinPL  -  cos  pL] } 

(8.48) 


C  I  TjP(aa  Cosh  aL  +  Sinh  aL)  +  Tja  (aP  Cos  PL  +  Sin  PL) 

D  (aaCoshaL  +  Sinh  aL)[(aL)2  +  (PL)2  -TjafaaSinh  aL  +  Cosh  aL  +  aP  Sin  PL  -  Cos  PL] 


(8.51) 

Now  canceling  ‘C’  terms  to  obtain  the  ratio  A/D- 

TiAa(ap  cos  pL  +  sin  pL)  +  Cp  (ap  cos  pL  +  sin  PL)  =  D  [aa  cos  pL  +  sin  pL)  (aa 
sinh  a  +  cosh  aL  +  ap  sinPL  -  cos  PL] 

AP(aa  cosh  PL  +  sinh  aL)  +  Cp  (ap  cos  pL  +  sin  pL)  =  Dp  [aa  sinh  aL  +  cosh  aL  + 
ap  sinpL  -  cos  pL] 

Canceling  ‘C’  terms  we  get, 

TiAa(ap  cos  PL  +  sin  pL)  +  Ap  (ap  cosh  aL  -  sinh  aL)  =  D  [ap  cos  PL  +  sin  pL) 
(aa  sinh  a  +  cosh  aL  +  ap  sinpL  -  cos  PL]  +  p[aa  sinh  aL  +  cosh  aL  +  ap  sinpL  - 
cos  PL]  (8.52) 

A  [-Tia(aP  Cos  pL  +  Sin  pL)  +  P  (aa  Cosh  aL  -  Sinh  aL)  =  -D  [aPa2  sinh  a  cos 
pL  +  ap  cos  pL  cosh  aL  +  a2  p2  cos  pL  sin  PL  -  ap  cos2  pL  +  aa  sin  PL  sinh  aL  + 
cosh  aL  sin  pL  +  ap  sin2  PL  -  cos  pL  sin  pL]  +  P  [(aa  sinh  aL  +  cosh  aL  ap  sin  pL 


-  cos  PL)] 


(8.53) 
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=  D{[aPa2  sinh  a  cos  PL  +  ap  cos  pL  cosh  aL  +  cos  pL  sin  pL(a2-l)  -  aP  (cos2 
pL  +  sin2  pL)  +  aa  sin  pL  sinh  aL  +  cosh  aL  sin  pL]  +  p  (aa  sinh  aL  +  cosh  aL  + 
ap  sin  pL  -  cos  pL] } 


A  _  I  T,a  (ap  Cosh  PL  +  Sin  PL)  +  P(aa  Cosh  aL  -  Sinh  aL) 

D  (aPa2Sinh  aL  CoshpL  +  apCosPL  CoshaL  +  Cos  PLSin  PL(a2  -1) 

-ap  +  aa  Sin  P  Sinh  aL  +  CoshaLSin  PL]  +  P[aa  Sinh  aL+  CoshaL  + 
apSin  PL  -  CosPL] 


D 


B 

D 


=  -1 


The  final  mode  shape  expression  can  be  written  as- 

y(£) 

- =  A'sinh  aL  -  B cosh  aL  +  C sin  BL  +  DcosBL 

D 

8.4  Results  and  Discussion 

The  above  frequency  solutions  for  rocking  mode  are  tested  by  considering  a 
bellow  with  the  following  geometrical  and  physical  parameters  - 
Bellows  length,  L  =  0.00693  m,  mass  moment  of  inertia  per  unit  length,  J  = 
0.001 153kgm,  El  =  5.078Nm2,  total  bellows  mass,  mtot  =  5.138kg/m  and  total 
connecting  pipe  mass,  m  p  +  me  =  5.0kg/m. 
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The  study  involves  the  effect  of  varying  the  equivalent  support  stiffness  and 
mass  on  the  frequency  of  vibration.  Using  expression  (7.11),  (7.12)  and  (8.12),  (8.13), 
(8.14),  (8.15)  the  frequency  equation  (8.27)  is  solved  using  bisection  method.  The 
first  fundamental  natural  frequencies  obtained  for  rocking  modes  are  given  in  Tables 
8.1  &  8.2  by  varying  spring  stiffness  k  s  and  support  mass  Ms  respectively. 

It  is  seen  that  as  stiffness  increases  from  a  value  of  102  to  1015,  frequency 
increases  from  936.10rad/s  to  3856.1719rad/s,  a  rise  by  almost  75%.  The  frequency 
attains  almost  a  constant  value  form  108  onwards.  The  effect  of  varying  the  support 
mass  on  the  frequency  of  vibration  is  also  studied  and  the  results  presented  in  Table 
8.2.  It  is  seen  that  there  is  a  steep  drop  of  87%  in  frequency  as  the  support  mass 
increases  form  1.0  to  100  for  N  =  1  and  by  a  small  drop  of  2%  for  N  =  2.  The  exact 
frequencies  obtained  by  the  bisection  method  are  compared  with  the  results  presented 
in  the  thesis  [1]  and  are  presented  in  Table  8.3. 

It  is  seen  that  the  percentage  error  obtained  in  the  frequency  of  vibration  for  a 
bellow  with  a  lateral  support  having  a  stiffness  of  ks  =  102  quite  small  and  comparable 


to  a  bellow  with  no  lateral  supports. 


Frequency,  rad/s 
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Table  8.2  Rocking  Natural  Frequency  for  N=l,  2  and  T=oc  &  ks=0 


Support  Mass 
ms 

Exact,  rad/s 

N=1 

Exact,  rad/s 

N=2 

1 

429.1701 

3480.089 

10 

164.6297 

3411.481 

20 

117.9629 

3405.954 

30 

96.75 

3404.053 

40 

83.977 

3403.092 

50 

75.2143 

3402.511 

60 

68.7234 

3402.122 

70 

63.667 

3401.844 

80 

59.584 

3401.635 

90 

56.197 

3401.472 

100 

53.33 

3401.341 

Fig  8.3  Exact  Rocking  Mode  Natural  Frequency  -Double  Bellows 
For  N=l,2  and  k  =0 


■  N=1 
•N=2 


Support  Mass,  Ms 


Table  8.3  Comparison  of  frequencies  at  T=oc 


Mode  N 

[28]  to,  rad/s 

Exact  [46]  co,  rad/s 

Percent  Error 

ks=102 

o 

II 

P* 

1 

935.686 

936.107 

0.01 

2 

3855.928 

3856.1719 

0.006 
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Figures  8.4a, b&c  represent  the  first  three  mode  shapes  of  bellows  in  rocking 
mode. 


w 

t 


a)  First  mode  shape 


w 


b) Second  mode  shape 


c)  Third  mode  shape 


Fig  8.5  First  Mode  Shape  of  Double  Bellows  with 
Elastically  Restrained  end  -(Rocking  Mode)  Support 
Mass,  Ms=0  &  N=1 ,  T=infinity 


4.E+09 
2.E+09 
0.E+00 
-2.E+09  0 


ks=100 

ks=infinity 


0.5  1 

zeta=x/L 
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Fig  8.6  Second  Mode  Shape  of  Double  Bellows  with 
Elastically  Restrained  End,  Support  Mass  Ms=0,  T=infinity  & 

N=2 


The  first  and  second  mode  shapes  for  the  double  bellows  in  rocking  mode  are 
shown  in  Figures  8.5and  8.6  respectively.  By  varying  c  between  0  to  1.0  for  the  full 
cycle  and  keeping  rotational  restraint  T  =cc,  the  response  is  observed  for  support 
stiffness  ks=100  and  ks=oc.  It  is  seen  that  as  the  support  stiffness  ks,  increases  the 
response  of  bellows  changes  signs.  This  should  be  a  good  indicator  of  controlling  the 
bellows  vibration  by  altering  the  stiffness. 
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CHAPTER  9 

SEISMIC  RESPONSE  OF  ELASTICALLY  RESTRAINED 
SINGLE  BELLOWS  EXPANSION  JOINT  IN  LATERAL  MODE 

9.1  Background 

It  is  seen  that  most  of  the  nuclear  and  thermal  power  plants  around  the  world 
are  located  in  the  seismic  zones.  The  piping  in  these  plants  contains  expansion 
bellows  and  needs  to  be  analyzed  from  seismic  design  point  of  view.  Hence,  while 
designing  an  expansion  bellow  from  structural  point  of  view,  the  seismic  design 
aspect  is  also  to  be  given  utmost  importance. 

The  present  work  therefore  attempts  to  derive  an  exact  solution  for  the  seismic 
response  of  U  type  of  single  bellows  that  are  considered  elastically  restrained  against 
rotation  to  classical  fixed-fixed  case  considered  by  Morishita  et  al.  Hence,  it  is 
observed  that  this  assumption  will  lead  to  an  over  estimation  of  natural  frequencies 
and  incorrect  determination  of  seismic  response. 

The  fundamental  equation  of  motion  of  the  equivalent  Euler-Bernoulli  beam  is 
used  to  calculate  the  spectral  response  of  the  bellows  to  a  seismic  excitation  in  the 
lateral  direction.  The  present  study  excludes  shear  and  inertia  terms  of  the  equivalent 
Timoshenko  beam  for  the  bellows.  The  first  lateral  mode  frequencies  are  obtained  for 
the  case  of  equal  rotational  spring  stiffness  on  either  end. 
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9.2  Simplified  Analysis  of  Lateral  Mode  of  Vibration 

Consider  a  uniform  Euler-Bernoulli  beam  with  its  both  ends  rotationally 
restrained  which  has  the  equivalent  dynamic  characteristics  as  a  bellows  and  is 
excited  in  lateral  direction  as  shown  in  Fig  9.1.  Here  we  use  the  exact  method  to 
obtain  the  natural  frequencies.  Let  the  deflection  of  the  beam  by  bending  be  y  (z,  t) 
and  assuming  sinusoidal  vibration- 


Fig  9.1  Single  Bellows  Restrained  Against  Rotation 

y  (z,  t)  =  Y  (z)  e  ltot  (9.1) 

Where‘co’  is  the  natural  frequency  of  vibration,  t  is  time  in  seconds,  Z=  z/L,  the  non- 
dimensional  length  of  bellows  and  Y  (z)  is  given  by- 

Y  (z)=  A  sinh  az  +  B  cosh  az  +  C  sinaz  +Dcosaz  (9.2) 

Where  A,  B,  C  and  D  are  arbitrary  constants 

Equation  (9.2)  is  considered  as  the  solution  to  the  differential  equation  - 


Where 


mtot  co 

El 


(9.4) 


Where  El  is  bending  stiffness,  R  m  -mean  radius  of  bellows,  mtot  is  the  total  mass  of 


bellows  per  unit  length  that  includes  mass  of  bellows  and  fluid  mass  [40]. 


The  physical  quantities  of  the  beam  such  as  El  is  given  in  terms  of  the  dimensions 
and  material  constants  of  the  bellows  and  are  as  follows- 


5  (dp.tp)qE 
24  h3.Cr 


(9.5) 


pA  =  27rpb 


C 

t.L 

l  P 


(9.6) 


Where  dp  is  the  pitch  diameter  of  bellows,  tp  is  the  wall  thickness  of  bellows,  h  is 


height  of  convolution,  p  <b,  f)  density  of  bellows  material  and  fluid  and  q  is  the  pitch. 
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The  frequency  equation  is  derived  and  the  first  fundamental  transverse  frequency  is 
obtained  using  bisection  method. 


9.4  Boundary  Conditions 


The  boundary  conditions  for  the  rotationally  restrained  Euler-Bernoulli  beam  /  bellow 

are  as  follows- 

i)  y  =0  at  z=0 

(9.7) 

ii)  o2y  /  dz~  =  +  T  i .  dy/dz  at  z=0 

(9.8) 

iii)  And  y  =  0  at  z  =  L 

(9.9) 

iv)  <92y  /  dzr  =  -  T2  Sy/Sz 

(9.10) 

Let  the  solution  for  the  differential  equation  (9.3)  be  written  as- 

y(z)  =  A  Sinh  a  z  +  B  Cosh  a  z  +  C  Sin  a  z  +  D  Cos  a  z 

(9.11) 

—  =  a  (A  Cosh  az  +  B  Sinh  az  +  C  Cos  az  -  D  Sin  az) 
dz 

(9.12) 

j2 

— \  =  a2  (A  Sinh  az  +  B  Cosh  az  -  C  Sin  az  -  D  Cos  az) 
dz2  v  ' 

(9.13) 

(i)  y  =  0  at  z  =  0;  B  +  D  =  0 

(9.14) 

(ii)  f?=+T,  «*;  at  z  =  0 
dz  dz 

(9.15) 

(a  L2  )(B  -  D)  =  +  T,  a(A  +  C) 

(9.16) 

B  +  D  =0  .-.  D  =  -B 
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.-.a2(B  +  B)+  T,  a  (A  +  C) 
2(a  L)(BL)  =  T,  (A  +  C) 

T‘  ;(A+C) 


B  = 


2(aL) 


D  =  -B  =  — 


2(aL) 


(A  +  C) 


(9.17) 

(9.18) 


(iii )  y  =  0  at  z  =  1 

A  SinaL+  B  Cosh  aL  +  C  Sin  aL  +  D  CosaL  (9.19) 

(aL)2  (A  sin  haL+  B  cos  haL  -  C  sin  aL  -  D  cos  aL)  = 

— T2  (aL)(A  Cosh  aL+B  Sinli  aL+  C  Cos  aL-D  SinaL)  (9.20) 

.'.  (aLXASinh  aL  +  B  Cosh  aL-D  Cos  a  L) 

+  T2  (A  Cosh  aL  +  B  Sinh  a  L  +  C  Cos  aL-D  Sin  a  L)  =  0 

(9.21) 

We  know- 

ASinh  aL  +  B  Cosh  aL+  C  Sinh  aL  +  D  Cos  aL  =0 

(a  L( A)  +  T2  (BL))Sinh  a  L  +  (a  LB  +  T2 A)  Cosh  a  L  —  (a  LC  +  T2D)  Sina  L 

-(aLD-T2c)CosaL  =  0  (9.22) 


A  Sinh  aL  +  B  Cosh  aL  +  C  Sin  aL-B  Cos  aL  =  0 


A  Sinh  aL+  C  Sin  aL+  B  (Cosh  aL—  Cos  aL)  =  0 


(9.23) 
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(aLA  +  T,B)  Sinh  aL  +  (aB  +  T,  A)Cosh  aL 

-(aLC  +  T2B)  Sin  aL  -  (-  aLB  +  T2C)  Cos  aL=  0  (9.24) 


(aLA  +  T,B)  Sinh  aL + (aB  +  T, A)Cosh  aL 

-(aLC -T2)  Sin  aL  +  (aLB  +  T2C)  Cos  aL=  0 

(A  Sinh  aL  +  C  Sin  aL) 

(Cosh  aL  -  Cos  aL) 

BIT,  Sinh  aL+  a  Cosh  aL  +  X,  Sin  aL  +  aCosaL  } 

+  A{a  Sinli  aL  +  T2  CoshaL}+C{T2  CosaL-aSinaL}  =  0 


B  =  - 


^  A  Sinh  aL+  C  Sin  aL^ 
v  Cosh  aL- Cos  aL  j 


B  = 


^  n  ^ 


\2aL  j 


(A+C) 


^  n  ^ 

2a 


|(a+c)  =  - 


A  SinhaL  +C  Sin  aL 


v 


(Cosh  aL  -  Cos  aL)  J 

7t(A  +  C)  (  CoshaL  -CosaL)  =  -2a(ASinh  aL  +CSin  aL) 
A{Tj (CoshaL  -  Cos  aL)  +  2aSinhaL] 

=  -{^(CoshaL  -  Cos  aL)  +  2aSinhaL}  C 


C 

A 


Tj (CoshaL  -  CosaL)  +  2a SinhaL 
Tj  (CoshaL  -  CosaL)  +2aSinaL 


Substituting  for  B  in  last  equation  -  we  get 


—  (A  +  C){Tt  SinhaL  +  aL  CoshaL  +  T,  SinaL  +  aL  CosaL} 
2a 


(9.25) 


(9.26) 

(9.27) 

(9.28) 

(9.29) 


(9.30) 


+  A(a  LSinh  aL  +  T2  Cosh  aL )  +  C(T,  Cos  aL  —a SinaL)  =  0 


(9.31) 
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.'.  T,(A  +  C)(aCosh  aL+ T,  sinh  aL+acosaL+T,  sin  aL) 

+  2aA(asinh  aL+ T2  coshaL)  +2ac(T2  cos  aL- a  sin  aL)  =  0  (9.32) 


A{T[(a  coshaL  +T2  sinh  aL+  acosaL  +  T2  sin  aL)+  2a(asinh  aL+T2  coshaL)} 

—  C{Tj(a  coshaL  +  T2  sinh  aL  +  acosaL  +T?  sin  aL)  —  2a  (asin  aL— T2  CosaL)} 

C  Tj(a  coshaL  +T2  sinh  aL+  acosaL  +  T2  sin  aL)+2a(asinh  aL+T2  coshaL) 

A  T,(a  coshaL  +  T2  sinh  aL  +  acosaL  +T2  sin  aL)  -  2a  (asin  aL-T2 CosaL) 


(9.33) 


Equating  the  two  expressions  for  —  we  obtain- 

A 

Tj(cosh  aL-cosaL)  +2asinh  La 
T,  (cosh  aL-cosaL) +  2  a  sin  aL 

T,  (a  cosh  aL  +  T2  sinh  aL + a  cos  aL + T,  sin  aL)  +  2  a  (a  sinh  aL  +  T2  cosh  aL) 
Tj(a  coshaL+  T2  sinh  aL+acosaL+T2sin  aL)- 2  a  (a  sinh  aL-T2  coshaL) 

(9.34) 


By  cross-multiplying  we  get- 
(Tj(cosh  aL-cosaL) +2asin  aL)| 


T,(a  coshaL+  T2  sinh  aL  +  acosaL +T2  sin  aL) 
+  2  a  (a  sinh  aL  +  T2  cosh  aL) 


=  (Tj  (cosh  aL  -  cos  aL)  +  2  a  sin  aL  j 


tTj(a  coshaL+  T2  sinh  aL+acosaL+T,  sin  aL) 
-2a(asin  aL-T2  cosaL) 


(9.35) 


f 

if  (cosh  aL-cosaL)  + 

V 


^a  cosh  aL  +  T2  sinh  aL  +  a  cos  aL  +  T2  sin  aL  +Tj  a  (cosh  aL 
v(asinh  aL  +  T2 coshaL 


J 
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+ 


2Tj  a  sin  a(a  cosha+  T2  sinh  a+acosa+T2  sin  a)+  4a2  sin  a(asinh  a  +  T2  cosh  a)) 


(9.36) 


cosha  +  T2  sinh  a  +  a  cosa  +  T2  sin  a)  +2 
.  v2Tj  a  (cosh  a- cos  a)(a  sinh  a-T2  cosh  a  y 

+  2Tj  a  sinh  a(a  cosha+ T2  sinh  a+acosa+T2  sin  a)- 4a2  sinh  a(asin  a -T2  cosh  a)) 

(9.37) 

We  get- 

2Tja  (cosh  a  -  cos  a)(a  sinh  a  +  T2  cosh  a  +  a  sin  a  -  T2  cos  a) 

+  2Tja(acosha  +  T2  sinh  a  +  a  cos  a  +  T2  sin  a)(sin  a  -  sinh  a) 

+  4a2  (a  sinh  a.  sin  a  +  T,consha  sin  a  +  a  sinh  a  sin  a  -T,  sinh  a  cos  a)  =  0 

(9.38) 

2a[2asinh  a  sin  a  +  T2(cosha  sin  a  -  sinh  a  cos  a)] 

+  T,  (sin  a  -  sinh  a)[a(cosh  a  +  cos  a)  +  T2  (sinh  a  +  sin  a)] 

+  Tj  (cosh  a  -  cos  a)[a(sinh  a  +  sin  a)  +  T2  (cosh  a  -  cos  a  )]=0 

(9.39) 

2a[2a  sinh  a  sin  a  +  T9  (cosh  a  sin  a  -  sinh  a  cos  a)] 

-  T,  (sinh  a  -  sin  a)[a(cosh  a  +  cos  a)  +  T2  (sinh  a  +  sin  a)] 

+  T,  (cosh  a  -  cos  a)[a(sinh  a  +  sin  a)  +  T2  (cosh  a  -  cos  a)]  =  0 


T, 2 (cosh  a -cos  a) 


+ 


(9.40) 


4a2  sinh  a  sin  a  +  2T,  (cosh  a  sin  a  -  sinh  a  cos  a) 

+  Tja(sinh  acosha-sinh  acosa  +  sinacosha-cosasina 
-sinh  a  cosha  — sinh  a  cosa  +  sin  a  cosh  a  +  cos  a  sin  a) 

+  TjT2  [cosh2  a  +  cos2  a  -  2cos  acosh  a  -  sinh  2  a  +  sin 2  aj  =  0  (9.41) 

T,T2[(cosh2  a -sinh  2  a)+(cos2  a  +  sin 2  a)- 2cos  acosh  aj 
+  2T[  a[sin  a  cosh  a  -  cos  a  sinh  a] 

+  2T2a(sin  a  cosha-cosa  sinh  a)  +  4a2  sinh  a  sin  a  =  0 

(9.42) 

i.e„  2T,T2(1  -  cos  a  cosh  a)  +  2a(T]  +  T2  )(sin  a  cosh  a  -cos  a  sinh  a) 

+  4a2  sinh  a  sin  a  =  0 

(9.43) 

2a2  sinh  a  sin  a  +  a  (T,  +  T2  )(sin  a  cosh  a  -  cos  a  sinh  a)  +  T^T)  (l  -  cos  a  cosh  a) 

(9.44) 


=  c 


sin  az  - 


2a  sin  a  +  T,  (cosh  a  -  cos  a) 
2a  sinh  a  +  T,  (cosh  a  -  cos  a) 


sinh  az 


T,  /  J  ,  2a sin  a  +  T. (cosha-cosa) 

H — ^(cosha-cosazl  1 - - 4 

2a  2a  sinh  a  +  T,  (cosh  a  -  cos  a 


(9.45) 


=  c<  sin  az  - 


2a  sin  a  +  T,  (cosh  a  -  cos  a) 
2a  sinh  a  +  T,  (cosh  a  -  cos  a) 


sinh  az 


2a 


2a  sinh  a  +  T)  (cosh  a  -  cos  a)  -  2a  sin  a  -  T)  (cosh  a  -  cos  a) 
2a  sinh  a  +  Tj  (cosh  a  -  cos  a) 


.  (cosh  az  -  cos  az)} 


(9.46) 
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J  cm  r/7 _ 

2a  sin  a  +  T,  (cosh  a  -  cos  a) 

S  olll 

2a  sinh  a  +  T,  (cosh  a  -  cos  a) 

sinh  az 


2a(sinh  a  -  sin  a) 


2a  sinh  a  +  T,  (cosh  a  -  cos  a) 


(cosh  az  -  cos  az) 


(9.47) 


y(z)  =  c 


J  cm  n  7 _ 

2a  sin  a  +  Tj (cosh  a  -  cos  a) 

S  ulll  KAiZ-i 

2a  sinh  a  +  T,  (cosh  a  -  cos  a) 

sinh  az 


Tj  (sinh  a  -  sin  a) 


2a  sinh  a  +  (cosh  a  -  cos  a) 


(cos  az  -  cosh  a) 


(9.48) 


The  final  frequency  equation  is  - 

2(aL2 )  sin  aL  sinh  aL  +  (aL^T,  +  T,  )(sin  aL  cosh  aL  -  cos  aL  sinh  aL) 

+  T,T2(l-  cos  aL  cosh  aL)  =  0  (9.49) 


The  expression  for  mode  shape  is  as  follows- 
y(z)  =  A  sinh  a  +  B  cosh  az  +  C  sin  az  +  D  cos  az 

We  have  B+D=0  &  B  =  —  (A  +  C) 

2a 

.'.  y(z)  =  Asinh  (aL)z  +  Csin  (aL)z  +  B(cosh(aL)z  -  cosh(aL)z) 

T 

=  Asinh  (aL)z  +  C  sin  (aL)z  H — 4  (A  +  C)(cosh(aL)z  -  cos  h(aL)z)  (9.50) 

2a 


=  A 


sinh  aLz  H — -  (cosh  aLz  -  cos  aLz) 
2a 


+  C 


T 

sin  aLz  H — 4  (Cosh  aLz  -  cos  aLz) 
2a 


(9.51) 


-c 


Tj  (cosh  aL  -  cos  aL)+2aL  sin  aL 
Tj  (cosh  aL- cos aL)+2aLsinh  aL 


sinh  aLz  H — -  (cosh  aLz 
2a  V 


-cos  a 


(9.52) 


=C 


T 

sin  aLz+— (cosh  aLz  -cos  aLz )- 


2a 


2aL  sin  aL  +  T,  (cosh  aL  -  cos  aL) 
2aLsinh  aL+T,  (cosh aL -cos aL) 


T 

sinh  aLz  + — (cosh  aLz  -  cos  aLz ) 
2a 


(9.53) 


y(z)  =  C 


sin  aLZ  - 


2aL  sin  aL  +  Tj  (cosh  aL  -  cos  aL) 
2aL  sinh  aL  +  T,  (cosh  aL  -  cos  aL) 


T1  (sinh  aL  -  sin  aL) 


2aL  sinh  aL  +  T,  (cosh  aL  -  cos  aL) 


(cos  aLz  -  cosh  aLz) 


(9.54) 


Or  the  above  equation  can  also  be  written  in  form  of  1/T1- 


=C 


|  sin  aLz 


(cosh  aL-cosaL)+(l/T1)aLsin  aL 
(cosh  aL- cos aL)+(l/T1)2aL sinh  aL 


(cos aLz -cosh aLz)  >(9.55) 


We  can  easily  see  that  above  equation  reduces  to  the  following  for  1/T,  =>  go  i.e.. 
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y(z)  =  C  {sinh  aL-sin  aL -  (cos azL -cosh  aLz)] 

l  (cosh  aL- cos  aL)v  'J 


(9.56) 


By  applying  the  boundary  conditions  we  obtain  the  exact  frequency  equation  (9.49) 
and  the  mode  shape  expression  (9.55)  and  at  T=oc  reduces  to  equation  (9.56) 


We  neglect  the  effect  of  internal  pressure,  P  and  rotatory  inertia,  J  and 
consider  the  influence  of  rotational  restraint  parameter,  T  only  on  seismic  response. 


We  know 


c  = 


'p^Rm~Jco2 

El 


Since  P  and  J  are  equal  to  zero,  c=0  and  if  the  roots  of  differential  equation  are  a,  p 
If  a=p,  then  substituting  c=0  in  (9.16)  we  get- 


(9.57) 

9.5  Exact  Analysis  -Seismic  Response 

The  various  dimensional  parameters  of  U-shaped  bellows  are  taken  to  be  same 
as  those  that  have  been  considered  by  Morishita  et  al.  [40]- 
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Bellows  inner  diameter  (DO  -545.8mm,  convolution  height,  h=30mm;  convolution 
pitch,  p=25mm;  bellows  thickness,  tb=0. 563mm:  Bellows  mean  diameter, 
Dm=550.8rnm:  number  of  convolutions,  N=20:  pitch  diameter  of  bellows, 
Dp=575.8mm  and  Bellows  length  L  =  0.5  m  respectively. 

The  total  mass  mtot  of  bellows  per  unit  length  is  found  out  using  the  expression  given 
in  chapter  3- 

m  tot  =  440.76kg/m. 

From  (9.4)  we  get  - 

A.  =  7.0016a/co  (9.58) 

The  rotational  restraint  parameters  Ti  and  T2  are  varied  from  a  minimum  value  of  0.01 
to  a  maximum  of  10 10.  However,  it  is  assumed  that  there  are  equal  rotations  on  either 
ends  i.e.T=Ti=T2. 

The  maximum  response  displacement  of  any  arbitrary  section  of  bellow  to  a  given 
response  spectrum  is  given  by- 


y  ■ 


P 


Dr 


(2nfy  2*N*p 

Where  [1  is  participation  factor  and  is  assumed  to  be  equal  to  1.0 


y(z) 


(9.59) 


y(z)  is  maximum  deflection  at  the  central  convoluted  portion  of  the  bellows  and 
reduces  gradually  from  crown  to  root.  Therefore,  y(z)  is  taken  as  1.0. 
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Sa  is  the  spectral  gravitational  acceleration  at  frequency  ‘f  and  has  units  of  min/s2. 

The  bending  stress  induced  by  lateral  vibration  of  bellows  is  not  caused  by  lateral 
displacement  of  crowns  directly,  but  by  axial  displacement  between  adjacent  crowns, 
which  is  the  deflection  of  neutral  axis  of  bellows.  So  we  must  relate  the  axial  and 
lateral  displacements  in  beam  mode  vibration  of  bellows. 

9.6  Results  and  Discussion 

Table  9.1  gives  a  comparison  of  the  exact  frequencies  obtained  using  the 
bisection  method  for  a  single  bellows  that  are  elastically  restrained  at  both  ends  vis-a-vis 
to  the  fixed- fixed  end  results  presented  by  Morishita  et  al.  [5]. 

Table  9.1  Comparison  of  Frequencies  at  T=  oc 


Mod 

e# 

Morishita 

[5] 

f,  Hz 

Exact 

f,  Hz 

1 

8.6 

8.05 

It  is  seen  that  the  frequencies  obtained  by  considering  the  ends  as  elastically 
restrained  are  lower  than  that  was  found  out  by  Morishita  et  al.  Therefore,  it  can  be 
concluded  that  the  seismic  response  by  exact  method  is  accurate. 


Table  9.2  presents  the  first  mode  frequency  at  T=Ti=T2.  The  time  period,  t  in  sec  and 
the  spectral  acceleration,  Sa  values  are  obtained. 
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Table  9.2  Fundamental  Frequency  for  various  values 
of  restraint  parameter  T=0.01  to  oc _ 


T 

f  (Hz) 

t  (sec) 

Sa 

(mm/s2) 

0.01 

0.8023 

1.25 

79 

0.1 

0.8061 

1.24 

105 

1.0 

1.6826 

0.6 

473 

10 

1.7134 

0.58 

479 

102 

1.8101 

0.553 

500 

103 

1.8114 

0.55 

498 

104 

2.6169 

0.382 

763 

105 

3.4224 

0.292 

1513 

10® 

4.2279 

0.237 

1315 

107 

5.0328 

0.12 

2000 

108 

5.6366 

0.177 

2210 

109 

6.4427 

0.155 

3000 

1010 

7.2476 

0.138 

2658 

1011 

8.0525 

0.124 

2263 
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Fig  9.2  Time  History  Curve  of  Seismic  Response 


It  is  observed  that  for  N=1  and  at  T=0.01  the  time  period  found  out  is  0.49sec, 
while  at  T=oc,  the  time  period  is  0.12sec.  It  is  seen  that  as  rotational  stiffness 
increases  and  approaches  infinity,  the  time  period  of  vibration  decreases.  A  standard 
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time  history  graph  of  seismic  wave  and  its  acceleration  response  is  used  for  obtaining 


the  seismic  accelerations  (Sa)  at  various  frequencies  as  shown  in  Fig  9.1. 


Table  9.3  gives  the  maximum  displacement,  ymax 


and  bending  stress  ab  values  for 


varying  rotational  restraint  parameter  T. 


Table  9.3:  Response  of  Bellows  to  Lateral 


Seismic  excitation 


T 

ymax(mm) 

CTb(kg/mm2) 

0.01 

3.1 

21.8 

0.1 

4.1 

28.9 

1.0 

4.2 

29.7 

10 

4.1 

28.8 

102 

3.86 

27.1 

103 

3.84 

27.0 

104 

2.82 

19.8 

105 

3.27 

23.0 

106 

1.86 

13.1 

107 

2.0 

14.0 

108 

1.76 

12.4 

109 

1.83 

12.9 

1010 

1.28 

9.0 

1011 

0.88 

6.22 

The  seismic  gravitational  acceleration  obtained  at  T=oc  and  time  period  of  t=0.124sec 
is  2263mm/s2.  Substituting  this  value  in  equation  (9.59),  the  maximum  lateral 
displacement  is  obtained  as  ymax  =  0.88mm  compared  to  0.84mm  by  Morishita  et  al 


[5]  at  T=oc. 
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Fig  9.3:  Influence  of  Rotational  Restraint  on  Seismic  Response 


£ 

s 


Rotational  Restraint,  T 


The  exact  method  developed  is  to  calculate  the  dynamic  characteristics  and 
seismic  response  of  bellows  for  its  lateral  vibrations.  The  theoretical  formulations  are 
based  on  the  equations  of  an  equivalent  Euler-Bernoulli  beam,  and  the  bending  stress 
due  to  these  vibrations  is  estimated.  In  case  of  lateral  vibrations  it  is  found  that  the 
influence  of  elastically  rotational  restraint  at  either  ends  to  its  natural  frequencies  is 
significant. 

Also,  it  can  be  seen  from  fig  9.2  that  there  is  a  considerable  influence  of 
rotational  restraint  on  seismic  response.  As  the  rotational  restraint  parameter  T 
increases  from  0.1  to  1.0,  there  is  a  steep  rise  in  seismic  response  followed  by  a  fall  at 
T=102.  Beyond  T=102,  the  seismic  response  falls  up  to  T=oc,  except  at  106  where 
there  is  a  rise  because  of  sudden  peak  in  acceleration. 
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CHAPTER  10 

CONCLUSIONS 


10.1  Conclusions 

It  is  important  to  know  the  natural  frequencies  for  predicting  the  dynamic 
response  of  structures.  The  system  under  consideration  in  the  thesis  is  the  bellows 
expansion  joints.  Several  types  of  expansion  joints  are  used  in  practice.  However,  the 
thesis  considers  two  types  -single  and  double  or  Universal  type  of  expansion  joints 
only.  Theoretical  models  for  the  investigation  of  axial  and  transverse  vibrations  in 
single  and  double  bellows  were  developed  in  the  thesis.  The  Euler-Bernoulli  and 
Timoshenko  Beam  differential  equations  were  used  to  derive  the  “exact  frequency 
equations”  with  necessary  boundary  conditions.  It  is  seen  that  other  researchers  had 
considered  classical  fixed-fixed  end  conditions  only  in  determining  the  natural 
frequencies. 

However,  in  practice  it  is  observed  on  many  occasions  that  the  pipes  are 
subjected  to  equal/unequal  rotations.  Hence,  the  effects  of  variation  of  the  rotational 
restraint  parameter,  internal  pressure  and  flow  velocity  on  natural  frequencies  were 
studied  by  considering  elastically  restrained  ends. 

The  exact  solutions  were  then  compared  quantitatively  with  approximate 


solutions  of  other  investigators  to  verify  their  accuracy.  The  finite  element  method 
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was  used  in  establishing  the  closeness  of  values  with  exact  solutions  and  also 
experimental  investigations. 

In  order  to  simplify  the  double  bellows  expansion  joint  problem-  it  was 
divided  in  to  two  separate  problems  for  lateral  and  rocking  modes  respectively.  These 
two  problems  are  governed  by  the  same  differential  equation  but  with  two  different 
boundary  condition  sets. 

Comparison  of  these  calculated  natural  frequencies  with  experiments 
conducted  by  other  researchers  showed  very  close  agreement.  This  can  be  explained 
by  the  simplistic  classical  approach  considered  by  others  and  the  approximate 
methods  used  in  finding  out  the  natural  frequencies. 

The  analysis  provided  in  the  thesis  demonstrated  that  by  considering  the 
bellows  as  a  “beam”,  the  shear  influence  is  negligible  because  of  the  transverse 
flexibility  of  bellows.  On  the  other  hand  the  rotatory  inertia  of  bellows  cross-section 
including  the  fluid  trapped  is  an  important  factor  for  the  natural  frequencies  of 
bellows,  as  does  the  internal  pressure.  Similarly,  the  effect  of  the  variation  of  support 
mass  and  support  stiffness  on  natural  frequencies  is  studied  in  case  of  Universal  type 
of  expansion  joint. 

The  exact  solutions  were  within  10%  of  error  with  results  of  finite  element 


analysis  and  other  investigators. 
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The  dimensions  of  bellows  chosen  in  the  thesis  are  same  as  those  considered 
by  Jakubauskas.V.F,  Li  Ting-Xin  and  Morishita  et  al,  -only  to  arrive  at  logical 
comparison  purpose.  The  bellows  considered  here  were  of  single  ply  construction. 

10.2  Concluding  Remarks 

The  investigations  reported  in  literature  demonstrate  that  the  expressions 
currently  used  for  the  computation  of  axial  and  transverse  frequencies  in  bellows  are 
not  adequate  for  reasonable  predictions  in  design  practice.  It  was  found  that  the 
elastically  restrained  end  conditions  in  the  presence  of  internal  pressure  and  flow 
velocity  was  ignored  by  previous  authors. 

10.3  Future  Scope  of  Studies 

The  criteria  for  vibrations  of  multi-ply  bellows  goes  beyond  the  scope  of  this 
thesis.  Hence,  for  establishing  the  same  and  for  better  understanding  of  axial  and 
transverse  vibrations  of  single  and  double  multi-ply  bellows  -further  studies  are 
recommended.  Research  on  composite  material  bellows  is  also  interesting  and  worth 


attempting. 
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